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1.0  INTRODUCTION 


This  volume  of  the  report  discusses  the  computer  program  OPTFORCE  II 
from  the  User's  standpoint.  Volume  I,  Ref.  1,  discusses  the  theoretical 
development  of  force  method  optimization,  as  well  as  the  governing  equa¬ 
tions  and  solution  procedure  logic  incorporated  into  the  OPTFORCE  II  code. 

The  reader  must  refer  to  Ref.  1  for  in-depth  discussions  of  these  items. 

Included  in  the  OPTFORCE  II  program  is  the  capability  for  analyzing 
(static  and  dynamic)  and  weight  optimizing  truss  and  aerospace  type  struc¬ 
tures.  Both  isotropic  and  orthotropic  material  behavior  can  be  considered. 
Four  finite  elements  have  been  included;  namely,  truss  (rod),  symmetric 
shear  web,  triangular  and  quadrilateral  shaped  membranes.  Only  the  membrane 
elements  are  capable  of  handling  composite  structures  by  utilizing  a  layer¬ 
ing  procedure.  Section  II  discusses  the  aforementioned  items  and  in 
addition  an  overview  of  program  capabilities  and  options  is  presented.  The 
weight  optimization  procedure  is  also  briefly  described  in  that  section. 

Description  of  the  OPTFORCE  II  program  itself  is  given  in  Section  3.0. 
Particular  emphasis  has  been  given  to  those  sub-sections  which  describe  the 
input  and  output  data.  It  is  noted  here  that  all  input  data  cards  are  com¬ 
patible  with  NASTRAN  formats.  Each  required  input  card  is  discussed  thor¬ 
oughly  and  illustrated  in  the  sample  problem  solution.  Output  options 
are  described  and  typical  output  format  is  given  with  the  sample  problem. 

A  programmer's  manual  is  given  in  Appendix  A  of  this  volume.  It's 
purpose  is  to  fully  describe  program  logic  and  the  required  peripheral 
storage.  Individual  sub-routine  write-ups  are  presented.  A  short  des¬ 
cription  of  each  routine  is  included  to  aid  the  reader's  understanding  of 
the  program.  This  manual  is  specifically  written  for  the  computer  programmer, 
however. 
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The  rapid  reanalysis  development  has  resulted  in  a  computer  program 
described  in  Appendix  B.  Program  capabilities,  limitations,  input  and 
output  data  descriptions  and  illustrative  examples  are  given.  The  de¬ 
tailed  theoretical  development  of  the  rapid  reanalysis  procedure  is 
presented  in  Volume  I  of  this  report,  Reference  1. 
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2.0  OPTFORCE  II  OVERVIEW 


The  governing  equations  and  solution  procedure  logic  for  the  force- 
method  weight  optimization  code,  OPTFORCE  II,  is  developed  "in-toto"  in 
Volume  I,  Ref.  1,  of  this  report.  Because  of  the  vast  amount  of  technical 
material  presented  in  that  volume  only  a  brief  description  of  program 
capabilities,  finite  elements  and  the  weight  optimization  procedure  is 
given. 

The  main  objective  of  the  work  reported  herein  and  in  Volume  I  was 
to  develop  a  method  of  weight  optimization  using  the  force  method  of 
structural  analysis.  As  work  progressed  it  became  apparent  that  addi¬ 
tional  analysis  capabilities  other  than  weight  optimization  were  readily 
available  and  these  should  be  made  accessible  to  the  structural  analyst. 

As  a  result,  options  in  OPTFORCE  II  are  available  to  conduct  several 
types  of  analyses.  These  are  shown  in  Figure  1  in  flow  chart  form. 

Each  of  the  analysis  blocks  are  further  expanded  in  Figures  2  to  5. 
Examination  of  these  figures  shows  the  scope  of  analyses  available.  All 
of  these  are  resident  in  OPTFORCE  II  except  the  rapid  reanalysis  capability. 
This  analysis  tool  is  described  in  Appendix  B  of  this  report.  The  matrices 
noted  in  the  flow  charts  are  defined  in  Section  2.0  of  Volume  1,  however 
a  few  definitions  are  given  to  clarify  the  interpretation  of  the  flow 
charts. 

The  matrix  [a]  relates  the  external  grid  point  forces  (p)  to  internal 
stress  state  {s}  and  load  reactions  {R}.  [f]  is  a  normalized  diagonal 

flexibility  matrix  for  the  entire  structure;  [4>] ,  [ip]  &  [fi]  are  defined  as: 


m  - 

Cbx]T  [?]  [bx] 

(i) 

[ip]  - 

[b/  [f]  [D] 

(2) 

[n]  « 

&>  3T  [f]  Cd] 

(3) 
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The  superscripts  refer  to  other  flow  charts  for  more  information. 


Figure  1  General  Program  Outline 
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Figure  3  Basic  Analyses 
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Figure  4  Optimization 
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[b^]  &  [d]  arise  from  applying  a  structure  cutter  procedure;  where  [b^] 
represents  unit  values  of  self-equilibrating  force  systems  and  [d]  re¬ 
presents  forces  which  are  in  static  equilibrium  with  the  external  forces 

I 

{P} .  The  global  flexibility  matrix  is  defined  by  [f]  which  relates  {P} 
to  (A), {A}  being  the  displacement  of  the  structure,  [m]  is,  of  course, 
the  lumped  mass  matrix  of  the  structure.  The  set  of  redundants  is 
defined  by  the  vector  {X}  and  is  a  basic  unknown  in  the  force  method  as 
described  in  References  1  and  2. 

The  weight  optimization  procedure  coded  and  described  in  detail  in 
Volume  1,  Section  2.3.3  is  a  modification  and  expansion  of  the  one  developed 
in  Ref.  2.  Only  the  salient  features  of  the  procedure  are  given  here. 
Reference  to  Figure  4  will  aid  in  the  cursory  description  of  the 
solution  algorithm.  The  initial  guess  vector  for  structural  member  sizes 
(design  variables)  can  be  obtained  in  one  of  four  ways.  First,  the  initial 
sizes  may  be  obtained  directly  from  input  data  (OPDVIR  cards)  or  they  may 
be  set  to  the  minimum  sizes  A*.  Each  of  these  may  be  further  used  in  a 
stress  ratio  method.  These  four  options  are  controlled  through  input  codes  on 

the  OPTIM  card.  Selection  of  the  input  guess  vector  permits  entry  into 
the  linear  programming  (L.P.)  scheme  which  solves  for  sets  of  Lagrange 
multipliers  y.  The  particular  set  of  y's  solved  for  imply  satisfaction 
of  the  constraints  associated  with  the  non-zero  multipliers.  The  con¬ 
straints  considered  in  OPTFORCE  II  are:  minimum  sizes  of  design  variables, 
member  stress,  structural  displacement  and  natural  frequency  of  vibration. 

The  linear  programming  phase  is  cycled  three  times  as  shown  in  Figure  2.4 
to  obtain  a  new  design.  If  a  new  design  is  not  converged  upon  within 
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three  iterations  the  program  exits  L.P.  and  produces  a  full  stressed  design 
(FSD).  Note  that  the  frequency  constraint,  if  desired,  is  imposed  at  this 
time  and  the  linear  programming  phase  is  re-entered.  The  Lagrange  multi¬ 
pliers  are  now  checked  for  non-negativity  and  the  remaining  constraints 
are  checked  for  violations.  If  all  the  p's  are  positive  and  no  constraint 
is  violated,  the  minimum  weight  solution  is  found  as  the  FSD  or  the  design 
implied  by  the  LP  scheme  had  it  converged.  If  a  particular  p  is  negative 
then  the  corresponding  constraint  is  "turned-off"  and  the  multiplier  is 
set  equal  to  zero.  If  a  constraint  is  violated  the  corresponding  sizes  of 
design  variables  are  changed  to  negate  this  violation.  The  current  design 
at  this  point  in  the  solution  procedure  is  expressed  as  a  system  of  un¬ 
knowns  in  element  sizes,  structural  redundants  and  the  non-zero  Lagrange 
multipliers.  This  system  is  non-linear  in  nature  and  is  solved  using  the 
Newton-Raphson  procedure.  Output  from  the  optimization  procedure  is,  of 
course,  minimum  structural  weight,  member  stresses,  values  of  the  design 
variables  and  the  structures  displacement  behavior. 

Five  finite  elements  are  provided  in  OPTFORCE  II;  namely,  truss  or 
rod,  plane  stress  triangle,  plane  stress  quadrilateral  and  two  shear 
panels.  The  truss  elements  permit  analysis  of  truss  type  structures  like 
transmission  towers  for  example  and  more  importantly  portions  of  aerospace 
structures  such  as  lifting  surfaces.  The  symmetric  shear  panel  is  of 
direct  use  here  since  it  allows  modeling  of  only  one-half  of  the  lifting 
surface.  The  plane  stress  elements  are  membrane  type  elements  and  conse¬ 
quently  only  accept  in-plane  loadings.  Figure  6  displays  the  above 
elements  and  their  characteristics.  Notable  among  these  are  the  design 
variables  (denoted  as  the  vector  f A) ) . 

Note  that  the  membrane  elements  can  be  used  to  model  composite  materials. 


Element  characteristics:  Two  grid  points,  length  &  design  variable 

cross-sectional  area  A.  Isotropic  materials 
only.  Stress  output  Ox. 

a)  Truss  (rod)  Axial  Force  Member  ~  ROD 


Element  Characteristics:  Three  grid  points,  surface  area  As,  design 

variable  plate  thickness  t,  angle  of  ortho- 
tropy  0.  Orthotropic  and  isotropic  materials 
Centroidal  stress  output  oxm,  o ym,  Txym. 

b)  Triangle  Membrane  Plate  -  TRMEM 


Figure  6  Finite  Element  Library 


Element  Characteristics:  Four  grid  points,  surface  area  As,  design 

variable  plate  thickness  t,  angle  of  orthotropy 
0.  Orthotropic  and  isotropic  materials.  Cen- 
troidal  stress  output  0xm,  aym,  Txym* 

c)  Quadrilaterial  Membrane  Plate  ~  QDMEM1 


Element  Characteristics:  Two 
grid  points,  surface  area  As, 
design  variable  web  thickness 
t.  Isotropic  materials  only. 
Centroidal  stress  output  Txy. 
N.B.  always  define  grid  points 
in  positive  global  Z  direction. 


d)  Symmetric  Shear  Panel  ~  WEB 


Figure  6  (cont'd). 


Element  Characteristics:  Four  grid  points,  surface  area  As, 

design  variable  web  thickness  t.  Iso¬ 
tropic  materials  only.  Centroidal  stress 
output  Txy.  N.B.  Always  define  grid  points 
in  positive  global  Z  direction. 

e)  Shear  Web  -  SHEAR 


Figure  6  (cont'd.) 


This  was  accomplished  by  providing  an  element  layering  capability  coupled 
with  orthotropic  material  behavior.  Each  layer  of  the  plate  is  idealized 
as  a  separate  membrane  plate  finite  element.  When  modeling  a  composite 
plate  consisting  of  many  layers,  the  gridpoints  of  all  the  membrane  elements 
or  layers  of  the  plate  are  selected  to  be  the  same.  This  restriction  ensures 
that  each  of  the  layers  of  the  composite  plate  are  subjected  to  the  same 
deformation  state  and  is  in  conformance  with  small  displacement  plate 
theory  in  which  it  is  assumed  that  straight  normals  to  the  plate  median 
surface  before  deformation  remain  straight  after  deformation.  To  ensure 
that  this  condition  is  satisfied,  the  layers  of  the  composite  plate  must 
be  oriented  symmetrically  with  respect  to  the  median  surface  of  the  com¬ 
posite  plate. 

Each  layer  in  a  composite  membrane  plate  was  assumed  to  be  composed 
of  fibers  embedded  in  a  matrix  material.  The  fibers  are  assumed  to  be 
oriented  in  such  a  way  as  to  result  in  their  characterization  as  an  ortho¬ 
tropic  material.  A  typical  orthotropic  layer  of  this  type  is  illustrated 

in  Figure  7.  The  coordinate  system  (x  ,y  )  are  the  material  axes  of 

m  in 

orthotropy  in  which  the  x  axis  is  oriented  parallel  to  the  fibers.  The 

m 

angle  9  gives  the  orientation  of  the  material  axis  with  respect  to  the 

local  element  x  axis.  This  axis  is  defined  by  the  side  of  the  element 
e 

connecting  grid  points  1  and  2.  Note  that  0  is  measured  in  a  counterclock¬ 
wise  direction,  from  the  x  axis  to  the  x  axis. 

e  m 

The  relationship  between  the  elastic  stresses  and  strains  for  an 
orthotropic  material  in  plane  stress  referenced  to  the  material  axis  of 
orthotropy  shown  in  Figure  7  is  given  by: 
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b)  Tyoical  Orientation  of  Membrane  Element  in 
Fiber-reinforced  Composite  Layer 


Figure  7  Fiber-reinforced  Composite  Idealization 
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These  stress-strain  relations  are  used  in  the  calculation  of  membrane 
finite  element  characteristics  and  furthermore  are  reducible  to  isotropic 
material  behavior.  As  a  result  it  is  possible  to  construct  a  composite 
plate  composed  of  both  fiber-reinforced  orthotropic  materials  and  isotropic 
metallic-type  materials.  It  is  noted  here  that  the  MAT1  and  MAT2  input 
data  cards  are  used  to  define  the  above  material  properties. 

The  Mises-Henchy  stress  failure  criteria  has  been  inserted  into  the 
OPTFORCE  II  program  as  follows: 

2.  k 

(5) 


2  2  2  k 

Y  =a.  (a  +0  -a  a  +3x  V 
i  i  x  y  x  y  xy 


where  a  .a  and  t  are  element  stresses  calculated  in  the  local  element 
x  y  xy 

x  -y  axis  system.  A  is  the  element's  design  variable;  for  a  bar  element, 
m  in  I 

for  example,  A^  would  be  its  cross-sectional  area.  The  above  criteria  was 


used  to  formulate  the  stress  constraint  expression: 

-•  1  <o 


Vi* 


(6) 


where  is  the  yield  stress  or  some  other  failure  stress  value  for  the 

i—  finite  element.  It  is  a  User  input  and  is  given  on  the  MAT1,  MAT2 

Yi 

cards.  The  quantity  —  is  output  with  the  local  stresses  for  each  element. 

Ai 
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3.0  0PTF0RCE  II  COMPUTER  PROGRAM  DESCRIPTION 
3.1  Input  Overview 


This  section  describes  the  input  data  required  to  execute  the  NASTRAN 
compatible  optimization  program  OPTFORCE  II.  This  input  consists  of  a 
deck  beginning  with  "BEGIN  BULK"  and  ending  with  "ENDDATA".  All  data 
cards  are  optional  except  those  defining  controls,  loads,  grid,  and  boun¬ 
dary  conditions.  Elements  must  also  be  defined. 


The  NASTRAN  input  feature  of  submitting  data  in  any  order,  either 
left  or  right  adjusted,  has  been  preserved.  The  GRID  cards,  SPC  cards, 

FORCE,  and  KAT1 ,  MAT 2  cards  are  similar  to  NASTRAN  input.  The  following 
is  a  summary  of  OPTFORCE  II  Input  Cards: 

BEGIN  Bulk  First  card  of  input 

CONROD  Axial  force  member  property  and  connection 

CQDMEM1  Quadrilateral  membrane  connection 

CROD  Axial  force  element  connection 

CS HEAR  Shear  panel  element  connection 

CTRMEM  Triangular  membrane  element  connection 

CWF.B  Shear  web  element  connection 

FORCE  Static  load  at  grid  point 

GRID  Grid  point  coordinates 

♦ICON  Displacement  constraints 

MAT1  Isotropic  material  properties 

MAT2  Anisotropic  material  properties 

*OPDVIR  Selected  element  design  variables 

OPLOADS  Loads  for  optimization 

OPTIM  Optimization  control  parameters 

P0DMEM1  Property  card  for  quadrilateral  membrane 

PROD  Property  card  for  axial  force  member 

PSHEAR  Property  card  for  shear  panel 

PTRMEM  Property  card  for  triangular  membrane 

PWEB  Property  card  for  shear  web 

SPC  Single  point  constraint 

SPC1  Sets  of  single  point  constraints 

TITLE  Title  card  Information 

ENDDATA  End  of  data  deck 

3.2  Description  of  Input  Data 

In  this  section,  detailed  OPTFORCE  II  input  data  card  descriptions 
are  presented.  Each  data  card  is  described  individually  in  alphabetical  order 


♦These  cards  are  optional 
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Input  Data  Card  BEGIN  Bulk 
Deacriptlon;  First  Card  of  Input 

123456789 

BEGIN  Bulk 
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Input  Data  Card  C0NR0D  Axial  Force  Element  Property  and  Connection 

Description:  Defines  an  axial  force  element  (R0D)  without  reference  to 
a  property  card. 


Format  and  Example: 


1 

C0NR0D 

C0NR0D 

Field 

EID 

Gl,  G2 

MID 

A 


2 

3 

4 

5 

6  7  8 

9  10 

EID 

Gl 

G2 

MID 

A 

2 

16 

17 

23 

2.69 

Contents 

Unique  element  identification  number  (Integer  >0) 

Grid  point  identification  numbers  of  connection  points 
(Integer  >0;  Gl  ^  G2) 

Material  identification  number  (Integer  >0) 

Area  of  member  (Real) 


Remarks :  1.  Element  identification  numbers  must  be  unique  with  respect 

to  all  other  element  identification  numbers. 


2.  C0NR0D  cards  may  only  reference  MAT1  material  cards. 


3.  A  is  used  as  the  design  variable  minimum  for  optimization 
runs,  or  as  the  design  variable  for  statics  ordynamic  runs. 
If  an  OPDVIR  card  group  is  included  in  the  bulk  data  deck 
then  A  is  ignored. 

4.  See  CR0D  for  alternative  method  of  rod  definition. 
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Input  Data  Card  CQDMEM1  Quadrilateral  Element  Connection 


Description:  Defines  a  quadrilateral  membrane  element  (QDMEM1), 


Format  and  Example: 


1  2  3  4  5  6  7 

CQDMEM1  EID  PID  G1  G2  G3  G4 

CQDMEM1  72  10  13  14  15  16 


Field 


Contents 


EID  Element  identification  number  (Integer  >0) 

PID  Identification  number  of  a  PQDMEM1  property  card  (Integer  >0) 

Gl,  G2,  Grid  point  identification  numbers  of  connection  points 

G3,  G4  (Integer  >0);  Gl  +  G2  t  G3  *  C4) 

TH  Material  property  orientation  angle  in  degrees  (Real). 

The  sketch  below  gives  the  sign  convention  for  TH. 


Remarks :  1.  Element  identification  numbers  must  be  unique  with  respect 

to  all  other  element  identification  numbers. 

2.  Gridpoints  Gl  through  G4  must  be  ordered  consecutively 
around  the  perimeter  of  the  element. 

3.  All  interior  angles  must  be  less  than  180  degrees. 
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Input  Data  Card  CR0D  Axial  Force  Element  Connection 


Description:  Defines  an  axial  force  element  (R0D)  with  reference  to  a 
property  card. 


Format  and  Example: 


12  3  4 

CR0D  EID  PID  Cl 

CR0D  12  13  21 


5  6  7  8 

G2  EID  PID  G1 

23  3  12  24 


9  10 

G2 
5 


Field 

EID 

PID 

Gl,  C.2 


Contents 

Element  identification  number  (Integer  >0) 

Identification  number  of  a  PR0D  property  card 

Grid  point  identification  numbers  of  connection  points 
(Integer  >0;  Gl  +  G2) 


Remarks : 


1.  Element  identification  numbers  must  be  unique  with 
respect  to  all  other  element  identification  numbers. 

2.  See  C0NR0D  for  alternative  method  of  rod  definition 

3.  One  or  two  R0D  elements  may  be  defined  on  a  single  card. 
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Input  Data  Card  CSHEAR 


Shear  Panel  Element  Connection 


Description; 


Defines  a  shear  panel 


element  (SHEAR) . 


Format  and  Example: 


1  2  3  4  5 

CSHEAR  EID  PID  G1  G2 
CSHEAR  3615 


6 

G3 

3 


7 

G4 

7 


8  9  10 


Field 

EID 

PID 

Gl,  G2 
G3,  G4 


Contents 

Element  identification  number  (Integer  >0) 

Identification  number  of  a  PSHEAR  property  card 

Grid  point  identification  numbers  of  connection  points 
(Integer  >0;  Gl  #  G2  t  G3  *  G4) 


Remarks :  1.  Element  identification  numbers  must  be  unique  with 

respect  to  all  other  element  identification  numbers. 

2.  Grid  points  Gl  thru  G4  must  be  ordered  consecutively 
around  the  perimeter  of  the  element. 
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Input  Data  Card  CTRMEM  Triangular  Element  Connection 


Description:  Defines  a  triangular  membrane  element  (TRMEM) . 


Format  and  Example: 


12  3 

CTRMEM  EID  PID 

CTRMEM  16  2 


4  5  6 

G1  02  G3 

12  1  3 


7  8  9  10 

TH 
16.3 


Field 

EID 

PID 


Contents 

Element  identification  number  (Integer  >0) 
Identification  number  of  a  PTRMEM  property  card 


Gl,  C2,  Gridpoint  identification  numbers  of  connection  points 

G3  (Integer  >0;  Gl  ^  G2  #  G3) 


TH 


Material  property  orientation  angle  in  degrees  (Real). 
The  sketch  below  gives  the  sign  convention  for  TH. 


Remarks : 


1. 


Element  identification  numbers  must  be  unique  with 
respect  to  all  other  element  identification  numbers. 
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Input  Data  Card  CWEB  Shear  Web  Element  Connection 

Description:  Defines  a  2  node  symmetric  web  element  (Web)  of  the 
structural  model. 


Format 

and  Example: 

1 

2 

3  4 

5  6  7 

8  9  10 

CWEB 

EID 

PID  Gl 

G2 

CWEB 

4 

2  1 

3 

Field 

Contents 

EID 

Element  identification 

number 

PID 

Identification  of  PWEB 

property  card 

Gl,  G2 

Gridpoint 

identification  numbers 

Remarks: 


1.  See  Figure  2.6 


Input  Data  Card  F0RCE  Static  Load 


Description:  Defines  a  static  load  at  a  grid  point  by  specifying  a 
vector. 


Format  and  Example: 

1  2  3  4 

F0RCE  SID  C 

F0RCE  2  5 


5  6 

F  N1 

2.9  0.0 


7  8  9  10 

N2  N3 

1.0  0.0 


Field 

SID 

C 

F 

Nl,  N2 
N3 

Remarks : 


Contents 

Load  set  identification  number  (Integer  >0) 

Gridpoint  identification  number  (Integer  >0) 

Scale  factor  (Real) 

Components  of  Vector  (Real) 

1.  The  static  load  applied  to  gridpoint  G  is  given  by 

f  =  F  N 

where  N  is  the  vector  defined  in  fields  6,  7,  and  8. 

2.  Load  set  is  selected  on  the  OPLOADS  card. 
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Input  Data  Card  GRID  Grid  Point 


Description:  Defines  the  location  of  a  geometric  gridpoint  of  the 

structural  model  and  its  permanent  single-point  constraints. 


Format  and  Example: 


1  234  5  6789  10 

GRID  ID  XI  X2  X3  PS 

GRID  2  1.0  2.0  3.0  13 


Field  Contents 

ID  Gridpoint  identification  number  (0<  Integer<  999999) 

XI,  X2  Location  of  the  grid  point 

X3 

PS  Permanent  single-point  constraints  associated  with  grid- 

point  (any  of  the  digits  1-3  with  no  imbedded  blanks) 
(Integer  £0  or  blank) 


Remarks:  1.  The  coordinate  system  defined  on  all  GRID  cards  is 

called  the  Global  Coordinate  System.  All  degress-of- 
freedom,  constraints,  and  solution  vectors  are  expressed 
in  the  Global  Coordinate  System. 
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Input  Data  Card  ICON  Displacement  Constraints 


Description:  Defines  the  displacement  constraint  limits  (both  lower  and 

upper)  and  grldpoint  components  which  require  constraints. 


Format  and 

Example : 

1 

2 

3  4 

5  6 

7 

8  9  10 

ICON  Components 

Lower  Upper 

Cl  G2 

G3 

G4  G5 

ICON 

13 

1.0  3.0 

1  2 

3 

6 

Field 

Contents 

Components 

1-3  signifies  component 

numbers 

Lower 

Lower 

limit 

Upper 

Upper 

limit 

Gl,  G2 

Gridpoint  numbers.  Up 

to  5  gridpoint  numbers  may  appear 

C3,  etc. 

on  one 

ICON  card. 

Remarks : 

1.  Each  ICON  card  contains  uo  to 

5 

gridpoints.  When  more 

than  5  are  desired. 

list  them 

on 

separate  ICON  cards. 

2.  For  example,  let  there  be  the  following  individual 
constraints:  Ul,  U2,  Wl,  W2  greater  than  1.0  and  less 
than  4.0.  VI,  V2,  V3  greater  than  1.0  and  less  than 
8.0. 

The  required  ICON  cards  would  be: 

1  2  3  4  5  6  7 

ICON  13  1.0  4.0  1  2 

ICON  2  1.0  8.0  1  2  3 


Input  Data  Card  MAT1  Material  Property  Definition 


Description:  Defines  the  material  properties  for  linear,  temperature- 
independent,  isotropic  materials. 


Format  and  Example :  (Consists  of  two  cards) 


1 

2 

3  4 

5 

6  ; 

7  8  9  10 

MAT1 

MID 

£  G 

NU 

RH0 

+abc 

MAT1 

17 

3. +7  1.9+7 

4.28 

ABC 

+abc 

SL 

SU 

+ABC 

20.  +4 

15. +4 

Field 

MID 

E 

G 

NU 

RHO 

SL.SU 


Contents 

Material  identification  number  (Integer  >0) 

Young's  modulus  (Real  £0.0  or  blank) 

Shear  modulus  (Real  £0.0  or  blank) 

Poisson's  ratio  (-1.0  <Real  10.5  or  blank) 

Material  density  (Real)  (lbs/cu.  in.)  or  (Ibs/cu.  ft.) 
Lower  stress  limit,  upper  stress  limit 


Remarks : 


1.  One  of  E  or  G  must  be  positive  (i.e.  either  E  >0.0 
or  G  >0.0  or  both  E  and  G  may  be  >0.0). 

2.  If  any  one  of  E,  G  or  NU  is  blank,  it  will  be  computed 
to  satisfy  the  identify  E  =  2(1+NU)G;  otherwise,  values 
supplied  by  the  user  will  be  used. 

3.  The  material  identification  number  must  be  unique  for 
all  MAT1  and  MAT 2  cards. 

4.  SU  is  used  as  the  upper  limit  for  the  Mises-Hencky 
criteria  of  the  element  referencing  this  material  during 
an  optimization  analysis. 
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Input  Data  Card  MAT 2  Material  Property  Definition 


Description:  Defines  the  material  properties  for  linear,  anisotropic 
materials . 


Format 

and 

Example:  (consists  of  2  cards) 

1 

2 

3  4  5  6  7 

8 

9 

10 

MAT  2 

MID 

Gil  G12  G13  G22  G23 

G33 

RHO 

+abc 

MAT  2 

13 

6.2+3  6.2+3 

5.1+3 

0.056 

BC 

+abc 

SL 

SU 

+BC 

20. +5  15. +3 

Field 

Contents 

MID 

Material  identification  number 

(Integer 

>0) 

Gij 

The  material  property  matrix  (Real) 

RHO 

Material  density  (Real)  (lbs/cu 

.in.)  or 

(lbs/cu. 

ft.) 

SL,  SU 

Lower  stress  limit,  upper  stress  limit 

Remarks 

: 

1.  The  material  identification 

numbers 

must  be 

unique 

for  ail  MAT1,  MAT2  cards. 


2.  The  convention  for  the  Gij  in  fields  3  through  8  is 
represented  by  the  matrix  relationship. 


3.  Only  TRMEM  and  CQDMEM  elements  may  use  MAT2  cards. 
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Input  Data  Card  OPDVIR 


Description:  Defines  design  variables  for  selected  elements.  These  design 
variables  are  used  as  starting  guesses  in  an  optimization  run 
or  else  as  specified  values  in  a  statics  or  dynamics  run. 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

OPDVIR 

El 

DVAk 
.1  1 

E2 

DVAR. 

E3 

DVAR. 

E4 

DVAR. 

OPDVIR 

1 

3 

.6  2 

4 

.50  3 

7 

.234 

Field  Contents 


El,  E2,  Element  numbers  for  which  design  variables  are  input, 

etc . 


DAVR^ ,  DVAR^,  Design  variable  starting  guesses  for  elements  1,  2,  etc. 
etc . 


Remarks :  1.  Up  to  4  elements  may  be  defined  on  one  OPDVIR  card. 

More  elements  may  be  defined  by  successive  OPDVIR 
cards . 

2.  The  use  of  these  values  is  dependent  on  the  control 
option  selected  on  the  OPTIM  input  card  -  field  #9, 
OPT. 


a.  If  OPT  =  GOPT  the  program  will  opt lmize 

using  the  design  variable  minimums  input  in  this 
section  as  a  starting  guess. 

b.  If  OPT  =  GSTA,  the  program  will  perform  a  statics 
solution  using  the  design  variables  input  in  this 
section. 

c.  If  OPT  ■  GDYN,  the  program  will  perform  a  dynamics 
solution  using  the  design  variables  input  in  this 
section. 

d.  If  OPT  -  GOPD,  the  program  will  optimize  with 
dynamic  constraints  using  the  design  variables 
input  in  this  section. 


Input  Data  Card  OPLOADS 


Description:  Defines  the  character  of  the  loads  used  in  the  optimization 
procedure . 


Format  and  Example: 


1 

OPLOADS 

OPLOADS 

2 

LOAD 

1 

3 

4  5 

Fi 

5 

b  7 

8 

9 

Field 

Contents 

LOAD 

LOAD 

number 

identification  (any 

number) 

F. 

SID 

number 

on  FORCF  eard(s) 

1 


Remarks :  1.  Da  eh  FORCE  card  defining  this  external  load  must  contain 

the  SID  number  entered  for  F.. 

1 

2.  Each  external  load  requires  an  OPLOADS  entry. 
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Input  Data  Card  OPTIM 


Description:  Defines  various  control  parameters  for  the  OPTFORCE  II 

program. 


Format  and  Example : 


1 

2 

3 

5 

6 

7 

8 

9 

10 

OPTIM 

PRI 

PRS 

ISTRES 

NMODE 

MAP 

CONVD 

WS 

OPT 

OPTIM 

No 

Yes 

Yes 

6 

50 

.000  5 

495.0 

OPT 

Field  Contents 


PRI  No/Yes 

PRS  No/Yes 

ISTRES  No /Yes 

NMODE 

MAX 

CONVD 

WS 

OPT  OPT 

COPT 
STA 
CSTA 
DYN 
GDI  N 
OPTD 

GOPD 


Print  intermediate  debug  information  (file  101) 

Pi  ini  stresses  and  displacements  every  iteration 

Use  stress-ratio  method  for  initial  guess  when  optimizing 

Number  of  modes  in  dynamic  analysis 

Maximum  number  of  iterations  permitted  (default=50) 

Convergence  criteria  (def ault= .0001 ) 

Frequency  constraint  limit  (in  Hertz) 

Optimize  using  design  variable  on  property  card 
Optimize  using  OPDVIR  starting  guess 

Statics  analysis  using  design  variable  on  property  card 

Statics  analysis  using  OPDVIR  starting  guess 

Dynamics  analysis  using  design  variable  on  property  card 

Dynamics  analysis  using  OPDVIR  starting  guess 

Optimize  with  dynamic  constraint  using  design  variable  on 
property  card 

Optimize  with  dynamic  constraint  using  OPDVIR  starting 
guess 


Remarks: 


CONVD  -  Convergence  Criteria  (design  variables,  eigen 
vectors,  redundants  and  Lagrange  multipliers).  If  two 
successive  iterations  of  all  variables  meet  the  criteria 
.  .  I  then  the  iteration  process  is  terminated. 

£  CONVD 
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Input  Data  Card  PQDMEM1  Quadrilateral  Membrane  Property 


Description:  Defines  the  properties  of  the  quadrilateral  membrane  as 
referenced  by  the  CQDMEMl  card. 


Format  and  Example : 


1 

PQDMEM1 

PQDMEM1 


2 

PID 

235 


3  456  789  10 

MID  T  PID  MID  T 

2  0.5 


Field 

PID 

MID 

T 


Contents 

Property  identification  number  (Integer  >0) 
Material  identification  number  (Integer  >0) 
Minimum  thickness  of  membrane  (Real  >0.0) 


Remarks :  1.  All  PQDMEM1  cards  must  have  unique  property  identifi¬ 

cation  numbers. 

2.  One  or  two  quadrilateral  membrane  properties  may  be 
defined  on  a  single  card. 
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Input  Data  Card  PR0D  Rod  Property 


Description:  Defines  the  properties  of  a  rod  as  referenced  by  the  CR0D 
card. 


Format  and  Example: 


12  3 

PR0D  PID  MID 

PR0D  17  23 


4  5  6  7  8  9  10 

A 

42.6 


Field 

PID 

MID 

A 


Contents 

Property  identification  number  (Integer  >0) 
Material  identification  number  (Integer  >0) 
Area  of  rod  (Real) 


Remarks :  1.  PR0D  cards  must  all  have  unique  property  identification 

numbers . 

2.  PR0D  cards  may  only  reference  MAT1  material  cards. 

3.  For  remarks  on  use  of  A,  see  CONROD  card  description. 
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Input  Data  Card  PSHEAR  Shear  Panel  Property 

Description :  Defines  the  elastic  properties  of  a  shear  panel  as 
referenced  by  the  CSHEAR  card. 


Format  and  Example : 


1 

2 

3 

4 

5  6 

7  8  9 

10 

PSHEAR 

PID 

MID 

T 

PID 

MID  T 

PSHEAR 

4.9 

14 

6  4.9 

Field 

Contents 

PID 

Property 

identif icat ion 

number  (Integer 

>0) 

MID 

Material 

identif icat ion 

number  (Integer 

>0) 

T 

Minimum 

thickness  of  shear  panel  (Real  ^ 

0.0) 

,  Remarks :  1.  All  PSHEAR  cards  must  have  unique  identification  numbers. 

2.  PSHEAR  cards  may  only  reference  MAT1  material  cards.  i 

i 

3.  One  or  two  shear  panel  properties  may  be  defined  on 

|  a  single  card. 

! 


Input  Data  Card  PTRMEM  Triangular  Membrane  Property 

Description:  Defines  the  properties  of  a  triangular  membrane  element 
as  referenced  by  the  CTRMEM  card. 


Format  and  Example: 


1  2  3  456  789  10 

PTRMEM  PID  MID  T  PID  MID  T 

PTRMEM  17  23  4.25 


Field  Contents 

PID  Propertv  identification  number  (Integer  >0) 

MID  Material  identification  number  (Integer  >0) 

T  Membrane  thickness  (Real  >0.0) 

Remarks:  1.  All  PTRMEM  cards  must  have  unique  property  identi¬ 

fication  numbers. 

2.  One  or  two  triangular  membrane  properties  may  be 
defined  on  a  single  card. 


35 


Input  Data  Card  PWEB  WEB  Property 


Description:  Define  the  properties  of  2  node  WEB  as  referenced  by  the 
CWEB  card. 


Format  and  Example: 


1 

PWEB 

PWEB 


2  3  4  5 

PID  MID  T 

235  2  0.5 


6  7  8  9 

PID  MID  T 


10 


Field 

PID 

MID 

T 


Contents 

Property  identification  number  (Integer  >0) 
Material  identification  number  (Integer  >0) 
Thickness  of  the  shear  web  (Real  >0.0) 


Remarks :  1.  All  PWEB  cards  must  have  unique  property  identification 

numbers . 

2.  One  or  two  WEB  properties  may  be  defined  on  a  single 
card. 


Input  Data  Card  SPC  Single-Point  Constraint 
Description:  Defines  sets  of  single-point  constraints 


Format  and  Example: 


1 

2 

3 

4 

5  6 

7 

8  9 

SPC 

SID 

G 

C 

G 

C 

SPC 

2 

32 

231 

5 

1 

Field 

Contents 

SID  Identification  number  of  single-point  constraint  set 

(Integer  >0) 


G  Grid  point  identification  number  (Integer  >0) 

C  Component  number  (any  unique  combination  of  the  digits 

1-3  (with  no  imbedded  blanks)  when  point  identification 
numbers  are  gridpoints) . 


Remarks :  1.  Single-point  constraint  sets  must  be  present  in  the 

input . 

2.  From  one  to  twelve  single-point  constraints  may  be 
defined  on  a  single  card. 

3.  SPC  degrees  of  freedom  may  be  redundantly  specified 
as  permanent  constraints  on  the  GRID  card. 

4.  The  SID  number  should  be  the  same  on  all  SPC  cards. 
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Input  Data  Card  SPCI  Single-Point  Constraint 


Description:  Defines  sets  of  single-point  constraints 


Format  and  Example: 


1 

2 

3 

4 

5 

6  7 

8  9 

SPCI 

SID 

C 

Gl 

G2 

G3  G4 

G5  G6 

SPCI 

3 

2 

1 

3 

10  9 

6  5 

Alternate  form: 

SPCI 

SID 

C 

GID1 

"THRU" 

GID2 

SPCI 

313 

12 

6 

THRU 

32 

Fi^ld 

Contents 

SID 

L 

dent  if icat ion 

number 

of  single-point 

constraint 

(Integer  >0) 

C  Component  number  (any  unique  combination  of  the  digits 

1-3  with  no  imbedded  blanks)  when  point  identification 
numbers  are  gridpoints. 


Gl,  GID1  Grid  point  identification  numbers  (Integer  >0) 


Remarks : 


1.  Single-point  constraint  sets  must  be  present  in  the 
input . 

2.  SPC  degrees  of  freedom  may  be  redundantly  specified 
as  permanent  constraints  on  the  GRID  card. 

3.  All  gridpoints  referenced  by  GID1  thru  GID2  roust 
exist . 

A.  No  NASTRAN  "continuation"  cards  are  allowed. 


38 


Input  Data  Card  TITLE  Title  Card 


Description:  Information  to  be  printed  at  the  beginning  of  the  computer 

listing . 


Format  and  Example: 

123456789  10 

TITLE 

Plate  Membrane  Prob .  No.  2 


Input  Data  Card  ENDDATA 


Description:  Defines  the  end  of  the  data  deck. 


Format  and  Example: 

123456789  10 

ENDDATA 

ENDDATA 


Remarks :  i.  This  card  required  even  if  no  physical  data  cards 

exist  in  the  deck. 

2.  ENDDATA  must  begin  in  columns  1  or  2. 

3.  Failure  to  include  this  card  will  result  in  an 
operating  system  termination  caused  by  input  end 
of  file  error. 
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3.3  Description  of  Output  Data 

The  format  of  the  output  was  designed  with  the  engineering  User  in 
mind.  Liberal  use  of  labels  was  emphasized  which  serve  to  identify  the 
output  data  and  guide  the  User  through  the  optimization  solution  procedure. 
Output  print  options  are  limited  to  those  available  on  the  OPTIM  input 
card.  This  input  permits  printing  of  output  for  "bebug"  purposes.  The 
User  should  use  this  option  when  solutions  are  troublesome.  A  second  op¬ 
tion  permits  the  printing  of  stresses  and  displacements  at  every  iteration 
in  the  optimization  process. 

Typical  output  is  given  in  Section  3. A.  An  echo  print  of  the  input 
data  is  always  presented  first.  This  is  followed  by  output  pertinent  to 
the  analysis  being  conducted  such  as  shown  in  Section  3.4  for  optimization 
solutions.  Execution  of  static  or  dynamic  analyses  yields  standard  struc¬ 
tural  analysis  data,  ie.  grid  point  displacements,  element  stresses  and 
structural  reactions.  Vibration  mode  shapes  are  also  given  as  required. 
These  same  data  are  given  at  the  end  of  an  optimization  analysis  in  con¬ 
junction  with  the  minimum  weight  value  and  associated  design  variable 
vector. 
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3.4  Illustrative  Example 

A  seventeen  bar  cantilever  truss  structure  is  used  to  illustrate  the 
input/output  features  of  OPTFORCE  II.  Table  1  displays  material  properties, 
minimum  sizes  of  design  variables,  allowable  stress  levels  and  loading  con¬ 
dition.  Figure  8  shows  the  geometric  layout  and  finite  element  modeling  of 
the  truss.  Input  data,  as  written  in  NASTRAN  format,  is  given  in  Figure  9. 
Each  input  card  is  described  in  detail  in  Section  3.2  of  this  report.  The 
input  data  deck  is  very  similar  to  that  used  in  conducting  any  structural 
analysis  using  the  NASTRAN  computer  code.  In  fact,  the  only  new  data  In 
this  application  of  the  OPTFORCE  II  code  are  the  OPLOA.DS,  FORCE  and  OPTIM 
cards.  The  OPTIM  card  is  particularly  important  since  it  controls  analysis 
type,  convergence  criteria  and  output  print  options. 

Problem  output  is  given  in  Figure  10.  The  User's  input  data  deck  is 
reproduced  for  reference  purposes.  Selected  output  from  the  solution  pro¬ 
cedure  is  given  next.  Output  statements  and  the  use  of  Figure  4  should  be 
sufficient  to  guide  the  User  through  the  optimization  procedure.  When  the 
final  solution  Is  obtained  the  minimum  weight  of  the  structure  is  given  along 
with  values  of  the  design  variables.  In  addition,  the  structures  displace¬ 
ment,  finite  element  stresses  and  reactions  are  given  to  complete  the  des¬ 
cription  of  the  anaylsis. 
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TABLE  1  MATERIAL  PROPERTIES,  CONSTRAINTS 
&  LOADS  -  SEVENTEEN  BAR  TRUSS 


(1)  Material  Properties 

E  =  30.0  x  106  psi 
p  =  .268  lbs.  per  cubic  inch 
0  =  .30 

(2)  Minimum  Si2e  (Size  Constraints) 

2 

Bar  element  area  »  .10  in 

(3)  Allowable  Stress  (Stress  Constraints) 

°1  =  -50000.0  psi,  °u  =  50000.0  psi 
all  elements  except  Nos.  2,  6,  10 

°1  =  -125,000.0  psi,  Ou  =  125,000.0  psi 
elements  Nos.  2,  6,  10 

(4)  Loading  Condition  (Single  Loading  Case) 

Grid  Pts •  3,  5,  7,  9 

Direction  +Y 

Value  -100,000.0  lbs. 


Loads  @  Grid  Pts. 

Direction 

Value 


3,  5,  7,  9 
+Y 

-100,000  lbs. 


Figure  8  Seventeen  Bar  Truss 
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LIST 

10 

T  I  TLE 

1  1 

ClOENT EEN  BAS 

TRUSS-FOUR  LOAPi 

t 

20 

grib 

1 

0.0 

0.0 

0.0 

123456 

30 

GK  I  [i 

7 

0.0 

100.0 

0.0 

123456 

40 

GRIB 

3 

100.0 

0.0 

0.0 

3456 

■  C 

GR  1  n 

4 

10  0.0 

100.0 

0.0 

3456 

60 

GRIB 

5 

700.0 

0.0 

0.0 

3456 

70 

GRIB 

6 

200.0 

100.0 

0.0 

3456 

80 

GR  I B 

7 

300 . 0 

0.0 

0.0 

3456 

90 

GRIB 

8 

..00.0 

100.0 

0.0 

3456 

100 

GRIB 

V 

400.0 

0.0 

0.0 

3436 

110 

OFLOAB$ 

1 

1 

150 

fit  IIM 

YES 

NU 

YES 

OFT 

160 

F  URLE 

1 

i 

1 .0+5 

0.0 

-1.0 

0.0 

170 

FORCE 

1 

5 

1 . 0  +  5 

0.0 

-1 .0 

0.0 
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F  ORCE 

1 

7 

1 .0  +  5 

0.0 

-1.0 

0.0 

190 

FORCE 

1 

V 

1.0+5 

0.0 

-1.0 

0.0 

•..’00 

CONROB 

1 

n 

4 

t 

0.10 

0 

0 

0 

710 

CONROB 

-) 

1* 

3 

o 

0.10 

0 

0 

0 

220 

CONKOB 

3 

1 

3 

1 

0. 10 

0 

0 

0 

230 

CONROB 

4 

3 

4 

1 

0.10 

0 

0 

0 

240 

CONROB 

5 

4 

6 

1 

0.10 

0 

0 

0 

250 

CONROB 

6 

4 

5 

3 

0.10 

0 

0 

0 

760 

CONROB 

7 

3 

5 

1 

0.10 

0 

0 

0 

270 

CONROB 

8 

j 

6 

1 

0.10 

0 

0 

0 

280 

CONROB 

9 

6 

8 

1 

0.10 

0 

0 

0 

290 

CONROB 

10 

6 

7 

4 

0.10 

0 

0 

0 

300 

CONROB 

U 

5 

7 

1 

0. 10 

0 

0 

0 

310 

CONROB 

12 

7 

a 

1 

0.10 

0 

0 

0 

370 

CONROB 

13 

8 

9 

1 

0.10 

0 

0 

0 

330 

CONROB 

14 

7 

9 

) 

0.10 

0 

0 

0 

340 

CONROB 

15 

1 

4 

1 

0.10 

0 

0 

0 

350 

CONROB 

16 

3 

6 

1 

0.10 

0 

0 

0 

360 

CONROB 

1  7 

5 

8 

1 

0.10 

0 

0 

0 

370 

MAT  1 

1 

30. +6 

0.3 

0.268 

+MA7A 

380 

+  MAT  A 

-5.0+4 

f  5  «  04 4 

390 

MAT1 

-> 

30.46 

0.3 

0.268 

+  MATF: 

400 

+  MATB 

-12.5+4 

412. 544 

410 

MAT  1 

3 

30.46 

0.3 

0.268 

+  MA1C 

420 

+  MATC 

-17.5+4 

412.544 

430 

ma  r  i 

4 

30.46 

0.3 

0.268 

+  MATB 

440 

+  MATB 

-12.5+4 

412.544 

450 

CNBBAI A 

Figure  9  OPTFORCE  II  Input  Data  -  Seventeen  Bar  Truss 
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C  AR  C  COUNT  Of  INPUT  DATA 
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i 
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GP  10 
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1 

0 
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FORCE 
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0 . 1 00E ♦ 06  0.0 

0.100E*C6  0.0 
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■1.0000  0.0 
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Figure  10  OPTFORCE  II  Output  -  Seventeen  Bar  Truss 
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Figure  10  (cont'd.) 
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A . 1  INTRODUCTION 


The  programming  aspects  of  OPTFORCE  II  are  described  in  this  Appendix. 
The  information  presented  here  is  geared  to  the  Programmer.  It  is  suffi¬ 
cient  to  fully  describe  the  general  program  logic  and  the  required  peri¬ 
pheral  storage.  All  elemenc  generated  data  is  stored  externally  to  reduce 
core  storage.  A  separate  section  is  devoted  to  the  description  of  these 
fiLes  so  that  I/O  time  may  be  optimized  through  efficient  buffer  descrip¬ 
tion.  Individual  subroutine  write-ups  are  presented  along  with  the  complete 
Fortran  source  listing. 

A  short  description  of  each  routine  is  included  to  aid  in  obtaining 
an  overall  familiarity  with  the  program's  components. 

Finally,  a  discussion  is  provided  for  the  concept  of  dynamic  storage 
which  allows  the  program  to  execute  in  a  variable  storage  environment. 


A.  2  GENERAL  PROGRAM  LOGIC 


The  general  organization  of  the  OPTFORCE  II  program  is  illustrated 
in  Figures  A.l  thru  A. 5.  The  computer  program  consists  of  a  control  pro¬ 
gram  which  includes  four  principal  phases: 

1.  Data  Interpretation  Phase 

2.  Initialization  Phase 

3.  Structure  Cutter  Phase 

4.  Calculation  Phase 

A . 2 . 1  Control  Program  ( Figure  A.l) 

The  "MAIN"  routine  controls  the  program  flow.  The  number  of  work 
storages  in  the  entire  program  is  determined  by  the  variable  NWORK .  This 
parameter  can  be  adjusted  _ ■,>  fit  the  storage  requirements  of  the  problem. 

LINK1  consists  of  phase  1  and  2  -  Data  Interpretation  Phases  and  Initiali¬ 
zation  Phase.  LINK2  consists  of  the  phases  3  and  4  -  Structure  Cutter 
phase  and  Calculation  Phase. 

A .  2 . 2  Data  Interpretation  Phase  (Figure  A. 2) 

The  principal  purpose  of  this  phase  is  to  read  the  NASTRAN  form  input 
data  and  prepare  the  information  for  the  next  three  phases.  Subroutines 
SORT,  ZZ,  and  OPINPT  perform  this  function.  An  echo  print  of  this  input 
is  also  provided  here.  The  SORT  and  ZZ  routines  sort  and  count  the  NASTRAN 
data  records  and  then  form  the  dynamic  storage  constants.  Each  NASTRAN 
label  card  is  processed  and  stored  in  core  by  "ZZ".  An  input  file  data 
unit  NSS1  is  then  written.  This  data  file  is  used  by  OPINPT  which  writes 
on  unit  NTAPE  all  the  data  needed  for  initialization. 

A . 2 . 3  Initializaton  Phase  (Figure  A. 3) 

This  phase  produces  a  variety  of  information  for  the  Structure  Cutter 
and  Calculation  Phase.  Upon  being  activated  by  the  0PTIM2  routine,  the 
"NEWS"  routine  determines  the  dynamic  storage  allocation  for  the  Initiali¬ 
zation  Phase.  The  "AONE"  routine  is  then  activated  and  performs  the 
following  tasks: 

a.  A  load  matrix  is  generated  in  reduced  form  and  written  on 

file  NSS1 . 

b.  Flexibility  and  FORCE  matrices  are  defined  by  the  called 
element  routines  (ELEM1-ELEM4 )  and  are  written  on  file  NSS1.  The  number 
of  forces  in  the  element  and  the  weight  of  the  element  with  a  unit  design 
variable  are  also  written  on  file  NSS1. 

A. 2. 4  Structure  Cutter  Phase  (Figure  A. 4) 

This  phase  produces  four  data  sets  (Ill,  112,  113  and  115)  for  the  Calcu¬ 
lation  Phase.  Routine  TAPE11  defines  the  structure  cutter  matrix  from 
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information  on  unit  NSS1.  As  this  process  continues,  unit  111  is  generated. 
I’pon  calling  routine  AA,  the  matrix  is  solved.  Finally,  routine  S241  controls 
the  generation  of  the  basic  structural  matrices  j) ,  >ji  and  f>  for  each  element  in 
the  structure.  These  matrices  and  information  about  the  mass  matrix  is  written 
to  unit:  113,  t  lie  b^  matrix  is  written  to  unit  112  and  the  matrix  is  written 
to  unit  115. 

A.  2.  5  Calculation  I’hase  (Figure  A.  5) 

The  purpose  of  phase  4  is  to  perform  all  computations  required  for 
static  analysis,  dynamic  analysis  or  structural  optimization.  Input  files  ' 
required  by  this  phase  are  FIFE  112,  113  and  115  while  formatted  output  is 
written  to  FILE  101  and  102  (line  printer).  Phase  4  (OPTFR)  is  called  by 
routine  "LINK2”  which  reads  all  element  and  control  data  from  FILE  Ill 
and  allocates  dynamic  storage  for  phase  4.  Subroutine  "OFTFR"  is  the,'* 
controL  and  initialization  routine  for  the  phase  4  routines. 

Selection  ct  static  or  dynamics  analysis  by  the  user  results  in 
output  of  the  design  variables  followed  by  either  a  static  analysis 
(BASIC)  or  dynamic  analysis  (BAS1CD)  and  output  oc  the  computed  results. 

Control  is  then  returned  to  routine  I.INK2. 

User  selection  of  structural  optimization  follows  the  flow  path  as 
depicted  in  Figure  A.  5  (optimization)  and  is  described  as  follows: 

a.  Initial  Guess  -  User  input  design  variables  may  be  used  as 
the  initial  guess  or  used  as  input  to  the  FSD  routine  to  compute  a  new  de¬ 
sign  using  the  stress  ratio  method.  The  initial  guess  is  printed  and  a 
basic  analysis  is  performed  by  subroutine  "BASIC". 

b.  Compute  Coefficients  -  Subroutines  S432,  S435,  S4310, 

S4323,  S4325  ar.d  S4310A  are  used  to  compute  all  9g/9A,  9g/9X  values  and 
assemble  them  in  a  coefficient  matrix  which  is  input  to  the  linear  pro¬ 
gramming  routine  "ITER". 

c.  Linear  Programming  -  The  function  of  subroutine  ITER  is 
to  determine  if  a  linear  solution  is  feasible  for  minimization  of  the 
objective  function  which  is  subjected  to  the  equality  constraints  9L/9A=:0 
and  9L/9X=0.  In  addition  the  Lagrange  multipliers  (p^,  Us»  Up  and  pw)  are 
constrained  to  be  >0.  If  a  feasible  solution  does  exist  the  active  con¬ 
straints  (p>0.0)  are  used  in  subroutine  S451  to  solve^for  a  new  design. 

The  linear  equations  (Ae+1-Ae)  ^  +  g=0  and  (Xe+|-Xe)  ^+gx=0  are  solved 
simultaneously  for  new  design  Ae+i  an^  redundant  Xe4^  using  routine  S1MQ. 
Subroutine  S451  is  re-iterated  until  convergence  or  a  maximum  number  of 
iterations  is  reached.  If  the  design  converges  and  all  constraints  are 
satisfied,  linear  programming  is  exited,  otherwise,  a  maximum  of  2  ad¬ 
ditional  linear  programming  iterations  are  made.  Whenever  3  linear  pro¬ 
gramming  iterat lore  have  been  completed  without  convergence  or  "ITER" 
indicates  an  infeasible  design  a  full  stressed  design  is  computed  and 
output . 

After  the  statically  constrained  design  has  converged  or  a  fully 
stressed  design  is  generated  any  requested  dynamic  constraints  are  included 
and  a  dynamic  analysis  is  performed  using  subroutine  "BASICD".  Linear 
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programming  is  then  re-entered  to  obtain  a  new  design  with  dynamic  constraints. 


d.  Compute  Lagrange  Multipliers  -  The  design  obtained  from 
linear  programming  is  used  to  compute  the  Lagrange  multipliers  for  all 
active  constraints  using  subroutine  "EQSOL".  If  all  of  the  Lagrange  multi¬ 
pliers  are  >0.0  and  constraint  equations  are  satisfied,  optimization  is 
complete  and  the  final  soltuion  is  printed. 

e.  Newton-Raphson  Procedure  -  In  the  cases  where  there  are 
negative  Lagrange  multipliers  or  constraint  violations,  a  full  Newton- 
Raphson  is  performed  by  Executing  subroutine  S461.  Subroutine  S461  is  an 
iterative  routine  which  takes  derivatives  of  9L/9A,  9L/9X  and  g  to  obtain 
coefficients  for  a  set  of  linear  equations.  These  equations  are  solved 
simultaneously  for  all  AA,  AX  and  Ap  by  subroutine  "GELS".  The  iterative 
process  is  terminated  when  convergence  of  all  variables  is  achieved  or  a 
maximum  number  of  it e  r  a  t i ons  has  been  reached. 

f.  Output  of  Final  Design  -  A  print  of  the  final  design  and 
engineering  parameters  is  accomplished  by  routines  PRIN1  ,  BASIC  and  if  dy¬ 
namic  constraint.-,  are  included  BASICD  is  executed. 

Subroutine  OPFTR  then  returns  to  LINK2  for  program  termination. 
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FIGURE  A.  4  STRUCTURE  CUTTER  PHASE 
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A. 3  EXTERNAL  DATA  SET  STRUCTURE 


This  program  uses  seven  data  sets  during  execution.  The  delivery 
version  of  the  OPTFORCE  II  program  comes  with  the  following  variable 
names  and  real  unit  designations. 


Unit  Name 
NTAPE 


NSS1 

NSS2 

NSS3 

NSS4 


Ill 


112 


113 


114 


Unit  Id 

Usage 

13 

(BINARY) 

Output  from  "MAIN"  and  "OPINPT" 
routines  and  input  to  "AONE"  routine. 

1 

Contains  load  matrix  and  element 
flexibility  matrices  with  assembly 
information. 

Contains  element  B-Force  matrices. 

ITOT  (no.  of  elements)  pairs  of 
information  from  element  matrix 
generation  routines  used  to  initialize 
the  calculation  phase. 

11 

(BINARY) 

Output  from  "TAPE11"  routine  and 
used  in  calculation  phase  contains 
general  description  of  structure 

12 

(BINARY) 

Output  from  "5241"  routine  contains 
b.,  matrix  resulting  from  structure 
cutter. 

13 

(BINARY) 

Output  from  "S241"  routine  contains 
all  bar  (element)  matrices  of  PHI, 

PS I  and  Omega. 

14 

(BINARY) 

Output  from  "S4310A"  routine,  contains 
matrix  element  row,  column,  data,  transpose/ 

symmetry  control  and  summing  control  - 
scratch  file. 


115 

L4 


L7 


15 

(BINARY) 


Output  from  "S241"  routine  contains  Dj 
matrix  resulting  from  structure  cutter. 


9 

(FORMATTED) 


Output  from  "SORT”  and  input  to  "ZZ". 
Contains  identification  record  and 
NASTRAN  input  card  image. 


7  Output  from  "ZZ"  and  input  to 

(FORMATTED)  "0PTIM2".  Contains  card  images  in 
OPTIM  MAGIC  format. 


NPIT/J5/L5  5  Standard  card  input  (80  column  card) 
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NOTE: 


NPOT/JE6/J6/ 

L6/I02  6  Standard  Line  Printer  (132  characters/line). 

101  10  Debug  print  (132  character/record). 

Since  units  5  and  6  are  standard  input  and  output,  they  will  not 
be  considered  in  this  discussion. 
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NTAPE  CONTENTS 


Note:  This  unit  is 
called  15  in 
routine  AONE. 


NTCDS 
Records 
(Title  Cards) 


One  Record 
(Coordinates) 


r. 


(X(I),Y(I),X  (I) ,  1  =  1 ,N2 


One  Record 
Per  Element 
(ITOT  Records) 
(Element  Data) 


IE,  I  ELK  I)  <  IBL'CKL  ,N5(I),N6(I),N7(I),N8(I),N11(I),N13(I)  , 
N15(1),N17(I),(EM(I1,I),I 1=1,11) ,JMAT(I) ,ANGLE(I)  I 


One  Record 
(DOF  Bounded) 


I  (KL(I),I=1, NBOU )  1 
L _ _ _ 1 


NN2L 

Records 

(Load  Components) 


IF 

NDL.NE.O 

Individual 

Constraint 

(DOF-ICON) 


IR, IC,C1 ,C2  j 


(NBDF(I),I=1, NDL 


One 
Record 
If  NDL.NE.O 
(Displacement 
Limits) 


(DISPU(I),DISPL(1) ,1=1, NDL 


If  IRST 
. EOl  or 
2  or  5 
One  Record 
(Input  Element 
Design  Paramets) 


|jALL(I),  1-1,1  TO  T) 
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A.  3.1  NSS1  CONTENTS 


1st 

Record 

(Problem 

Control 

Data) 

2nd 

Record 

(Stress  Limits 
Min  Size 
Element  Type) 

Dummy 
Read 
If  NDL 
NE  0 

(Displacement 

Constraint 

DOF's) 

If  NDL 
NE  0 

(Displacement 

Limits) 

Program 
Flow  Controls 

Load 

Matrix 

1  Record  Per 
Element 
(Element 
Mat  r  ices 
and  Data) 

If  IRST 
NE  0 
or  NE  3 
or  NE  A 
or  NE  6 


NRDF, NLOAD , ITOT , NALD , NDL , NSG , NAA , NBOV , NDOF, N2 , NELI , IRST . 
NMODES 


(D,SIGU(I),ALL(1) ,1=1, ITOT) , I BS , (D , 1=1 , IBS ) 
, ( (D, I=1,11),J=1, ITOT), ( LELT ( I ) , 1=1 , ITOT ) 


(DISPU(I),DISL(l),l=l,NDL 

(N22I(I)=I~1,NDL 

(DISPU(I),DISL(I),I=1, NDL 

IREST , IPRI , IPRS, ISTRM, 100,MAX50 ,NLINK1 , NICON 

NGRUPS , NODES 4 , NOSYMM , NANTI , NNONS , C0005 , WISTAR , IRST , NMODES 


J1 , (DELTA( I) , 1=1, Jl) ,NBOU, (KL(I) , I=l,NBOU) 


KS , KEL , AREA ( NEL ) , (L 1 3 ( I ) , 1=1 , KEL 
, ( (B( I, J ) , 1=1 ,KEL) , J=1 ,KS) , 

D,D,D,D, ( (F (I , J) ,1=1 ,  KS) , ,J=1,KS) ,D,NNE, ( (D, 1=1 ,NNE) ,J=1,NNE) 
(D, 1=1, KEL) , (D, 1=1 ,KEL) , NNNO, D, D, D,D,D,D,D,D,D, 

WT(NEL) ,XC(NEL) ,YC(NEL) _ 


ALL2(I), 1=1, ITOT 


Input  Starting 
Design  Parameters 
(OPDVIR  Card) 


Written  by  Routine  AONE 

Read  by  Rc  it  ines  TAPE11  and  S241 
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A. 3. 2  Ill  CONTENTS 


1st 

Record 

(Problem 

Constants) 


2nd 

Record 

(Element  Data 
Size  and 
Stress  Limits) 


3rd 

Record  I - 

(Load  Matrix  (DELTA (I) ,T=1,J1),(DL(J),J=1 ,NBOU) 

and  Assembly  ' — — - 

Data) 


4th 

Record 

(Displacement 
Constraints 
and  DOF 
Affected ) 


(DiSPU(I) ,DISPL(I) ,N221(I) , 1=1 ,NDL 


(  IELT ( I )  ,  ALL  ( I )  ,WT ( 1  )  ,AREA(I)  ,XC(I)  ,  YC  ( 1 )  ,ALI. 2  (I ) , 
SIGU( 1 3 , 1=1 , ITOT) 


M , J 1 , NR  ED , NBOU , NDO  F , N2 , NEL I , NDTNX , NDL , IRST , NMODES 
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A. 3. 3  112  CONTENTS 


One 

Record 

Per 

Element 

(P1I 

Matrix) 

One 

Record 
(B1r 
Mat r  Lx) 


(Bll(J),  3=1,  ( 5*NX) ) 


(BIR(.I)  ,J=1,  (NBOl!*NX)  1 


Written  bv  S241  routine 


A. 3.4  113  CONTENTS 


One 

Record 

Per 

Element 
(Normalized 
Structural 
Bar  Matrices 
and  Mass  Matrix 
Data) 


((PHIB(l.J) . 1  =  1 ,NX) . 1=1 .NX) , 
((PSIB(I.l), 1=1 .NX)  , J  =  1 , NDN)  , 

( (OMEGAB( I ,J) . 1=1, NDN) ,J=1,NDN) 
(2(1), 1=1. NDN) 


Written  by  S241 
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A. 3. 5  114  CONTENTS 


One 

Record 
Per 

Matrix 
Element 
(Matrix  -  Row 
.‘lumber.  Column 
Mumber,  Transpose/ 

Symmetry  Control, 

Element  Data, 

Summary  Control) 

Scratch  file  written  by  S4  31C  \  routine 
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A. 3.6  115  CONTENTS 


One 

(Dll(J) ,J=1» (5*NDN) ) 

Record 

Per 

Element 

(Dlj. 

Matrix) 

One 

[  (DIR(J) , J=l» (UBOU*NDN) ) 

Record 

(DV  , 

Matr  ix  ) 


Written  by  S241  routine 


A. 3.  7  L4  CONTENTS 


(80  Character 
card  images) 


These  records  can  be  in  any  order  on  the  file.  It  represents  data  in 
"pairs"  of  records.  Each  pair  represents  one  specific  type  of  data 
items  determined  by  the  value  "NKIND"  and  "NCOUNT"  which  are  written 
on  the  first  record  of  .he  pair.  R 

Record  one  of  the  pair:  NKIND,  NCOUNT 

Record  two  is  determined  by  NKIND  as  shown  in  the  following  table: 
NKIND  Second  Record 


1 

2 

4 

6 

8 

(2 

13 

16 

19 

20 
21 

23 

24 

27 

28 

30 

31 
33 

39 

40 
43 

46 

47 

48 

49 

50 

51 


BULK  data 

TITLE 

CONROD 

00DMEM1 

I  I 

CROD 

CSHEAR 

CTRMEM  " 

END* 

ENDDATA 

FORCE 

GRID 

ICON 

MAT1 

OPLOAD 

0PDV1R 

0PT1M 

P0DMEM1 

PROD 

PSHEAR  " 

PTRMEM 

SPC  " 

SPC1 

SPC1  THRU  " 

CWEB 

PWEB 

BEGIN  BULK 
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A. 4  Subroutine  Write-Ups 


Each  subprogram  of  OPTFORCE  il  is  described  in  this  section.  In¬ 
cluded  with  each  description  is  a  statement  declaring  the  size  of  the 
subprogram.  This  number  is  intended  as  a  guide  only,  as  it  reflects 
the  storage  requirement  on  an  IBM/ 370/ 3031 ,  FORTRAN  compiler. 

In  this  manual  the  subroutine  write-up  ire  presented  in  alpha¬ 
betical  order,  separately  for  T.1NK1  routines  and  LINK.2  routines. 


\ 


1 
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Routine 

ADJUST 

AONE 

EGMAT 

EMAT 

ELEM1 

ELEM2 

ELEM3 

I 

ELEM4  j 

EXIT 

FOMO 

1NSPC 

ISOTR 

LINK! 

MINV 

NEWS 

OPINPT 

OPTIM2 

READI 

SORT 

SPCSUB 

TMAT 

WRITEL 

XTRAK 

ZZ 


A. 4.1  BRIEF  ROUTINE  DESCRIPTIONS,  LINK1 

Name  Purpose 

Adjust  each  NASTRAN  form  field  so  that  it  is  either  left  or 
right  adjusted. 

Control  initialization 

Generates  EG  matrix  for  triangular  and  quad  plates 

Generate  elasticity  matrix  for  either  orthotropic  or  isotropic 
properties,  membrane  triangle,  quad 

Axial 

Element  Shear  web 

Matrix 

Generators  |  Triangle 

[  Quad 

Provide  for  final  stop  in  program 
Process  Force  Cards 

Stores  data  obtained  from  input  card  into  working  storage 
Given  any  two  of  G,  E  and  p  computes  third  quantity 
Control  program  flow  for  PHASE1  and  PHASE2 
Matrix  Inverse 

Allocate  storage  for  initialization 

Interpret  report  from  input  and  place  results  on  NTAPE 
Interprets  input  and  initializes  data  -  NEWS,  OPINPT,  EXIT 
Reads  and  modifies  input  data 

Sort  and  count  data  based  on  LABEL  information 
Process  SPC  (single  point  constraint)  cards 

Generate  transformation  matrix  for  triangular  plate  orthotropic 
angle 

Tests  character  of  element  connection  card  and  writes  element 
information  on  to  file 

Interpret  degree  of  freedom  informations 

Generates  OPTIM  data  which  is  input  by  NASTRAN  form  input  cards 
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1.  Subroutine  Name 


ADJUST 


2.  Purpose 


Adjust  each  NASTRAN  form  field  so  that  it  is 
either  left  or  right  adjusted. 


3.  Procedure; 


Each  character  of  the  word  is  tested  to  see  if  it 
is  blank.  When  a  non-blank  is  met,  it  is  shifted 
to  the  end  of  the  word. 


4.  Input  Arguments 


6.  Error  Returns 


7.  Calling  Sequence 


8.  Subroutine  User 


9.  Storage  Required 


WORD  -  Input  word  to  be  either  right  or  left  adjusted 
Icode  -  I  code  =  8  specifies  that  the  word  should  be 
right  adjusted 

-  Icode  =  1  specifies  that  the  word  should  be 
left  adjusted 

Word  right  adjusted,  word  is  stored  back  into  WORD 


CALI.  ADJUST  (Word,  ICODE) 


READI 


(642  Bytes)  161  words 
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1.  Subroutine  Name: 


AONE 


2.  Purpose:  Assemble  all  initialization  information  for  the 

calculation  phase  of  the  program. 

3.  Equations  and  a)  Print  out  title  cards. 

Procedures:  b)  Retrieve  (13)  coordinate+element  data.  If 

requested  print  this  data  out. 

c)  Recalculate  boundary  conditions  if  symmetry 
plane  nodes  are  specified. 

d)  Create  load  matrix  from  input  loads. 

e)  Define  individual  constraints. 

f)  Call  element  routines. 

g)  One  unit  of  data  in-  written  with  information 
for  the  calculation  program. 


4. 

Input 

Arguments : 

15 

(Unit  13)  tape 

NPOT 

(Unit  6)  printer 

NSS1 

(Unit  1  I/O) 

N2 

No.  of  nodes 

NDL 

No.  of  Individual  Constraints 

NSYM 

No.  of  load  conditions 

ITERN 

Max.  No.  of  iterations 

NTCDS 

No.  of  title  cards 

5. 

Output 

Arguments: 

C1INP 

C2INP 

|  Convergence  limits 

NRDF 

No.  of  reduced  DOF 

NBOU 

Total  No.  of  constrained  DOF 

IRST 

Calculation  control 

6. 

Error 

Returns : 

None 

7. 

Calling  Sequence: 

CALL  AONE(X,Y,Z,N5,N6,N7,E,AMU,N8,N17 

S IGL , R , N 1 1 , N1 3 , N1 5 , NOAL , N0AL2 , KL , KL2 , LNOD , LN0D2 , 
DISPU , DISPL,NBDF,NBDF2,NSE, DELTA, 15 ,NP0T,NSS1,NSS2 , 
NSS  3 , NSS4 , IBUKL , IELT, IBK , N1 , N2 , NAA , NDL , NSG , NNZL , 

I ELI , NREACT , Cl I NP , C2 INP , KNOAL , IBKGP , EM , JMAT , ANGLE) 


8. 

Input  Tapes: 

NT APE  (unit  8) 

9. 

Output  Tapes: 

NSS1  (unit  1) 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(17840  bytes)  4460  words 

12. 

Subroutine  User: 

NEWS 

13. 

Subroutine  Required: 

EL  EMI 

ELEM2 

ELEM3 

ELEM4 

ISOTR 

14.  Remarks:  None 
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1. 

Subroutine  Name: 

EGMAT 

2. 

Purpose : 

Generate  Eg  matrix  for  triangular  and  quad  plates 

3. 

Equation  and 
Procedures: 

Eg  =  TETT 

4. 

Input  Argument 

E  =3x3  matrix 

T  =3x3  transformation  matrix 

5. 

Out  Arguments 

Eg  =  matrix  defined  by  equations  above 

6. 

Error  Returns 

None 

7. 

Calling  Sequence 

Call  EGMAT  (EG,  E,  T) 

8. 

Subroutine  User 

ELEM3,  ELEM4 

9. 

Storage  Required: 

(716  bytes)  179  worm; 
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1. 


Subroutine  Name: 


EMAT 


2.  Purpose: 


To  generate  elasticity  matrix  for  either  orthotropic 
or  isotropic  properties 


3.  Equations  and 
Procedures : 


Orthotropic 
(IMAT  =  2) 

E  = 

G11 

C12 

G13 

C2L 

G22 

G23 

C3i 

C32 

C33 

Isotropic 
(IMAT  =  1) 

E  = 

E 

X 

M 

XX 

o 

M  Ey 
xy 

Ey 

0 

0 

0 

C 

xy 

4.  Input  Arguments:  E^,  Ey,  -  material  properties 

EM  -  contains  G..  values 

IMAT  -  Orthotropy  code  =  2  orthotropic 
Orthotropy  code  =  1  isotropic 


5. 

Output  Arguments: 

E  -  output  3x3  matrix  defined  by  equations  above 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  EMAT  (E  ,  E  M  ,  G  ,  E,  EM,  IMAT) 

x  y  xy  xy 

8. 

Subroutine  User: 

ELEM3,  ELEM4 

9. 

Storage  Required: 

(532  bytes)  144  words 
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1.  Subroutine  Name: 


EL  EMI 


2.  Purpose: 


To  compute  the  elemental  force  and  flexibility 
matrices  for  an  axial  force  element. 


3.  Input  Arguments:  X,  Y,  Z 

E 

RHO 

N5,  N6 
NEL 
N1 
N2 

4.  Output  Arguments:  M 

M(25) 

M(26)=1“KKK 
S  ( 1 , 1 ) 

F  ( 1 , 1  ) 

WT 

ALEN 

ENEL 


Node  point  coordinates 

Modulus  of  elasticity 

Density 

Element  nodes 

Element  number 

Degrees  of  freedom/node 

Total  number  of  nodes 

Unreduced  direction  numbers  for  degrees 

of  freedom  of  element 

Order  of  B-Force  matrix 

Number  of  force  components 

Start  of  P.-Force  matrix 

Stait  of  flexibility  matrix 

Weight  of  element  with  unit  design  variable 
Length  of  element 
Modulus  of  elasticity 


5.  Error  Returns:  None 


6. 

Calling  Sequence: 

CALL  EL EMI  (X, Y , Z , E, AMU , RHO , N7 ,N5 ,N6 ,N8 ,N9 , N10, Nil , 
N12,R,M,NEL,N1,N2,B,F,ALEN,KKK,ENEL,WT) 

7. 

Input  Tapes: 

None 

8. 

Output  Tapes: 

None 

9. 

Scratch  Tapes: 

None 

10. 

Storage  Required: 

(1234  bytes)  321  words 

11. 

Subroutine  User: 

AONE 

12.  Subroutine  Required:  None 


13:  Remarks: 


The  force  matrix  provides  one  force  component  S, . 


1.  Subroutine  Name: 


ELEM2 


2.  Purpose: 


To  compute  the  element  force  and  flexibility  matrices 
for  a  shear  web  element. 


3.  Input  Arguments:  X,Y,Z 

E 

AMU 

RHO 

N5,N6,N7,N8 

NEL 

N1 

N2 

4.  Output  Arguments:  M 

M(25) 

M(26)-1=KKK 

8(1,1) 

F(l,l) 

WT 

AREA 

ENEL 


Node  point  coordinates 
Modulus  of  Elasticity 
Poisson's  Ratio 
Element  Density 
Element  nodes 
Element  number 
Degrees  of  freedom/ node 
Total  number  of  nodes 

Unreduced  direction  numbers  for  degree 

of  freedom  element 

Hows  of  force  matrix 

Number  o<  stress  components 

Start  of  force  matrix 

Slai  1  of  flexibility  matrix 

Weight  of  element  with  unit  design 

variable 

Surface  area  of  element 
Modulus  of  elasticity 


3. 

Error  Returns: 

None 

6. 

Calling  Sequence: 

CALL  ELEM2  (X,Y,Z,E,AMU,RHO,N7,N5,N6,N8,N9,N10,N11 
N12, R,M, NEL, N1,N2,B,F, AREA, KKK, ENEL, WT) 

7. 

Input  Tapes: 

None 

8. 

Output  Tapes: 

None 

9. 

Scratch  Tapes: 

None 

10. 

Storage  Required: 

(3648  bytes)  912  words 

11. 

Subroutine  User: 

AONE 

12.  Subroutine  Required:  None 

13.  Remarks:  A.  This  element  can  be  defined  by  2  or  4  nodes. 

B.  The  force  matrix  provides  one  force  component. 


83 


L.  Subroutine  Name: 


ELEM3 


2.  Purpose: 


3.  Input  Arguments: 


4.  Output  Arguments: 


5.  Error  Returns: 


To  compute  the  element  force  and  flexibility  ortho¬ 
tropic  matrices  for  a  triangular  piate  element,  with 
properties  and  material  angle  variation. 


EE 

AMU 

RHO 

N5.N6.N7 

NEI. 

N1 

N2 

EM 

IMAT 

ANGLE 


Node  point  coordinates 
Modulus  elasticity 
Poisson's  ratio 
P.-nsi :  y 
Element  nodes 
Element  number 
Degrees  of  freedom/node 
Total  number  of  nodes 
Material  properties 
Orthotropic  code 
Material  angle 


M 

M(25) 

M(26)=3=KKK 

B(1,D 

F(l,l) 

WT 

AREA 

E 


Reduced  direction  number  for  degrees 

of  freedom  of  element 

Rows  of  force  matrix 

Number  of  force  components 

Start  of  matrix  force 

Start  of  matrix  flexibility 

Weight  of  element  with  unit  design 

variable 

Surface  area  of  element 
Elasticity  matrix  (3  x  3) 


None 


6 . 

Calling  Sequence: 

CALL  ELEM3  (X , Y , Z , EE, AMU , RHO, N7 , N5 , N6 ,N8 , N9 , N10, Nil , 
N12, R,M,NEL,N1,N2, CM, 1MAT, ANCLE, B,F, AREA, KKK.E.WT) 

7. 

Input  Tapes: 

None 

8. 

Output  Tapes: 

None 

9. 

Scratch  Tapes: 

None 

10. 

Storage  Required: 

(3510  bytes)  878  words 

11 . 

Subroutine  User: 

AONF. 

12.  Subroutine  Required:  EGMAT, EMAT,MINV,TMAT 

13.  Remarks:  Three  components  of  force  at  the  midpoint  are  provided. 
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1.  Subroutine  Name: 


ELEM4 


? 


3. 


4. 


5. 


6. 


7. 

8. 

9. 


10. 


11. 


12. 


13. 


Purpose:  This  routine  computes  the  element  force  and  flexi¬ 

bility  matrices  for  a  quadrilateral  plate  element 
with  property  and  material  angle  variation. 


Input  Arguments: 


X ,  Y ,  Z 

Note  point  coordinates 

E 

Modulus  of  elasticity 

AMU 

Po isson' s  ratio 

RHO 

Den si  tv 

N5,N6,N7,N8 

Element  nodes 

NEL 

Element  number 

N1 

Degrees  of  f reedotn/node 

N2 

Total  number  of  nodes 

I  MAT 

Orthotropic  code 

ANGLE 

Material  angle 

EM 

Material  properties 

Output  Arguments:  M 

M<25) 

M(26)=5 

B (1 , 1 ) 

F(l,l) 

AREA 

WT 

VYo 

NNE=3 

E 


Unreduced  direction  numbers  for  degrees 

of  freedom  of  element 

Rows  of  force  matrix 

Number  of  force  components 

Start  of  force  matrix 

Start  of  flexibility  matrix 

Surface  area  of  element 

Weight  of  element  with  unit  design 

variable 

Coordinates  of  centroid  (local  coordinates 
Elasticity  matrix  order  (3  x  3) 


Error  Returns: 

Calling  Sequenc-  CAL  ELEM4  (X,Y,Z,E,AMU,RHO,N7,N5,N6,N8,N9,N10,N11, 

N12, R,M.NEL,N1,N2, EM, LMAT, ANGLE, B,F, AREA, NNE,E,WT, 
XC.YC) 


Input  Tapes: 

Output  Tapes: 

Scratch  Tapes: 
Storage  Required: 
Subroutine  User: 
Subroutine  Required: 
Remarks : 


None 

None 

None 

(5688  bytes)  1422  words 
AONE 

EGMAT , EMAT , MINV , TMAT 

Five  components  of  force  are  defined  at  the  midpoint. 
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Subroutine  Name : 

EXIT 

Purpose : 

To  p 

rovide  for 

Equations  and 
Procedures : 

Ente 

r  and  stop 

Input  Arguments: 

None 

Output  Arguments: 

Nolle 

Error  Returns: 

Nolle 

Calling  Sequence: 

CALL 

EXIT 

Input  Tapes: 

None 

Output  Tapes: 

None 

Snatch  Tapes: 

None 

Storage  Required: 

Subroutine  User:  MAIN,0PTIM2 


in  the  program. 


Subroutine  Required:  None 


A 128  235  FORCE  METHOD  OPTIMIZATION  II  VOLUME  II  USER'S  MANUAL 
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JNCLASSlFIED 


t 


p 


! 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 196 J-A 


1.  Subroutine  Name 


FOMO 


2.  Purpose  Process  force  cards 

3.  Equations  and  The  cards  are  read  from  file  L5  and  counted  by  a 

Procedures:  counter  KFORCE.  The  input  values  are  (Force  1(1). 

i  ■  1,  2),  F,  and  (F2(i),  i  =  1,  2,  3);  (Force  2u) 
F  -  Fn(i),  i  -  1,  2.  3). 

4.  Input  Arguments:  OUT  =  storage  lor  reading  data  from  label  field 

( '  'i  8  columns) 

FORCE  1  •=  «.-*L  id  number,  grid  point  number 

1,5  “  input  file  number 

L6  =  output  file  number 

NFO  =>  number  of  loads 

5.  Output  Arguments:  Force  2(i)  for  each  component 

6.  Error  Returns:  None 

7.  Calling  Sequence:  Call  FOMO  (FORCE  1,  FORCE  2,  L5,  NFO,  K  FORCE,  OUT, 

L6) 

8.  Subroutine  User  ZZ 

9.  Storage  Required:  (712  bytes)  178  words 


Note:  F  is  a  scale  factor  for  force  components  read  from  cards. 
FN(i),  i  ■  1,  2,  3  component  of  force  read  from  cards. 
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1.  Subroutine  Name: 


INSPC 


2. 

Purpose : 

Stores  data  obtained  from  input  card  into  working 
storage. 

3. 

Equations  and 

SPCINF  -  Out  storage,  SPCTYP 

Procedures: 

4. 

Input  Argument: 

SPCTYP  *  Special  code  for  type  of  card 

OUT  *  Storage  for  card  data 

NSPCS  =  Maximum  number  of  special  input  cards 

5. 

Output  Argument 

SPCINF 

NSCARD  -  Counter  for  no.  special  input  cards. 

6. 

Error  Returns: 

7. 

Calling  Sequence: 

Call  INSPC  (NSCARD,  OUT,  NSPCS,  SPCTYP,  SPCINF) 

8. 

Subroutine  User : 

ZZ 

9. 

Storage  Required: 

(418  bytes)  105  words 
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1.  Subroutine  Name: 

2 .  Purpose : 

3.  Equations  and 

Procedures: 

4.  Input  Arguments: 

5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Input  Tapes: 

9.  Output  Tapes: 

10.  Scratch  Tapes: 

11.  Storage  Required: 

12.  Subroutine  User: 


ISOTR 

Given  any  two  quantities  of  E,  G  and  p  compute  the 
third  if  not  defined  («  0.0) 

E  -  2G*(1  +  p) 

G  -  E*(l  +  p)/2 

p  -  (E/2G)  -  1.0 

E  »  Elastic  modulus 

G  *  Shear  modulus 

p  =  Poisson's  ratio 

See  purpose. 

None 

Call  ISOTR  (E,NU,G) 

None 

None 

None 

(360  bytes)  90  words 
AONE 


13.  Subroutine  Required:  None 


1.  Subroutine  Name: 


LINK1 


2.  Purpose:  Control  program  flow  for  the  data  interpretation 

and  initialization  phases. 

3.  Equations  and  Call  sort  routine  which  reads  input  NASTRAN  data 

Procedures:  and  adjusts  fields  and  writes  on  Unit  9. 

'•’he  ZZ  routine  is  called  next  where  the  data  is 
interpreted. 

Next  the  0PTIM2  routine  is  activated  which  controls 
the  initialization  phase. 

Control  is  then  returned  to  the  main  routine. 


4. 

Input  Arguments: 

W  -  Dynamic  storage 

WORK  -  Maximum  size  of  dynamic  storage 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  LINK1  (W, NWORK) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage.  Required: 

(2352  bytes)  588  words 

12. 

Subroutine  User: 

Main 

13. 

Subroutine  Required: 

SORT,  0PTIM2,  ZZ,  EXIT 
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1.  Subroutine  Name: 


MINV 


2.  Purpose:  Invert  a  matrix 

3.  Equations  and  Uses  the  standard  Gauss-Jordan  method  in  which  the 

Procedures:  inverted  matrix  is  stored  back  on  itself. 

4.  Input  Arguments:  A  -  Matrix  to  be  inverted 

N  =*  Order  of  matrix 

D  =  Determinant  value 

L  «  Work  vector  of  length  N 

M  =  Work  vector  of  length  N 

5.  Output  Arguments:  A  -  Contains  the  inverted  matrix 

6.  Error  Returns:  If  D=0,  matrix  is  singular 

7.  Calling  Sequence:  CALL  MINI’  (A,N,D,1.,M) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 


10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1440  bytes)  360  words 

12. 

Subroutine  User: 

ELEM3,  ELEM4 

13. 

Subroutine  Required: 

None 

14. 

Remarks : 

None 
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1 .  Subroutine  Name : 


NEWS 


2 .  Purpose ; 


3.  Equations  and 
Procedures: 


To  determine  max  core  available  vs  problem  size  and 
set  up  dynamic  storage  subscripts  for  the  initiali¬ 
zation  phase  of  the  program. 

Storage  requirements  for  the  arrays  used  by  the  AONE 
(utilization  routine)  are  computed  from  the  input 
parameters. 

If  the  storage  r-.'juired  exceeds  MAXCOR(the  size  of 
the  W  array)  then  IER  is  set  to  one,  a  message  is 
written  describing  the  additional  storage  required 
and  the  routine  returns  to  the  main  routine. 

If  the  storage  available  is  adequate,  locations  in 
the  W  array  are  computed  for  the  required  arrays  and 
routine  AONE  is  called. 


Calling  Arguments: 
(*  Input) 


*N2 

*ITOT 

*NBOU 

*NSYM 

*NDL 

*NTAPE 

*NPQT 

*NSS1 

*NALD 
*ITERN 
*NTCDS 
*IREST 
I  ELI 


NREACT 


*C1INP 

*C2INP 

NRDF 

*IRST 

*MAXCOR 


Control  for  degrees  of  freedom  of  a 
node  set  to  3  internally. 

No.  of  nodes  in  problem. 

Number  of  elements  in  problem. 

Total  number  of  constrained  degrees  of 
freedom  for  a  symmetric  load  condition. 
Number  of  load  conditions. 

Number  of  individual  constraints. 

Input  unit  for  AONE  routine  (Unit  13). 
Print-out  data  set  (Unit  6). 

Output  data  set  for  initialization 
phase  of  program  (Unit  1). 

Total  number  of  load  conditions  (NSYM) . 
Maximum  number  of  iterations. 

Number  of  title  cards. 

Print  control  for  printing  input. 

Integer  array  indicating  number  of  each 
type  element,  e.g.,  IELI(4)«N  type  4 
elements. 

Integer  array  six  elements  long  indi¬ 
cating  total  number  of  constraints,  e.g. 
NREACT(l)-  number  of  constrained  U 
components. 

Convergence  control. 

Convergence  control. 

Number  of  reduced  degrees  of  freedom  for 
summetric  load  condition. 

Calculation  control  which  determines 
use  of  OPDVIR  section. 

Size  of  ’W  array  available  for  dynamic 
storage. 

Array  used  for  dynamic  storage. 

Error  control. 


5.  Error  Returns: 


lERjtQ 


Not  enough  core  available  to  complete 
initialization. 


6. 


7. 


8. 


9. 


10. 


11. 


12. 


13. 


Calling  Sequence: 

Input  Tapes: 

Output  Tapes: 

Scratch  Tapes: 
Storage  Required: 
Subroutine  User: 
Subroutine  Required: 
Remarks: 


CALL  NEWS (N1,N2,  NSG,NAAALL,NDL,NAA,OTAPE,NPOT, 

NSS1 , NSS2 , NSS3 , NSS4 , IBUKL , NNZI. , IELI , NREACT , Cl INP , 

C2INP, KNLMAX.MAXCOR.W) 


None 


None 

(3078  bytes)  770  words 

0PTIM2 

AONE 


None 
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1.  Subroutine  Name:  OPINPT 

2.  Purpose:  To  complete  NTAPE  (Unit  13);  input  to  the  AONE 

routine . 

3.  Equations  and  Process  input  needed  by  the  AONE  routine  in  this  order. 

Procedures: 

COORD 

EL  EM 

OEXTERN 

LINKS 

BOUND 

OLOADS 

ICON 

GCON 

OPDVIR 


4.  Input  Arguments: 


5.  Output  Arguments: 


Place  all  information  or  NTAPE  (Unit  13).  Print 
diagnostic  messages  if  necessary. 


NPIT 

NPOT 

NTAPE 

N2 

NDL 

NNZL 

NSS1 

NREF 

NGRD 

NDOFPN 

NALD 

ITOT 

NIBCP 


Unit  5  (card  reader) 

Unit  6  (printer) 

Unit  13  (tape) 

Number  of  nodes 

Number  of  individual  constraints 
Number  of  input  cards  for  loads 
Unit  1  (tape) 


Info  from  System  section  input 


NREACT  Array  indicating  number  of  bounded  DOF. 

IER  Error  indicator. 


6.  Error  Returns:  IER*0 


7.  Calling  Sequence:  CALL  OPINPT(NPIT, NPOT, NTAPE, NREF, NGRD, NDOFPN, NALD, 

ITOT , NIBCP , N1 , N2 , NSGIN , NDL , NAA , NDLA , NREACT , NNZL , IELI , 
IRST , NAANUM , NSG , NS  S 1 , KNLMAX , I ER , GRID , IBOUND , N BOUND , 
INODE, NBDF,DISPU,DISPL,LINKB,LINKN,ILINK,DVIR,IBUCKL, 
NAM, IBO, IG1, IG2 , IG3, IBGP) 


8.  Input  Tapes:  NSS1 

9.  Output  Tapes:  NTAPE 

10.  Scratch  Tapes:  NSSL  (Unit  1)  used  to  merge  ELEM  and  OEXTRN  sections. 

11.  Storage  Required:  (14922  bytes)  3731  words 

12.  Subroutine  User:  0PTIM2 
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13.  Subroutine  Required:  None 

14.  Remarks:  None 
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1.  Subroutine  Name: 


0PTIM2 


2.  Purpose: 


To  control  the  program  execution  for  phase  2. 


3.  Equations  and  Identify  all  external  files  and  rewind  them. 

Procedures:  Read  REPORT,  TITLE,  SYSTEM  and  0PT1M  sections  to 

define  variables  needed  for  dynamic  storage. 

If  there  are  any  title  cards  read  them  and  put  them 
on  NTAPE  (Unit  8). 

Print  message  if  any  sections  are  out  of  order. 

Call  OPINPT  routine  to  read  rest  of  input  sections. 
If  there  is  any  error  (lER^O)  call  exit. 

Call  NEWS  routine  to  perform  dynamic  storage  alloca¬ 
tion  and  call  routines  to  perform  initialization. 

If  there  is  any  error  (IER^O)  call  exit. 


4. 

Input  Arguments: 

None 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

(IER^O)  call  exit. 

7. 

Calling  Sequence: 

Call  0PTIM2  (WORK,  NPIT) 

8. 

Input  Tapes: 

WORK  ■  Work  storage 

NPIT  =  Input  file  number 

9. 

Output  Tapes: 

NTAPE  (Unit  13) 

10. 

Scratch  Tapes: 

NSS1  (Unit  1) 

11. 

Storage  Required: 

(3784  bytes)  946  words.  The  work 
dynamically  locate  arrays. 

12. 

Subroutine  User: 

LINK1 

13. 

Subroutine  Required: 

OPINPT 

NEWS 

EXIT 

14. 

Remarks : 

None 
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1. 

Subroutine  Name 

READI 

2. 

Purpose: 

Reads  and  modifies  input  data. 

3. 

Equations  and 

Reads  the  card  data  from  file  J5,  right  adjust  the 

Procedures: 

data  fields. 

counts  each  data  type,  prints  the  data 

and  finally 

stores  the  modified  data  on  file  J6. 

4. 

Input  Argument: 

LABEL 

--  array  of  BCD  label  codes 

I  LAB 

-  array  of  label  integers 

ISPECL 

=  array  of  special  labels 

NILAB 

=  total  no.  of  labels 

J6 

=  output  file  number 

NSPECL 

=  total  no.  special  labels. 

5. 

Output  Arguments: 

NCARDS 

=  array  of  card  counters 

I.7CASE 

=  code  to  indicate  that  file  is  =  7. 

6. 

Error  Returns : 

None 

7. 

Calling  Sequence: 

CALL  READI ( LAB  EL , I LAB , I SPECL , NC ARD  S , NILAB , NSPECL , 

J6.L7CASE) 

8. 

Subroutine  User: 

SORT 

9.  Subroutines  Used:  ADJUST 

10.  Storage  required:  (1940  bytes)  485  words 
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1.  Subroutine'  Sume:  SORT 

2.  Purpose:  Sort  and  count  data  based  on  LABEL  information. 

3.  Equations  and  Using  READI,  the  data  is  read  and  counted.  The 

Procedures:  final  counters  are  then  modified. 

4.  Input  Arguments:  1,4  =  file  number  for  storage  of  sorted  input  deck. 

5.  Output  Arguments:  MEL  Total  no.  Elements 

MGR  total  no.  Grid  points 

MMT  Total  no.  Materials 

MOGCON  Total  no.  Generalized  Constraints 

MIC0N8  Total  no.  Individual  Constraints 

MILINKS  Total  r.o.  Links 

MFO  Total  no.  Forces 

'WO  Total  no.  Moments 

MLOAOS  Total  no.  Loads 

6.  Error  Returns:  None 

7.  Calling  Sequence:  Call  SORT(MEL,MGR,MMT, MOGCON, MI C0N8 .MILINKS , MFO, 

MMO , MLOADS , MGROUP , MSPCS , L4 , L7CASE) 

8.  Subroutine  User:  LINK1 

9.  Subroutines  Used  READI 


10.  Storage  Required: 


(1578  bytes)  395  words 


1. 

Subroutine  Name: 

SPCSUB 

2. 

Purpose : 

Process  SPC 

(single  point  constraint)  cards. 

3. 

Equations  and 

Boundary  information  is  processed  as  read  in  OUT  (I) 

Procedures: 

and  NOUT  (I) 

.  This  information  is  interpreted  and 

stored  in  LBOUND. 

4. 

Input  Arguments: 

NOUT 

Input  data  storage 

OUT 

Input  data  storage 

NSPC 

No.  SPC  cards 

KWORI) 

Work  storage 

NUM 

Work  storage 

NKIND 

Type  of  boundary  information  available 

NGR 

Total  no.  grid  points 

5. 

Output  Arguments: 

'.BOUND 

Boundary  array  information 

"BCARD 

Counter  of  boundary  information 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  SPCSUB 
NKIND, NGR) 

(  NOUT , OUT , NSPC , KWORD , NUM , LBOUND , NBCARD , 

8. 

Subroutine  User: 

ZZ 

9. 

Subroutines  Used: 

XTRAK 

10. 

Storage  Required: 

(1602  bytes) 

401  words 
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1.  Subroutine  Name: 

2.  Purpose: 

3.  Equations  and 

Procedures: 

4.  Input  Arguments 

5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence 

8.  Subroutine  User: 

9.  Subroutines  Used: 

10.  Storage  Required: 


TMAT 


Generate  transformation  matrix  for  triangular  plate 
orthotropic  material  angle. 


T 


„  2 
Cos  6 

Sin2e 

-2  SineCose 


2 

Sin  e  SineCose 

2 

Cos  e  -SineCose 

2  2 

2  SineCose  Cos  e-Sin  e 


THETA  =  material  angle 
T  *  3x3  matrix 


None 

CALL  TMAT  (THETA, T) 
ELEM3,  ELEM4 
None 

(420  bytes)  105  words 
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1.  Subroutine  Name:  TMAT 

2.  Purpose:  Generate  transformation  matrix  tor  triangular  plate 

orthotropic  material  angle. 

2  2 

e  Sin  s  SineCose 

2  2 

s  Cos  0  -SineCose 

2  2 

SineCose  2  SineCose  Cos  e-Sin  e 

J 


4. 

Input  Arguments: 

THETA  =  material  angle 

5. 

Output  Arguments: 

T  =  3x3  matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  TMAT  (THETA, T) 

8. 

Subroutine  User: 

ELEM3,  ELEM4 

9. 

Subroutines  Used: 

None 

10. 

Storage  Required: 

(420  bytes)  105  words 
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/ 


1.  Subroutine  Name: 

2.  Purpose: 


WRITEL 

Tests  character  of  element  connection  card  and 
writes  element  information  into  file. 


3.  Equations  and 

Procedures: 

4.  Input  Arguments: 


If  and 
and  are  0, 

L7 

ICR 

M 

CL »  C2  9  Cp  c 


are  0,  C2  and  are  stored.  If  C2 
and  Cj  are  stored. 

Output  file  number 
i.rid  point  number 

Position  of  A  array  to  be  restored. 
^  Input  codes 


5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Subroutine  User: 

9.  Surboutines  Used: 
10.  Storage  Required: 


IGR  and  A  array  are  stored  on  file  L7. 
None 

CALL  WRITEL  (Cj  .C^.C^C^.M.IGR.L?) 

ZZ 

None 

(578  bytes)  145  words 
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1.  Subroutine  Name: 


XTRAK 


2.  Purpose: 


Interpret  degree  of  freedom  information 


3.  Equations  and 
Procedures: 


Interprets  NWORD  and  breaks  this  down  into  6 
individual  components.  These  components  are 
then  stored  in  KWORD  array. 


4.  Input  Arguments: 


NWORD 

NP 


No.  components  input  word 
Control  word 


5.  Output  Arguments: 


KWORD 

NUM 


Output  data  array 
Total  no.  DOF  recognized 


6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  XTRAK  (NWORD, KWORD, NUM, NP) 

8. 

Subroutine  User: 

SPCSUB.ZZ 

9. 

Subroutines  Used: 

None 

10. 

Storage  Required: 

(714  bytes)  179  words 
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1.  Subroutine  Name: 


2.  Purpose: 

3.  Equations  and 

Procedures: 

4.  Input  Arguments: 

5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Subroutine  User: 

9.  Surbroutines  Used: 


10.  Storage  Required: 


ZZ 


Generates  OPTIM  data  which  is  input  by  NASTRAN 
format  input  cards. 

Each  card  is  read  and  interpreted  based  on  content 
and  use  in  the  OPTIM  program. 

L5  *  Input  file  tape  number 
L7  *=  Output  file  tape  number 

All  of  the  grid  point,  boundary  condition,  element, 
material  property,  load,  constraint  and  buckling 
information  arrays  needed  by  OPTIM. 


CALL  ZZ  (NAST , NOPT ,MATN0 , NPROP , NBUCK , NNODES , NOID , 

REF , OPDVIR, LBOUND , COOR , AMAT , MID , EYEC , LINKS , NEL , NGR , 
NMT , NICON8 , NLINKS , NOGCON , NP , FA , NAP , GA , NFO , NMO , NLOADS , 
NSPCS, IGRID, SPCINF.GCOND, FORCE3, ANGLE, FORCE1 , F0RCE2 , 
M0MNT1 , M0MNT2 , OPLOAD , CROUP , L5 , L7 ) . 

LINK1 

FOMO 

INSPC 

SPCSUB 

WRITEL 

XTRAK 


(18018  bytes)  4505  words  ' 

i 
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A. 4. 2  LINK2  SUBPROGRAMS 


Each  subprogram  of  LINK2  is  described  preceded  by 
glossary  of  argument  definitions. 
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A. 4. 2.1  BRIEF  SUBROUTINE  DESCRIPTIONS 


106 


MFSD  Factor  a  given  symmetric  positive  matrix 


Subroutine 

Name  Purpose 


S316  Global  flexibility  matrix 


Subroutine 


A. 4. 2. 2  GLOSSARY  OF  ARGUMENT 
DEFINITIONS  FOR  I.INK2  ROUTINES 


A 

ALAMBD 

AMIN 

AREA 

COND 

IAREA 

IC 

ICC 

ICHECK 

I  DEL 

IDYN 

IMUAL 

IMUSL 

IMUDL 

IMUXL 


Vector  containing  the  current  design  parameters  (length  =  NE) 

Vector  for  storage  of  Lagrange  multipliers  (length  -  NE  + 
NSNL+NDCNL+NXNI.+NWNL ) 

Vector  for  storage  of  the  minimum  permissible  design  vari¬ 
ables  (length  =  NE) 

Vector  for  storage  of  surface  areas  of  quad,  triangle  and  shear 
panel  type  elements  (length  =  NE) 

IA  i  ^  “  A  I 

- DS-H  <C0ND) 

Aold  1 

Vector  identifying  unknown  design  variables  and  redundants. 

This  vector  is  used  for  row  or  column  labels  when  assembling 
the  partials  into  the  coefficient  matrix  for  solution  of  a  set 
of  linear  equations  (length  =  NE+NXNL) ;  IAREA(I)  =  0;  Variable 
known  -  IAREA (I) >0;  Variable  unknown  and  Value  =  l.abel 

Dummy  vector  containing  input  matrix  column  labels 

Iteration  counter 

Vector  for  intermediate  output,  identifying  the  convergence 
status  of  unknowns,  0  =  has  not  converged,  1  =  has  converged 

2  =  not  an  unknown 

(Length=N  E+N  SNL+NDCNL+NXNL+N WNL ) 

Pointer  to  starting  position  in  vector  W  at  which  the  dis¬ 
placements  are  stored 

Switch  for  inclusion  of  dynamic  constraints 
0  -  No  dynamic  constraints  in  problem 

1  -  Yes  dynamic  constraints  in  problem,  but  not  yet  included 
2  -  Yes  dynamic  constraints  in  problem 

Pointer  to  starting  position  in  vector  IMU  for  storage  of 
design  constraint  "active/not  active"  switches 

Pointer  to  starting  position  in  vector  IMU  for  storage  of 
stress  constraint  "active/not  active"  switches 

Fointer  to  starting  position  in  vector  IMU  for  storage  of 
displacement  constraint  "active/not  active"  switches 

Pointer  to  starting  position  in  vector  IMU  for  storage  of 
redundants  "active/not  active"  switches 


IMUWL  Pointer  to  starting  position  in  vector  IMU  for  storage  of 

frequency  constraint  "active/not  active"  switches 

IMU  Vector  identifying  "active"  constraints.  This  vector  is  used 

for  column  or  row  labels  when  assembling  the  partials  into  the 
coefficient  matrix  for  solution  of  a  set  of  linear  equations 
(length=NE+NSNL+NDCNL+NXNL+NWNL) 

IMU(I)  =  0;  constraint  is  not  active 

IMU(T)  >  0;  constraint  active,  value  of  IMU(I)  =  Label 

IRST  Control  for  execution  of  statics,  dynamics  or  optimization 

ISTRES  Stress-ratio  method  control  variable 

IR  Dummy  vector  containing  input  matrix  row  labels 

ISI  Pointer  to  starting  position  in  vector  "W"  at  which  the  S  values 

are  stored 

1YI  Pointer  to  starting  position  in  vector  "W"  at  which  the  Y  values 

are  stored 

IX  Pointer  to  starting  position  in  vector  "W"  at  which  the 

redundants  are  stored 

IPK  Pointer  to  starting  position  in  vector  "WDYN"  at  which  the 

PR  values  are  stored 

IXK  Pointer  to  starting  position  in  vector  "WDYN"  at  which  the 

XK  values  are  stored 

IZ  Pointer  to  starting  position  in  vector  "WDYN"  at  which  the  mass 

matrix  is  stored 

KL  Vector  containing  unreduced  degrees  of  freedom  numbers  for 

reactions  (length  *  NBOU) 

LCOL  Initial  column  number  of  matrix  waray  for  assembly  of  a  group 

of  coefficients.  Used  by  routine  assem  to  load  the  column 
label  vector 

LROW  Initial  row  number  of  matrix  waray  for  assembly  of  a  group  of 

coefficients.  Used  by  routine  assem  to  load  the  row  label 
vector 

LOW  Vector  containing  lower  displacement  limits  (length  =  NDC) 

MAXIT  Maximum  permissible  iterations  to  be  performed  in  all  con¬ 

vergence  loops 

NASTAR  The  number  of  active  ('on')  design  variable  (minimum  size) 

constraints 
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NBOU 


NC 

NCOL 

NCOL1 

ND 

NDC 

NDCNL 

NI)N 

NDNNL 

NDSTAR 

NDT 

NDTNX 

NF. 

NLL 

NODE 

NR 

NROW 

NROW1 

NSE 

NSENL 

NSET 

NSNL 


Length  of  vector  KL 

Dummy  argument  indicating  total  number  of  columns  in  the 
coefficient  matrix. 

Number  of  columns  in  the  A-matrix  from  the  structure  cutter 
rout  ine  NCOL  =  NSE+NBOU+1 

Number  of  columns  in  matrix  waray  into  which  a  group  of 
coef f ir ients  is  to  be  assembled.  Used  by  ASSEM  routine  to 
load  column  label  wctoi 

Vector  containing  degree  of  freedom  for  displacement  constraints 
(length  =  NDC) 

Number  of  displacement  constraints 

Number  of  displacement  constraint  variables  -  NDC  x  number  of 
loads 

Number  of  reduced  degrees  of  freedom 
NDN  x  number  of  loads 

The  number  of  active  ('on')  displacement  constraints 

Total  number  of  degrees  of  freedom  in  problem 

Total  number  of  columns  in  structure  cutter  matrix 

Number  of  elements  in  structure 

Number  of  load  conditions 

Total  number  of  nodes  in  structure 

Dummy  argument  containing  total  number  of  rows  in  the 
coefficient  matrix 

Number  of  rows  in  the  A-matrix  from  structure  cutter  routine 

Number  of  rows  in  matrix  WARAY  into  which  a  group  of  coeffi¬ 
cient  is  to  be  assembled.  Used  by  ASSEM  routine  to  load 
row  label  vector 

Number  of  element  forces  in  the  structure 
NSE  x  number  of  loads 

Dummy  vector  for  identifying  the  active  constraints  as 
determined  by  the  FSD  routine  (length  =  NSNL) 

Number  of  stress  constraint  variables  =  NE  x  number  of  loads 
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NS STAR 


NT RAN S 


NTYPE 

NW 

NWNL 

NWSTAR 

NX 

NXNL 

NXSTAR 

OMEGA 

P 

PHI 
PRD 
PRR 
PRS 
PRI 
PS  I 
S 

SIG 

UP 

W 

WARAY 

WDYN 


The  number  of  active  ('on')  stress  constraints 

Control  to  transpose  output  of  a  partial  computation  routine 
prior  to  assembly  in  matrix  waray. 

1  =  no  transpose 

2  =  transpose 

-N  =  assemble  in  lower  triangle 
Vector  identifying  element  type  (length  *  NE) 

Number  of  modes 

Number  of  frequency  constraint  variables 
Number  of  active  ('on')  frequency  constraints 
Number  ol  rcdundants 

Number  of  redundant  constraint  variables  -  NX  x  number  of  loads 
Number  of  active  ('on')  redundant  constraints 
Assembled  structural  matrix 

Vector  containing  static  loading  of  structure  (length  =  NBOU) 
Assembled  structural  matrix 

Control  for  printing  of  displacements  every  iteration 
Control  for  printing  of  reactions  every  iteration 
Control  for  printing  of  stress  every  iteration 
Control  for  printing  of  intermediate  data 
Assembled  structural  matrix 

Vector  containing  parameters  describing  stress  state  of 
elements  (length  =  NSENL) 

Vector  containing  stress  limits  (length  =  NE) 

Vector  containing  upper  displacement  limits  (length  =  NDC) 

Vector  containing  all  redundants,  displacements,  S  and  Y 
values,  (length  =  NXNL+NDNNL+NSENL+NSNL) 

Dynamic  storage.  Vector  for  storage  of  computed  partlals  and 
any  work  storage  required  by  Phase  4  subroutines. 

Vector  for  storing  dynamic  variables  DEL,  P^,  and  Z. 

(Length  =  2xNDN*NWNL+NX*NWNL+NE*NDN ) 
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ws 


Frequency  limit 

Vector  containing  element  normalized  weight  (length  =  NE) 
Vector  containing  redundants  (length  =  NXNL) 


VTT 
X 

XC  Vector  containing  the  X-coordinates  of  centroid  of  quad 

elements  (length  =  NE) 

Y  Vector  containing  element  forces  (length  =  NSNL) 

YC  Vector  containing  the  Y-coordinates  of  centroid  of  quad 

elements  (length  =  NE) 
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1. 

Subroutine  Name: 

AA 

2. 

Purpose: 

Structure  cutter  control  stage 

3. 

Equations  and 
Procedures : 

The  de  bug  print  control  determines  whether  the 
matrix  to  be  processed  is  printed  before  and  after 
the  structure  cutter  routine  MFGR. 

4. 

Input  Arguments: 

A , M , N , LL , AO , U , WO kK ,  WOHKt ,BO, IROW, ICOL 

These  arguments  reserve  storage  for  internal  cal¬ 
culation  by  structure  cutter  routines  AONE  and 

MFGR. 

5. 

Output  Arguments: 

A, IROW, ICOL  are  returned  after  the  structure 
cutter  process  and  written  on  FILE  112  by  the 

TAPEII  routine. 

6. 

Error  Returns: 

if  the  rank  of  the  matrix  does  not  equal  the  no. 
of  rows  an  error  message  is  printed  and  the  routine 
returns  to  TAPEII. 

7. 

Calling  Sequence: 

Call  AA(A,M,N,CC,A0,U,W0RK,W0RK1,B0, IROW, ICOL) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1314  bytes)  425  words 

12. 

Subroutine  User: 

TAPEII 

13. 

Subroutine  Required: 

SUBA , MFGR , MRPRNT 
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Subroutine  Name: 

ASSEM 

Purpose: 

Compute  row  and  column  vectors  for  matrix 
assembly. 

Equations  and 
Procedures : 

a)  Load  row  vector  with  row  numbers  (sequential) 
for  assembly  of  an  input  matrix. 

b)  Load  column  vector  with  column  numbers 
(sequential)  for  assembly  of  an  input  matrix. 

Input  Arguments: 

LR0W,LC0L,NR0W1 ,NC0L1 , 1R , 1C 

See  Glossary  for  definitions. 

Output  Arguments: 

IR.IC 

See  Glossary  for  definitions 

Error  Returns: 

None 

Calling  Sequence: 

Call  ASSEM  (LROW.LCOL.NROWl.NCOLl, IR.IC) 

Input  Tapes: 

None 

Output  Tapes: 

None 

Scratch  Tapes: 

None 

Storage  Required: 

<442  BYTES)  111  words 

Subroutine  User: 

OPTFR , S451 , S461 , EQSOL 

Subroutine  Required:  None 


1.  Subroutine  Name:  BASIC 

2.  Purpose:  Statics  analysis  computes  X,S,Y  and  Delta 

3.  Equations  and  a)  Compute  working  storage  requirements 

Procedures:  b)  Compute  PHI  inverse,  PSI  and  Omega.  Call  S311. 

c)  Compute  redundants  if  requested.  Call  S314. 

d)  Compute  displacements  (call  S315)  and  print 
if  requested. 

e)  Compute  work  storage  requirements. 

f)  Read  row  labels,  col  labels  and  A-matrix 
from  structure  cutter  output. 

g)  Compute  S(I).  Call  S231. 

h)  Compute  stress  and  1/(1)  for  each  element. 

Call  S433. 

i)  Compute  reactions  if  requested.  Call  S231. 

j)  Print  reactions  if  requested. 

4.  Input  Arguments:  PRD,PRS,PRR,NTYPE,AM1N,WT,AREA,XC,YC,A,KI.,P,NE, 

NDN , NX , NBOU , NSE , NODE , NDT , NDTNX , W 
See  Glossary  for  definitions. 

5.  Output  Arguments:  NCOL,W 

See  Glossary  for  definitions. 

6.  Error  Returns:  None 

7.  Calling  Sequence:  Call  BASIC  (PRD, PRS,PRR,NTYPE, AMIN, WT, AREA, XC, 

YC,A,KL,P,NE, NDN, NX, NBOU, NSE, NODE, NDT, NDTNX, HCOL, 

8.  Input  Tapes:  NTAPE,  UNIT  12  (output  of  structure  cutter) 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (3770  BYTES)  943  words 

12.  Subroutine  User:  FSD,0PTFR,S451,S461 

13.  Subroutine  Required:  S314,S315,S231,S433,S311 
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1.  Subroutine  Name: 


BASICD 


2.  Purpose: 


3.  Equations  and 
Procedures : 


4 .  Input  Arguments : 


5 .  Output  Arguments : 


To  perform  a  dynamic  analysis  and  return 
data  that  will  be  used  in  a  dynamically 
constrained  optimization  problem. 

S222  is  called  and  from  Tape  113  and  the 
element  design  parameters,  matrices  <t>,  41,  ft 
and  M  are  formed.  Routine  S316  is  then  called 
to  form  the  flexibility  matrix.  Routine  S321 
then  replaces  the  flexibility  matrix  by  the 
eigenvalue  matrix.  The  EIO  routine  then  solves 
for  the  eigenvalues  and  eigenvectors.  The 
inertial  load  modes  PK  are  computed  in  routine, 

5323.  The  displacement  modes  are  computed  by 

5324.  The  mode  shapes  are  normalized  on  the 
max  component  by  MAXVAL  and  printed.  The 
redundant  modes  XK  are  computed  by  S314.  If 
IPRT*=1  (set  at  the  beginning  of  the  routine) 
the  routine  returns,  otherwise  the  inertial 
stress  and  inertial  reaction  calculation  and 
print  is  done  by  S231  and  S433. 

PRD, PRS, PRR.NTYPE, AREA, XC.YC, A, KL.NE, NDN, NX, 
NBOU.NSE, NODE , NDT , NDTNX , NCOL , NW 
See  Glossary  for  definition. 

U  -  working  storage 

NOI  -  maximum  number  of  iterations 

CRIT  -  convergence  criteria 

In  the  W1  array  DEL(NDN.NW) ,PK(NDN,NW) , 
XK(NX,NW) ,  Z(NE,NDN)  are  contained. 


6.  Error  Returns:  None 


7.  Calling  Sequence: 

8.  Input  Tapes: 

9.  Output  Tapes: 

10.  Scratch  Tapes: 

11.  Storage  Required: 

12.  Subroutine  User: 


Call  BASICD (PRD, PRS, PRR.NTYPE, AREA, XC.YC, A, KL,NE, 
NDN , NX , NBOU , NSE , NODE , NDT , NDTHX , NCOL , W1 , W , NW , 
NOIT.CRIT) 
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None 


None 

(4526  BYTES)  1142  words 
OPTFR 


13.  Subroutine  Required: 


EIG,  S222.S231, S314 ,S316,S321 ,S323,S324 ,S433, 
MAXVAL .MRPRNT 


1. 


Subroutine  Name: 


DYNERR 


2. 

Purpose : 

Compares  storage  required  against  storage 
allotted. 

3. 

Equations  and 
Procedures : 

a)  Compute  total  storage  required. 

b)  Call  EXIT  if  required  storage  greater 
that  dimensioned  in  LINK2. 

4. 

Input  Arguments: 

NEND  -  Previous  storage  requirement 

NENUl  -  Additional  storage  requirement 

5. 

Output  Arguments : 

I  ERR 

6. 

Error  Returns: 

Writes  dynamic  storage  allocation  error... 
call  EXIT. 

7. 

Calling  Sequence: 

Call  DYNERR  (NFND.NENDl , IERR) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(538  BYTES)  135  words 

12. 

Subroutine  User: 

LINK2 

13. 

Subroutine  Required: 

EXIT 
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1.  Subroutine  Name:  EIG 

2.  Purpose:  Generates  eigenvalues  and  eigenvectors  by 

power  method. 

3.  Equations  and 

Procedures : 


11; i s  equation  is  repeated  until  and 

fall  within  the  criteria  and  X  and  X 

n  n+1 

converge . 


The  matrix  is  deflated  removing  the  higher  eigenvalue  and  then  the  previous 
iteration  is  repeated  to  get  the  next  value  and  vector. 

4.  Input  Arguments:  A  -  Eigenmatrix  less  than  OR  «=  (N94.N94) 

N  -  Order  of  matrix  A  less  than  OR  «=  N94 
IPRINT  -  Print  iteration  control 
NEIG  -  No.  of  eigenvalues  requested  less  than 
OR=N20 

GUESS  -  Guess  vectors  (N94,N20)  usually  1.0 
CRIT  -  If  *=0  then  .001  relative  error  criteria 

NOIT  -  If  =0  then  54  maximum  no.  of  iterations 
NEL  -  Selects  1ST  eigenvalue  requested  usually  =1 

N94  -  Order  of  a  matrix 

N20  -  No,  of  eigenvalues 

16  -  Print  output  unit 

XIN.XI .XIMIN.XIP.XIMINP.XINP-storage  for  successive 

eigenvector  iterations 

ICOUNT  -  Total  No.  of  eigenvalues  found 
VECTOR  =  Eigenvectors  ROOT  *  Eigenvalues 

NERR  «  0  NO  ERROR  =  1  eigencols 
Do  not  converge  ■  2  Eigenrows  do  not 
converge  *  3  valves  do  not  converge 

Cal 1  EIG (A ,N .IPRINT , NEIG .GUESS , ROOTS .VECTOR , 

NERR , CRIT ,N0IT , ICOUNT , NEL , N94 , N20 , 16 ,  XIN , 

XI , XIMI N , XIP , XIMINP , XINP) 

None 

None 

None 

(4438  BYTES)  1110  words 
BASICD 

13.  Subroutine  Required:  None 


5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Input  Tapes: 

9.  Output  Tapes: 

10.  Scratch  Tapes; 

11.  Storage  Required: 

12.  Subroutine  User: 
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1.  Subroutine  Name:  EQSOL 

2.  Purpose:  Compute  Lagrange  multipliers  for  "on" 

constraints . 

a)  Compute  matrix  dimensions  and  initialize 
arrays, WARAY  &  IAREA. 

b)  Compute  and  assemble  the  partials  C*A/A, 
G*S/A,G*S/X,G*0/A,G*D/X,GX/A,GX/X  and 
G*W/A  in  array  WARY  by  calling  routines 
ASSEM,S432,S435,S4310,S4323  and  S4325. 

c)  Solve  linear  equations  for  Lagrange 
multipliers.  Call  SIMQ. 

d)  Store  multipliers  in  vector  ALAMBD. 

e)  Reset  vector  IAREA. 

4.  Input  Arguments:  W, WARAY , IAREA, IMU ,NASTAR,NSSTAR,NDSTAR, NX, 

NWSTAR . IMUAL , IMUSL , IMUDL , IMU XL , IMUWL , AMIN , A , NE , 
IYI, ISI, NTYPE, SIC, AREA, XC.YC.NSE, NOT, NCOL.NOTNX, 
NBOU , NDC , P , NDN , UP , LOW , ND , IX , NW , WT , NSNL , NDCNL , 
NWNL , WD YN , PRI , WS , I PK , I XK , I Z 
See  Glossary  for  definitions. 

5.  Output  Arguments:  ALAMBD  -  See  Glossary 

6.  Error  Returns:  Write  EQSOL  matrix  is  singular . .STOP. 

7.  Calling  Sequence:  Call  EQSOL ( ALAMBD, W, WARAY, IAREA, IMU, NASTAR, 

NSSTAR , NDSTAR , NX , NWSTAR , IMUAL , IMUSL , IMUDL , 

IMUXL , IMUWL , AMIN , A, NE , IYI , ISI .NTYPE , SIG .AREA , 

XC , YC , NS  E , NDT , NCOL , NDTNX , NBOU , NDC , P , NDN , UP .LOW , 
ND , IX , NW , WT , IERR , NSNL , NDCNL , NXNL , NWNL ,WDYN , PRI , 
WS,IPK,IXK,IZ,ARAY) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (4516  BYTES)  1137  words 

12.  Subroutine  User:  OPTFR 

13.  Subroutine  Required:  SIMQ,S432,S435,ASSEM,S4310,S4323,S4325,MRPRNT, 

SHIFTA.S4310B 


3.  Equations  and 
Procedures: 
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1. 

Subroutine  Name: 

EXIT 

2. 

Purpose: 

To  provide  one 

3. 

Equations  and 
Procedures : 

This  routine  re 
monitor.  (Same 

4. 

Input  Arguments: 

None 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  EXIT 

8. 

Input  Tapes: 

None 

9, 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

12. 

Subroutine  User: 

DYNERR 

13. 

Subroutine  Required: 

None 

1. 


Subroutine  Name: 


FSD 


2 .  Purpose : 


Fully  stressed  design  -  loops  on  BASIC 
routine . 


3. 

Equations  and 

a) 

Initilize  multiple  load  variables  and 

Procedures : 

iteration  counter 

b) 

Convergence  loop 

1.  Perform  a  statics  analysis  based  on 
new  design  (compute  Y(I)  (Call  BASIC) 

2.  for  each  design  variable,  compute  a 
new  design  for  each  load  condition 
and  store  the  largest  value  and  the 
load  condition  number  ANEN=Y( I ) /SIG(1 ) 

3.  Set  design  variable  co  minimum  value  if 
smaller  than  minimum. 

4.  Check  for  convergence  and  store  new 
design. 

5.  Reset  the  stress  constraint  vector 
based  on  load  condition  number. 

6.  If  not  all  design  variables  converged 
and  iteration  number  is  less  than 
maximum,  repeat  convergence  loop. 

c)  Reset  stress  constraint  vector  to  0  for 

any  design  which  is  minimum. 

d)  Compute  structure  weight. 


4.  Input  Arguments: 


5.  Output  Arguments: 


PRD , PRS , PRR , NTYPE , AMIN , WT , AREA ,XC,YC,A,SIG,KL,P, 
NE , NDN , NX , NBOU , NSE , NODE , NDT , NDTNX , W , NSET , COND , 
MAXIT  (See  Glossary  for  definitions) 

A,W  -  See  Glossary  for  definitions 
NSET  -  Vector  containing  active  constraint 
identification 


6.  Error  Returns:  None 

7.  Calling  Sequence:  Call  FSD (PRD, PRS, PRR, NTYPE, AMIN, WT, AREA, XC.YC, A, 

SIG,KL,P,NE, NDN, NX, NBOU, NSE, NODE, NDT, NDTNX, W, 
NSET, COND, MAXIT) 


8. 


9. 


10. 


11. 


12. 


13. 


Input  Tapes: 

Output  Tapes: 

Scratch  Tapes : 
Storage  Required: 
Subroutine  User. 
Subroutine  Required: 


None 

None 

None 


(2766  BYTES)  692  words 
OPTFR 

BASIC, TRPRNT 
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1. 

Subroutine  Name: 

GELS 

2. 

Purpose: 

To  solve  a  system  of  simultaneous  linear  equations 
with  symmetric  coefficient  matrix. 

3. 

Equations  and 
Procedures : 

Solution  is  obtained  by  means  of  Gauss-elimination 
with  pivoting  in  main  diagonal,  in  order  to  pre¬ 
serve  symmetry  in  remaining  coefficient  matrix. 

a)  Search  for  greatest  main  diagonal  element. 

b)  Elimination  loop. 

1)  Test  usefulness  of  symmetric  algorithm 

2)  Pivot,  row  reduction,  and  row  interchange 
in  right-hand  side  matrix 

3)  Row  and  column  interchange  and  pivot  row 
reduction  in  coefficient  matrix 

4)  Save  column  interchange  information 

5)  Element  reduction  and  search  for  next 
pivot 

c)  Back  substitution  and  row  interchange 

4. 

Input  Arguments: 

R  -  M  by  N  right-hand  side  matrix 

A  -  Upper  triangular  part  of  the  symmetric 

M  by  N  coefficient  matrix 

M  -  The  number  of  equations 

N  -  The  number  of  right-hand  side  vectors 

EPS  -  Relative  tolerance  for  test  on  loss  of 
significance 

AUX  -  M-l  auxiliary  storage  vector 

5. 

Output  Arguments: 

R  -  M  by  N  matrix  containing  solution  to 

equations 

6. 

Error  Returns: 

IER 

7. 

Calling  Sequence: 

Call  GELS (R, A,M,N,EPS , IER , AUX) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

1730  bytes  (433  words) 

12. 

Subroutine  User: 

S461 

13. 

Subroutine  Required: 

None  j 

■  i'-S".  - 
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1. 

Subroutine  Name 

INTPR 

2. 

Purpose: 

Prints  linear  programming  tableau  for  routine 
ITER 

3. 

Equations  and 

a)  Print  vector  LABC 

Procedures: 

b)  Print  vector  LABR  and  ARRAY  "A" 

4. 

Input  Arguments : 

A  -  Array  of  dimension  MM*N  to  be  printed 

LABR  -  Vector  of  length  M  to  be  printed 
(row  labels) 

LABC  -  Vector  of  length  N  to  be  printed  (column 
labels) 

MM  -  First  dimension  of  "A"  in  calling  program 

M  -  Number  of  rows  in  ARRAY  "A" 

N  -  Number  of  columns  in  ARRAY  "A" 

IT  -  Number  of  iterations 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  INTPR(A, LABR, LABC, MM, M,N, IT) 

8. 

Input  Tapes: 

None 

9. 

Output  File: 

Unit  101  (printer) 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(656  BYTES)  164  words 

12. 

Subroutine  User: 

ITER 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Name: 


ITER 


2.  Purpose: 


3.  Equations  and 
Procedures : 


4.  Input  Arguments: 


Solve  linear  equations  with  provision  for  the 
automatic  treatment  of  any  or  all  of  the  follow¬ 
ing:  Equality  constraints,  greater  than  or 

equal  constraints,  less  than  or  equal  constraints, 
positive  variables,  free  variables  and  either 
maximization  or  minimization  of  the  objective 
function. 


Method 

The  method  is  a  primal-dual  algorithm  (some¬ 
times  referred  to  as  the  'criss-cross'  method). 

The  condensed  tableau  is  pivoted  in  the  follow¬ 
ing  manner.  First  artificial  variables  (which 
correspond  only  to  equality  constraints)  are 
pivoted  out  of  the  basis  and  free  variables 
(if  any)  are  pivoted  in.  If  there  are  still 
artificial  variables  in  the  basis  after  all 
free  variables  have  been  pivoted  in,  then 
the  remaining  artificial  variables  are  pivoted 
out  and  positive  variables  are  pivoted  in.  If 
there  are  more  free  variables  than  artificial, 
then  the  remaining  non-basic  free  variables  are 
next  pivoted  in  and  positive  variables  out. 

Now  alternate  dual  and  primal  iterations  are 
taken  until  the  solution  becomes  primal  or 
dual  feasible.  From  this  point  on  only  pri¬ 
mal  or  only  dual  iterations  are  taken  until  the 
optimal  solution  is  attained. 

In  the  event  the  solution  is  primal  feasible 
but  not  dual  feasible  and  no  primal  iteration 
may  be  taken,  the  solution  is  primal  unbounded. 
Similarly,  if  the  solution  is  dual  feasible  but 
not  primal  feasible  and  no  dual  iteration  may 
be  taken,  the  solution  is  dual  unbounded.  If 
the  solution  is  neither  primal  nor  dual  feasible, 
and  neither  a  primal  nor  a  dual  iteration  may 
be  taken,  then  the  solution  is  either  unbounded 
or  infeasible  in  either  the  primal  or  the  dual 
sense . 

For  reference  see,  'Linear  and  Integer  Program¬ 
ming,'  by  Stanley  Zionts,  Prentice-Hall,  1974. 

A  Array  of  dimension  MM  X  N 

A(l, 2)-A(l ,N)  Contain  the  coefficients 
of  the  cost  vector. 

A( 2 , 1) -A(M, 1)  Contain  the  RHS's  of  the 
constraint  equations. 

Rows  2-M  contain,  starting  in  column  2, 

the  coefficients  of  the  constraint  equations. 

A(l,l)  is  irrelevant 
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5.  Output  Arguments: 


6.  Error  Returns: 

7.  Calling  Sequences: 


8.  Input  Tapes: 

9.  Output  Tapes: 


LABR  Vector  of  length  M  containing  row  labels 
LABR(I)  indicates  the  type  of  constraint 
corresponding  to  row  I  as  follows: 

-1=LE,  0=Equality,  i=GE 
LABR(l)  is  irrelevant 

LABC  Vector  of  length  N  containing  column  labels 
LABC(I)  indicates  the  type  of  activity 
variable  corresponding  to  column  I  of 
array  A  as  follows: 

1-positive  variable,  0=free  variable 
LA3C(1)  is  irrelevant 

MM  First  dimension  of  array  "A"  in  calling 

program.  If  "A"  is  singly  subscripted  MM-'-' 

M  Number  of  constraints  +.1 

N  Number  of  activity  variables  +1 

NMAX  Code  to  maximize  (1),  or  minimize  (-1) 

TOL  Minimum  permissible  absolute  value  for 

non-zero  elements  of  array  "A".  Elements 
falling  below  TOL  are  set  to  0.0 

T0L1  Minimum  permissible  absolute  value  for  ar 
element  to  be  used  as  a  pivot 

IPFLAG  Print  code  for  array  "A" 

~3  =  no  printout 

-2  =  print  after  final  iteration 

-1  =  print  before  initial  and  after  final 
iteration 

0  =  print  before  initial  iteration 

1  =  print  before  initial  iteration  and 

after  all  succeeding  iterations 

A  Array  of  dimension  MMX  N 

(A(l,l)  x  MMAX)  is  the  optimal  value  of  the 
cost  function 

A( 2 , 1) -A(M, 1)  contain  corresponding  values 
of  basic  values 

LABR  Vector  of  length  M  containing  row  labels 
LABR(1)=0 

The  value  of  activity  variable  number 
(ABS (LABR(I) )  is  contained  in  A(I,1) 

IR  Number  of  iterations 

KOPT 

Call  ITER  (A, LABR, LABC, TOL, T0L1 ,MM,M,N,MMAX, IPLAG , 
IT, KOPT) 

None 

Noi.e 


127 


10.  Scratch  Tapes: 

11.  Storage  Required: 

12.  Subroutine  User: 

13.  Subroutine  Required: 


None 

2778  bytes  (695  words 
OPTFR 

PIVOT  INTPR 
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1.  Subroutine  Name: 


LINK2 


2.  Purpose:  Main  calling  routine  for  structure  cutter  and 

calculation  phases, 

a)  Initialize  I/O  unit  variables. 

b)  Compute  dynamic  storage  requirements  and 
check  against  maximum  allocated  (call 
DYNERR)  tor  TAPEII  routine. 

c)  Read  phase  2  output  and  execute  structure 
cutter  routines  call  TAPEII. 

d)  Compute  dynamic  storage  requirements  for 
storage  of  phase  3  output  and  check  against 
maximum  aJ located, call  DYNERR. 

e)  Read  phase  3  output  -  from  file  111. 

f)  Initialize  optimization  variables. 

g)  Compute  dynamic  storage  rquirements  for 
phase  4  and  check  against  maximum  allo¬ 
cated,  call  DYNERR. 

h)  Execute  phase  4  computations  (statics  or 
dynamics  or  optimization),  call  OPTFR. 


4. 

Input  Arguments: 

W  dynamic  storage 

NWORK  maximum  size  of  dynamic  storage 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  LINK2  (W, NWORK, AAA) 

8. 

Input  Tapes: 

NSS2  Unit  1 

Ill  Unit  U 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(4468  bytes)  1117  words 

12. 

Subroutine  User: 

MAIN 

13. 

Subroutine  Required: 

DYNERR,  TAPEH,  OPTFR 

3.  Equations  and 
Procedures: 


1.  Subroutine  Name;  LOC 

2.  Purpose:  Compute  a  vector  subscript  for  an  element  in 

a  matrix  of  specified  storage  mode. 

3.  Equations  and 

Procedures:  NS**0  Subscript  is  computed  for  a  matrix  with 

N*M  elements  in  storage  (general  matrix) . 

NS*1  Subscript  is  computed  for  a  matrix  with 
N*(N+l)/2  in  storage.  (Upper  triangle 
of  symmetric  matrix.) 

NS=2  Subscript  is  computed  for  a  matrix  with 
N  element  in  storage  (diagonal  elements 
of  a  diagonal  matrix.) 

4.  Input  Arguments:  I  Row  number  of  element 

J  Column  number  of  element 

M  Number  of  runs  in  matrix 

N  Number  of  columns  in  matrix 

NS  One  digit  number  for  storage  model  of  matrix 

0  -  General 

1  -  Symmetric 

2  -  Diagonal 


5. 

Output  Arguments: 

IR-Resultant  vector  subscript 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  LOC  (I,J,IR,N,M,NS) 

8, 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(492  bytes)  123  words 

12. 

Subroutine  User: 

MPRD.MSTR 

13. 

Subroutine  Required: 

None 

1. 

Subroutine  Name: 

MAXVAL 

2. 

Purpose: 

General  maximum  absolute  value 

of  Eigenvector. 

3. 

Equations  and 

Procedures : 

The  vector  is  searched  to  find 
component . 

maximum  absolute 

4. 

Input  Arguments: 

DISP=Vector 

NON=Length  of  vector 

5. 

Output  Arguments: 

VAL=Actual  algebraic  value  of  maximum  component 

6. 

Error  Returns : 

None 

7. 

Calling  Sequence: 

Call  maxual  (DISP,NDN,VAL) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(324  bytes)  81  words 

12. 

Subroutine  User: 

BASICD 

13. 

Subroutine  Required: 

None 
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1. 

Subroutine  Name: 

MFGR 

2. 

Purpose: 

3. 

Equations  and 
Procedures: 

Gaussian  elimination  technique  is  used  for 
calculation  of  the  triangular  factors  of  a 
given  matrix. 

Complete  pivoting  is  built  in. 

4. 

Input  Arguments 

A  Given  matrix 

M  No.  of  rows  of  A 

N  No.  of  cols,  of  A 

EPS  Test  value  of  zero  affected  by  roundoff 

noise 

5. 

Output  Arguments 

A, IRANK, IROW, ICOL 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  MFGR(A,M,N,EFS, IRANK, IROW, 1C0L) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1974  bytes)  494  words 

12. 

Subroutine  User: 

AA 

13. 

Subroutine  Required: 

None 
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1. 

Subroutine  Name: 

MFSD 

2. 

Purpose: 

Factor  a 

given  symmetric  positive  definite  matrix 

3. 

Equations  and 
Procedures: 

Solution 

obtained  using  the  square  root  method 

of  Cholesky.  The  given  matrix  is  represented 
as  product  of  two  triangular  matrices,  where 
the  left  hand  factor  is  the  transpose  of  the 
right  hand  factor. 

4.  Input  Arguments:  A  -Upper  triangular  part  of  the  given  symmetric 

positive  definite  N  by  N  coefficient  matrix. 
N  -The  number  of  rows  (columns)  in  given  matrix 
EPS-An  input  constant  which  is  used  as  relative 
tolerance  for  test  on  Joss  of  significance. 

5.  Output  Arguments:  A  -Contains  resultant  upper  triangular  matrix 

IER-Error  return  variable 


6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Input  Files: 

9.  Output  Files: 

10.  Scratch  Files: 

11.  Storage  Required: 

12.  Subroutine  User: 


If  input  parameter  N  is  wrong  or  some  radicand 
is  non-positive  then  argument  IER  is  set  to 
-1  and  return  is  made  to  calling  routine. 

Call  MFSD  (A,N,EPS,IER) 

None 

None 

None 

(732  bytes)  183  words 
SINV 


13.  Subroutine  Required:  None 
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1. 

Subroutine  Name: 

MINV 

2. 

Purpose: 

Invert  a  matrix 

3. 

Equations  and 

Procedures: 

The  standard  Causs-Jordan  method  is 
The  determinant  is  also  calculated. 

A  determinant  of  zero  indicates  that 
matrix  is  singular. 

4. 

Input  Arguments: 

A  Matrix 

N  Order  of  matrix 

|  Work  storage  each  of  length  N 

5. 

Output  Arguments: 

A  Resultant  inverse 

D  Resultant  determinant 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  MINV  (A,N,D,L,M) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1874  bytes)  469  words 

12. 

Subroutine  User: 

S311D,  S222,  SUBA,  S311 

13.  Subroutine  Required 


None 


1. 

Subroutine  Name: 

M1NV2 

2. 

Purpose : 

Control  routine  to  insert  a  positive  definite 
symmetric  matrix. 

3. 

Equations  and 
Procedures: 

The  upper  triangle  elements  of  a  general 
matrix  are  used  to  form  a  symmetric  matrix 
(subroutine  MSTR) .  The  matrix  is  inverted 
(subroutine  S^NV).  The  symmetric  matrix  is 
expanded  to  form  a  general  matrix  (subroutine 
MSTR) . 

4. 

Input  Argument  s : 

PHI  Symmetric  matrix  to  be  inverted 

NX  Number  of  rows  (columns)  in  matrix 

W  Work  storage  area 

5. 

Output  Arguments: 

PHI  Resultant  inverted  matrix 

6. 

Error  Returns : 

None 

7. 

Calling  Sequence: 

Call  MINV2  (PHI , NX, D, W, Wl) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(642  bytes)  161  words 

12. 

Subroutine  User: 

S311 

13. 

Subroutine  Required: 

MSTR, S IN V 
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1.  Subroutine  Name: 


MPRD 


2.  Purpose:  Multiply  two  matrixes  to  form  a  resultant 

matrix . 

3.  Equations  and 

Procedures:  The  M  by  L  matrix  B  is  premultiplied  by  the 

N  by  M  matrix  A  and  is  stored  in  the  N  by  L 
matrix  R.  This  is  a  row  into  column  product 
R  Is  always  output  as  a  general  matrix  except 
when  A  and  B  are  both  diagonal  matrices  then 
R  is  output. 

4.  Input  Arguments:  A  First  input  matrix 

B  Second  input  matrix 
N  No.  of  vows  in  A  and  R 

M  No.  of  columns  in  A  and  rows  in  B 

MSA  Storage  mode  of  A 
0-General 

1- Symmetr ic 

2- Diagonal 

MSB  Same  as  MSA  except  for  B 
L  No.  of  columns  in  B  and  R 


5. 

Output  Arguments: 

R  Output  matrix  cannot  have 

as  A  or  B. 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  MPRD(A,B,R,N,M,MSA,MSB,L) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(944  bytes)  238  words 

12. 

Subroutine  User: 

SUBA 

13. 

Subroutine  Required: 

LOC 
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1. 

Subroutine  Name: 

MRPRNT 

2. 

Purpose: 

Matrix  print  routine 

3. 

Equations  and 
Procedures : 

Print  an  NR  x  NC  matrix  (row 

by  row) 

4. 

Input  Arguments: 

NR  number  of  rows 

NC  number  of  columns 

AMTRTX  an  NR  x  NC  matrix  to 
unit  101 

be  output 

5. 

Output  Arguments: 

None 

6. 

Error  Returns : 

None 

7. 

Calling  Sequence: 

Call  MRPRNT  (NR.NC.AMTRIX) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

101  Unit  10 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(502  bytes)  126  words 

12. 

Subroutine  User: 

S461 ,  EQS0L ,  S241 ,  S311,  TAPEII,  S316, 
AA,  OPTFR,  S222,  S316,  S321,  S451,  S314 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Name: 


MSTR 


2.  Purpose:  Change  storage  mode  of  a  matrix 

3.  Equations  and 

Procedures:  The  input  matrix  is  restructured  to  form  an 

output  matrix. 

4.  Input  Arguments:  A  Input  matrix 

N  Number  of  rows  and  columns 

MSA  One  digit  code  for  storage  mode  of 

input  matrix 

0  -  general,  1  -  symmetric,  2  -  diagonal 
MSR  Same  as  MSA  except  for  output  matrix 


5. 

Output  Arguments: 

R  Output  matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  MSTR  (A,R,N,MSA,MSK) 

8. 

Input  Files: 

None 

9. 

Output  Files: 

None 

10. 

Scratch  Files: 

None 

11. 

Storage  Required: 

(542  bytes)  136  words 

12. 

Subroutine  User: 

MINV2 

13. 

Subroutine  Required: 

IOC 

135a 


1.  Subroutine  Name: 


OPTFR 


2.  Purpose: 

3.  Equations  and 

Procedures: 


Control  program  for  optimization 


a)  Initialize  variables  and  arrays. 

b)  If  BASIC  analysis  call  BASIC  or  BAS1CD  and 
return. 

c)  Compute  initial  guess  design  -  use  input 
design  or  approximate  fully  stressed  design 
from  input  design  (call  FSD,  BASIC). 

d)  Initialize  row  and  column  vectors,  and  work 
array . 

e)  Compute  row  and  column  to  store  partials  in 
work  array  (call  assera) . 

f)  Compute  8gA/9A,  9gs/9A,  ggs/9X,  ggD/9A, 

9gD/9X,  9gx/9A,  9gx/9x,  9gw/9A  and  store  in 

work  array  (call  S432,  S435,  S43IO,  S4323,  S432S), 

g)  Initialize  vectors  for  linear  programming 
rout ine . 

h)  Determine  which  constraints  will  be  active  for 
optimum  design  and  if  linear  design  is  feasible 
(call  ITER). 

i)  Initialize  row  and  column  vectors  for  partial 
Newton-Raphson  procedure. 

j)  Perform  the  partial  Newton-Raphson  procedure 
until  design  converges  or  maximum  allowable 
iterations  has  occurred  (cal  S451). 

k)  If  design  has  not  converged  repeat  steps 
D-J-maximum  of  3  iterations  for  linear  pro¬ 
gramming  phase. 

kl)  If  dynamic  constraints  are  to  be  included 
perform  a  basic  analysis  (call  BASICD)  and 
repeat  steps  D-K. 

l)  If  design  did  not  converge  in  linear  program¬ 
ming  phase  compute  a  fully  stressed  design 

.  (call  FSD)  and  initialize  row  and  column  vectors, 
check  for  displacement  or  frequency  constraint 
violations,  initial ize  ALAMBD. 

m)  If  design  converged  in  linear  programming  phase 
or  design  is  fully  stressed  with  no  displace¬ 
ment  or  freq  constraint  violations  then  compute 
lambda  for  all  ’6n"  constraints  (call  EQSOL) . 

n)  Verify  that  all  values  of  lambda  are  iO.O. 

If  not  then  set  to  0.0  and  turn  constraint  off. 

o)  If  constraint  violations  or  any  lambda  was  <0.0 
then  perform  the  full  Newton-Raphson  procedure 
(call  S461). 

p)  Print  final  results  of  optimization  and  return 
(call  PRIN1 ,  BASIC,  BASICD). 


4.  Input  Arguments: 


See  glossary  for  definition. 


136 


5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

Writes  -  Program  terminated  due  to  singular  matrix 
from  SIMQ  in  OPTFR  STOP 

7. 

Calling  Sequence: 

Call  OPTFR(PRD,PRS,PRR,NTYPE, AMIN, WT, AREA, XC,YC, 

A , KL , P , NE , NDN , NX , NBOU , NSE , NODE , NDT , NDTNX , COND , 
ISTRES , NDC , NW , SIG , ND , IDYN , UP , LOW , IAREA , IMU , 

ALAMBD , W , WARAY ,N5NL, NDCNL , NXNL , NWNL , NDNNL , 
NSENL,IRST,ICHECK,WUYN,PRI,MAXIT,WS,ARAY) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(16118  bytes)  4030  words 

12. 

Subroutine  User: 

LINK2 

13. 

Subroutine  Required: 

FSD , S4  32 , PRIN I , S4  35 , EQSOL , S4  310 , 

ITER , S451 ,ASSEM,S4323, MRPRNT , REGEND , 

SIMQ, S461, BASIC, S4325,BASICD,SHIFTA,S4310B 
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r. 

Subroutine  Name: 

PIVOT 

2. 

Purpose: 

Perform  a  simplex  pivot  about 
element . 

a  g iven  mat  rix 

3. 

Equations  and 
Procedures: 

a)  Interchange  labels  on  row 

and  column  being 

pivoted . 

b)  Update  pivot  element. 

c)  Update  pivot  column  excluding  pivot  element . 

d)  Update  general  element  excluding  those  in 
pivot  row. 

e)  Update  pivot  row  excluding  pivot  element. 

f)  Set  matrix  element  to  zero  if  less  than 
the  minimum  tolerance. 

g)  Pirnt  intermediate  Tableau. 


4. 

Input  Arguments: 

A 

Input  matrix 

IPR 

Pivot  row 

IPC 

Pivot  column 

LABR 

Vector  containing  row  labels 

LABC 

Vector  containing  column  labels 

MM 

First  declared  dimension  of  "A"  in 
calling  program 

M 

Number  of  rows  in  matrix 

N 

Number  of  columns  in  matrix 

TOL 

Minimum  absolute  value  allowed  for  element 
in  matrix 

1PFLAG  Intermediate  printout  code 

IT 

Previous  number  of  iterations 

5. 

Output  Arguments: 

A 

Output  matrix 

LABR 

Vector  containing  updated  run  labels 

LABC 

Vector  containing  updated  column  labels 

IT 

Updated  iteration  number 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call 

PIVOT  (A ,  IPR ,  IPC . LABR , LABC  ,MM  ,M , N  ,T0L ,  IPFLAC. ,  IT) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1110  bytes)  278  words 

12. 

Subroutine  User: 

ITER 

13.  Subroutine  Required:  INTPR 
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1.  Subroutine  Name:  PRIN1 

2.  Purpose:  Prints  structure  design  variables  and  computes 

total  weight. 

3.  Equations  and  a.  Print  out  design  variables. 

Procedures:  b.  Compute  total  structure  weight  and  pr  t. 

4.  Input  Arguments:  A,WT,NE 

See  glossary  for  definitions. 

5.  Output  Arguments:  None 

6.  Error  Returns:  None 

7-  Calling  Sequence:  CALL  PRINl  (A,WT,NE) 

8.  Input  Tapes:  None 

9.  Output  File:  102  -  Unit  6  -  Printer 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (638  bytes)  160  words 

12.  Subroutine  User:  OPTFR, S451 , S461 

13.  Subroutine  Required:  None 
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1.  Subroutine  Name:  REGEND 

2.  Purpose:  Regenerate  (inertial  load  inodes)  and  (redun¬ 

dant  modes)  for  recomputation  of  dynamic  constraints. 

The  mass  matrix  is  computed  first  using  the  Z  matrix 
and  the  A-vector. 

Then  S324R  is  called  to  compute  Pk  in  the  W  array. 
Return  if  NX  =  0. 

If  not,  call  S311D  to  assemble  <J>-^  and  (p. 

Finally  call  S314  to  generate  X^. 

4.  Input  Arguments:  W  -  Contains  DEL(NDN,NW)  room  for  P^(NDN,NW)  and 

Xjt(NX,NW)  and  last  Z  matrix  (NE.NDN)  as  output  from 
BASICD  routine. 

NDN,  NW,  NX,  NE,  A  -  See  Glossary. 

WI  -  Work  storage 


5. 

Output  Arguments: 

Pk  and  Xfc  are  inserted  in  the  respective  positions 
in  the  W  array. 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  RECEND  (W,NDN,NW,NX,NE,A,WI) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1048  bytes)  262  words 

12. 

Subroutine  User: 

OPTFR,S451,S461 

13. 

Subroutine  Required: 

S314  S311D  S324R 

3.  Equations  and 
Procedures: 
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1. 

Subroutine  Name: 

SBTFD 

2. 

Purpose : 

Compute  (B1)1  *  F  *  D1  =  PSI  ~ip. 

3. 

Equations  and 

If  NX=0  return. 

Procedures : 

4. 

Input  Arguments: 

B1  -  Element  matrix  from  structure  cutter  matrix 

element . 

F  -  Flexibility  matrix 

NP  -  No.  of  forces  defined  for  a  particular  element. 
W  -  Work  storage  (5  words  max.) 

D1  -  Element  matrix  from  structure  cutter  matrix. 
NX,NDN-See  Glossary 


5. 

Output  Arguments : 

PSI  (NX.NDN)  Resultant  matrix  (structural  matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  SBTFD  (B1 ,F,NP, NX.NDN, PSI ,W,D1) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1026  bytes)  257  words 

12. 

Subroutine  User: 

S241 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Name: 


SB1IR 


2.  Purpose:  Extract  element  and  reaction  B-matrices  from  struc¬ 

ture  cutter  matrix. 

3.  Equations  and  NS  thru  NF  are  used  to  extract  rows  from  the  A  matrix 

Procedures:  which  are  located  by  the  ICOL  vector. 

fbj^J  is  of  order  (NPj  X  Nx)  where  NPj  is  the  no.  of 
force  components  defined  for  the  element. 

[b^K]  is  of  order  (NBOU.NX)  where  NBOU  is  the  no.  of 
reactions  for  the  problem. 

4.  Input  Arguments:  NDT,NX,NDTNX  -  See  Glossary 

ICOL  -  Row  decoding  vector  for  a  matrix. 

A  -  Structure  cutter  matrix. 

NS  Starting  COL  position  of  element  or 

reactions  in  original  A  before  struc¬ 
ture  cutting. 

NF  Final  COL  position  of  element  or  reac¬ 

tions  in  original  A  before  structure 
cutting. 

NCOL  Length  of  ICOL 

5.  Output  Arguments:  W  contains  the  [b^]  or  [bj^]  matrix  by  rows. 

6.  Error  Returns:  None 

7.  Calling  Sequence:  CALL  SB1IR  (ICOL,A,NDT,NX,NDTNX,NS,NF,NCOL,W) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (780  bytes)  195  words 

12.  Subroutine  User:  S241 

13.  Subroutine  Required:  None 
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1.  Subroutine  Name:  SD1IR 

2.  Purpose:  Extract  element  and  reaction  D-matrices  from 

structure  cutter  matrix. 

3.  Equations  and  NS  thru  NF  are  used  to  extract  rows  from  the  A 

Procedures:  matrix  which  are  located  by  the  ICOL  vector. 

[Dl1]  is  of  order  (NPj  x  NDN)  where  NP^  is  the 
No.  of  force  components  defined  for  the  element. 

[Dj1]  is  of  order  (NBOU  X  NDN)  where  NBOU  is  the 
No.  of  reactions  for  the  problem. 

4.  Input  Arguments:  NDT,NDTNX,NDN,KL,NBOU-see  glossary  for  definitions 

NS,NF-see  procedure  for  definitions 
A-structure  cutter  matrix 

IROW, ICOL-logical  vectors  locating  rows  and  cols. 
NCOL- length  of  ICOL 

5.  Output  Arguments:  W  contains  the  [D^*-]  or  [D-j^]  matrix  by  rows. 

6.  Error  Returns:  None 

7.  Calling  Sequence:  CALL  SD1IR  ( IROW , I COL , A , NDT , NDTNX , NS , NF , NDN , KL , 

NBOU, NCOL, W) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (938  bytes)  235  words 

12.  Subroutine  User:  S241 

13.  Subroutine  Required:  None 
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1.  Subroutine  Name: 


SHIFTA 


2.  Purpose:  Shift  elements  in  a  general  matrix  and 

initialize  elements  in  working  matrix. 

3.  Equations  and 

Procedures:  W1(I)=W(I),  W(I)-0.0 

4.  Input  Arguments :  W  Matrix  containing  elements  to  be  shifted 

to  output  matrix 

M  Number  of  elements  to  shift 

MN  Number  of  elements  in  matrix  W  to 

initialize 


5. 

Out  put  Argument  s : 

W1  Output 

matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  SHIFTA 

(W,W1,M,MN) 

8. 

Input  Files: 

None 

9. 

Output  Files: 

None 

10. 

Scratch  Files: 

None 

11. 

Storage  Required: 

(334  bytes) 

84  words 

12. 

Subroutine  User: 

EQSOL,  OPTFR,  S461 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Name: 


SIMQ 


2.  Purpose: 


Obtain  solution  of  a  set  of  simultaneous  linear 
equations . 


3.  Equations  and  Method:  Elimination  using  largest  pivotal  divisor. 

Procedures:  a)  Forward  solution 

1.  Search  foi  maximum  coefficient  in  column  (pivot). 

2.  Test  for  pivot  less  than  tolerance  (singular 
matrix) . 

3.  Interchange  rows  if  necessary. 

4.  Divide  equation  by  leading  coefficient. 

5.  Eliminate  next  variable, 
b)  Perform  Back  solution. 


4. 

Input  Arguments: 

A  -  Matrix  of  coefficients  stored  columnwise  size 
(N  x  N).  Destroyed  in  computation. 

B  -  Vector  of  original  constants  (length  =  N). 

N  -  Number  of  equations  and  variables. 

5. 

Output  Arguments: 

B  -  Vector  containing  final  solution  values 

(length  =  N) 

6. 

Error  Returns: 

KS;  0  =  normal  solution,  1  =  singular  set  of 

equations . 

7. 

Calling  Sequence: 

CALL  SIMQ  (A,B , N , KS) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1286  bytes)  322  words 

12. 

Subroutine  User: 

OPTFR,  S451 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Name: 


SI  MV 


2.  Purpose:  Invert  a  given  symmetric  positive  definite 

matrix. 

3.  Equations  and 

Procedures:  Solution  is  obtained  using  the  factorization 

by  subroutine  MFSD. 

4.  Input  Arguments:  A  Upper  triangular  part  of  the  given 

symmetric  positive  definite  N  by  N 
coefficient  matrix. 

N  The  number  of  runs  and  columns  in 
matrix  A 

EPS  An  input  constant  which  is  used  as 

relative  tolerance  for  test  on  loss  of 
significance. 

5.  Output  Arguments:  A  Resultant  upper  triangular  matrix 

IER  Error  return  variable 

6.  Error  Returns:  If  input  parameter  N  is  wrong  or  some  radicard 

is  non-positive  the  argument  IER  is  set  to  -1 
and  return  is  made  to  the  calling  routine. 

7.  Calling  Sequence:  Call  SINV  (A, N, EPS, IER) 

8.  Input  Files:  None 

9.  Output  Files:  None 

10.  Scratch  Files:  None 

11.  Storage  Required:  (826  bytes)  207  words 

12.  Subroutine  User:  MINV2 

13.  Subroutine  Required:  MFSD 
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1.  Subroutine  Name:  SIMQ 

2.  Purpose:  Obtain  solution  of  a  set  of  simultaneous  linear 

equations. 

3.  Equations  and  Method:  Elimination  using  largest  pivotal  divisor. 

Procedures:  a)  Forward  solution 

1.  Search  for  maximum  coefficient  in  column  (pivot). 

2.  Test  for  pivot  less  than  tolerance  (singular 
matrix) . 

3-  Interchange  rows  if  necessary. 

4.  Divide  equation  by  leading  coefficient. 

5.  Eliminate  next  variable, 
b)  Perform  Back  solution. 

4.  Input  Arguments:  A  -  Matrix  of  coefficients  stored  columnwise  size 

(N  x  N).  Destroyed  in  computation. 

B  -  Vector  of  original  constants  (length  =  N) . 

N  -  Number  of  equations  and  variables. 

3.  Output  Arguments:  B  -  Vector  containing  final  solution  values  (length  =  N). 

6.  Error  Returns:  K.S;  0  =  normal  solution,  1  =  singular  set  of  equations. 

7.  Calling  Sequence:  CALL  SIMQ  (A,B,N,KS) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (1286  bytes)  322  words 

12.  Subroutine  User:  OPTFR,  S451 


13.  Subroutine  Required:  None 


1. 

Subroutine  Name: 

STDB 

2. 

Purpose : 

Compute 

[B1]T  [F]  [Bi]  =  4> 

3. 

Equations  and 
Procedures : 

If  !IX=0 
Note  B1 

return . 

as  input  is  actually  [BjJ  . 

4. 

Input  Arguments: 

B1  -  Element  matrix  from  structure  cutter  matr 

F  -  Element  flexibility  matrix. 

NP  -  No.  of  force  components  for  element. 
NX  -No.  of  redundant s. 

W  -  Working  Storage. 


5. 

Output  Arguments: 

PHIB 

(NX, NX) 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL 

STDB  (B1,F,NP,NX,PHIB,W) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(962 

bytes)  241  words 

12. 

Subroutine  User: 

S241 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Name:  SUBA 

2.  Purpose:  Generate  [Aq]-^  from  structure  cutter  matrix. 

Extract  lower  triangular  and  upper  triangular  factor 
matrices  and  determine  product. 

This  product  is  Aq. 

Next  get  inverse  of  Aq. 

Store  Ao“^  back  into  structure  cutter  matrix. 

M 

|  M  — |  1  j  « -  >  | 


M 


4.  Input  Arguments:  A  -  Structure  cutter  matrix  after  routine  MFGR. 

M  -  No.  of  rows  in  matrix  A. 

N  -  No.  of  COLS  in  matrix  A. 

WORK,  W0RK1  -  Work  storage 

5.  Output  Arguments:  A 

L  -  Lower  triangle  ractor 
U  -  Upper  trianglf  i'aitor 
BO-  A~'x  load 


6.  Error  Returns:  None 

7.  Calling  Sequence:  CALL  SUBA  (L,U,AO,WORK,WORK1,A.BO,M,N,WA) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (1392  bytes)  448  words 

12.  Subroutine  User:  AA 

13.  Subroutine  Required:  MINV  MPRD 


3.  Equations  and 
Procedures : 
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1.  Subroutine  Name: 


S222 


2.  Purpose: 

3.  Equations  and 

Procedures : 


Compute  basic  structural  matrices  for  dynamic  analysis. 
Read  <p,  ip  and  ft  from  data  set  113. 

Use  A  and  find 


NE  _  NE  _  NE  _ 

■l  t  •  ♦'?  r  •  »*?  r  »i 


Get  inverse  of  (}>  store  on  top  of  $  (M1NV). 


Form  mass  matrix  from  the  Z  matrix. 


4. 

Input  Arguments: 

A,  NE,  NX,  NDN,  -  See  Glossary  for  Definitions 
W  -  work  storage. 

5. 

Output  Arguments: 

PHI,  PSI,  OMEGA  -  See  Glossary  for  Definitions 
M  -  Mass  matrix 

Z  -  Z-matrix 

6. 

Error  Returns: 

7. 

Calling  Sequence: 

CALL  S222  (PHI,PSI,OMEGA,W,A,NE,NX,NDN,M,Z) 

8. 

Input  Tapes: 

113 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(2748  bytes)  687  words 

12. 

Subroutine  User: 

BAS I CD 

13. 

Subroutine  Required: 

MINV  MRPRNT 
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1.  Subroutine  Name: 

S231 

2.  Purpose: 

Compute  element 

forces  or 

reactions  as  function  of 

redundant  X  and 

load  P. 

3.  '  Equations  and 

(Si)  =  [b^HX] 

+  [DiHP] 

(2.3.1) 

Procedures : 

{  R  }  [b ! R]  [X] 

+  [Dr](P] 

(2.3.1) 

4.  Input  Arguments: 

NTYPE,NE,NDT,NDTNX,X,NX,NDN,P,NBOU,KL,NCOL,NLL,NSE 

See  Glossar''  for  Definitions. 


NW  -  Length  of  S  output. 


5. 

Output  Arguments: 

S  -  Either  forces  or  reactions  (NSR >  0  or  NSR <  0 
respectively) . 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  S231  (NTYPE.Il  2  ,NE,NDT  ,NDTNX,X,NX,NDN, 

P , NSR, NBOU , KL , W, NW , NCOL , S , NLL , NSE) 

8. 

Input  Tapes: 

112,  115 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(2012  bytes)  503  words 

12. 

Subroutine  User: 

BASICD,  BASIC 

13. 

Subroutine  Required: 

S231A 
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1.  Subroutine  Name: 

2.  Purpose: 


3.  Equations  and 
Procedures : 


4.  Input  Arguments 


5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Input  Tapes: 

9.  Output  Tapes: 

10.  Scratch  Tapes: 

11.  Storage  Required: 

12.  Subroutine  User: 

13.  Subroutine  Required: 


S231A 

Computes  start  and  final  columns  of  element 
B  matrix  in  a  (structure  cutter  matrix). 


Take  into  account  the  following  facts  about 
the  structure  cutter  matrix. 

1)  The  reactions  are  i:i  the  first  NBOU  columns. 

2)  Axial  and  shear  web  elements  have  only  1 
force  column  in  B~element  matrix. 

3)  Triangle  has  3  columns  and  quadrilateral 
element  has  5  columns  in  their  B  element 
matrices . 

NTYPE  Element  type  (1-4) 

NE  No.  of  elements 

NSR  Element  No. 

NBOU  No.  of  reactions 

NS.NF  starting  and  final  column  no's,  respectively . 

None 

Call  S231A  (NTYPE, NE, NSR, NS.NF, NBOU) 

None 

None 

None 

(574  bytes)  144  words 

S241.S231 

None 
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1. 

Subroutine  Name: 

S24 1 

2. 

Purpose : 

Compute  element  normalized  matrices 

and  put  on 

tape. 

Write  b^  and  D^ 

matrices  to 

tape . 

3. 

Equations  and 

Cii]  - 

[b/i1  tv  [b; 

Procedures : 

]  (NX  x  NX) 

(2.4.1) 

[$i]  = 

[b;  iT  [^3  [D. 3 

(NX  X  Hon) 

(2.4.2) 

[fy]  = 

CD^1  [?,]  [0^ 

(ndn  x  ndn 

)  (2.4.3) 

BMG21  is  used  to  compute  the  Z  matrix  a  row  of 
which  is  put  on  each  bar  matrix  record  of  unit 
113. 


4.  Input  Arguments:  NTYPE , KL ,  NE , NBOU ,  NDT , NDTNX , NDN , NX , NCOL-see  glossary 

IROW  Row  label  vector  (structure  cutter) 

ICOL  Column  label  vector  (structure  cutter) 

A  A  matrix  (structure  cutter) 

NTAPE  Output  file  unit  number 

W  Work  array 

BMG21  Element  mass  contribution  vector 


5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S241, (NTYPE, IROW, ICOL, A, KL,NE, NBOU, NDT, 
NDTNX , NDN , NX , NTAPE , NCOL , W , BMG21 ) 

8. 

Input  Tapes: 

NSS1 

9. 

Output  Tapes: 

113,  115,112 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(3228  bytes)  807  words 

12. 

Subroutine  User: 

TAPE1I 

13. 

Subroutine  Required: 

STDB , SBTFD ,SB1IR ,SD1IR, 

S 2 3 LA , MRPRNT , TRPRNT 
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1.  Subroutine  Name: 


S311 


2.  Purpose: 

Assembles 

bar 

matrices 

into  4>,^  and  SI 

3.  Equations  and 
Procedures : 

[0>] 

NE 
=  Z 
i  =1 

1 

AT 

c^] 

(3.1.1) 

0] 

NE 

=  Z 

i=l 

1 

Ai 

[^] 

(3.1.2) 

NE 
=  Z 
i=l 

1 

AT 

[^] 

(3.1.3) 

4. 

Input  Arguments: 

A,NE,NX,NDN  -  see  glossary 

W-  work  storage 

5. 

Output  Arguments: 

PHI , PSI .OMEGA  -  see  glossary 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S311(PHI, PSI, OMEGA, W.A.NE, NX, NDN) 

8. 

Input  Tapes: 

113 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(2146  bytes)  537  words 

12. 

Subroutine  User: 

BASIC 

13. 

Subroutine  Required: 

MINV2  .MRPRNT 
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1.  Subroutir 2  Name: 

2.  Purpose: 


S311D 

Compute  PHI-1  and  PSI  using  (3.1.1)  and  (3.1.2) 
for  dynamics  analysis. 


3.  Equations  and 
Procedures : 


M  -  I  ^  [$,]  (NX  *  NX) 

i=l 

HE  , 

[<J>]  =  J-  at  (NX  x  Ndh) 


4. 

Input  Arguments: 

A,NE,NX,NDN  -  see  glossary 

V-vork  storage 

5. 

Output  Arguments: 

PHI.PSI  -  see  glossary 

6 . 

Error  Returns: 

None 

7. 

Calling  Sequences: 

Call  S311D  (PHI , PSI , W, A,NE ,NX , NDN) 

8. 

Input  Tapes: 

113 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1298  bytes)  325  words 

12. 

Subroutine  User.' 

REGEND 

13. 

Subroutine  Required: 

MINV 
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1. 

Subroutine  Name: 

S314 

2. 

Purpose: 

Computes  redundants 

3. 

Equations  and 

fx>  =  -  or1  o]  fp}  (3.1.4) 

Procedures: 

4. 

Input  Arguments: 

P ,NX,NDN ,NLL , PHI , PSI ,NLL-see  glossary 
W-work  storage 

5. 

Output  Arguments: 

X  Vector  containing  computed  redundants 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S314  (PHI,PS1,P,X, NX.NDN , W,NLL) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1164  bytes)  291  words 

12. 

Subroutine  User: 

REC END, BASICD, BASIC 

13. 

Subroutine  Required: 

MRPRNT 
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1. 

Subroutine  Name: 

S315 

2. 

Purpose : 

Compute  reduced  displacement  vector  -  static 
analysis . 

3. 

Equations  and 
Procedures : 

{&}  =  i>]T  (X)  +  [fl]  {p}  (3.1.5) 

4. 

Input  Arguments: 

NLL, PSI, OMEGA, X,P,NDN, NX, NLL-see  glossary 

5. 

Output  Arguments: 

DISP  Vector  containing  displacements 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S315  (PSI ,X ,0MEG , A, P, DISP , NDN , NX , NLL) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(882  bytes)  221  words 

12. 

Subroutine  User: 

BASIC 

13. 

Subroutine  Required: 

None 
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1. 

Subroutine  Name: 

S316 

2. 

Purpose: 

Calculate  flexibility  matrix  for  dynamics 
analysis . 

3. 

Equat ions  and 
Procedures : 

[f]  •  [n]  -  M1  [o'1]  M  (3.1.6) 

(NDN  x  NON) 

4. 

Input  Arguments: 

W  Work  storage 

OMEGA, PS I, PHI, NX, NDN  -  see  glossary 

5. 

Output  Arguments: 

OMEGA  Replaced  by  (f) 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S316  (OMEGA, PSI, PHI, NX, NDN, W) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(1032  bytes)  256  words 

12. 

Subroutine  User: 

BASICD 

13. 

Subroutine  Required: 

MRPRNT 
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1. 

Subroutine  Name: 

S321 

2. 

Purpose : 

Calculate  Eigenvalue  matrix  for  a  ■ 
analysis . 

3. 

Equations  and 
Procedures : 

Q.  .  -  F.  .  (M.M.  )**'  (3.2.1) 

U  iJ  i  3 

4. 

Input  Arguments: 

F  Flexibility  matrix 

M  Mass  matrix  (only  diagonal 

NDN  Reduced  DOF 

5. 

Output  Arguments: 

F  Replaced  by  Q 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S321  (F,M,NDN) 

8. 

Input  Tapes : 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(680  bytes)  170  words 

12. 

Subroutine  User: 

BAS  I CD 

13. 

Subroutine  Required: 

MRPRNT 
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1.  Subroutine  Name: 


S323 


2-  Purpose: 

3.  Equations  and 

Procedures: 

4 .  Input  Arguments : 

5.  Output  Arguments: 

6.  Error  Return: 

7.  Calling  Sequence: 

8.  Input  Tapes: 

9.  Output  Tapes: 

10.  Scratch  Tapes: 

11.  Storage  Required: 

12.  Subroutine  User: 

13.  Subroutine  Required: 


Calculate  inertia  load  mode  for  dynamic  analy: 

[pki]  =  fV*  V  (3-2J) 

W  =  1/X  (3.2.2) 

M  Mass  matrix 

Q  q  (Eigen  vector) 

NDN  Red  DOF 

NW  No  of  modes 

Lamda  X 

Q  -  replaced  by  Pv 
None 

Call  S323  (M,Q, NDN, NW, LAMDA) 

None 

None 

None 

(540  bytes)  135  words 

BASICD 

None 
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1. 

Subroutine  Name: 

S324 

2. 

Purpose: 

Calculate  displacement  modes  in  dynamic  ana 

3. 

Equations  and 

i  i  /  2 

Procedures: 

4k  -pk/“k  Mi  l5-2-4) 

4. 

Input  Arguments: 

P  Inertial  load  mode 

M  Mass  matrix  (diagonal) 

NDN  Red  DOF 

NW  No  of  modes 

W  Work  storage 

5. 

Output  Arguments: 

DEL  -  Displacement  modes 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S324  ( F , M , NDN , DEL , NW , V ) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(528  bytes)  132  words 

12. 

Subroutine  User: 

BASICD 

13. 

Subroutine  Required: 

None 
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1. 

Subrout ine  Name : 

S324R 

2. 

Purpose: 

Solves  for  Pj^  from  for  dynamics  reanalysis 

(optimization). 

3. 

Equations  and 

?k  =  ak  *  Mi 

Procedures : 

4. 

Input  Arguments: 

DEL  Displacement  modes 

M  Mass  matrix 

NW  No.  of  modes 

NDN  No.  of  reduced  DOF 

5. 

Output  Arguments: 

PK  Inertial  lode  mode. 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  5324R  (DEL.M.PK.NW.NDN) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(470  bytes)  118  words 

12. 

Subroutine  User: 

REGEND 

13. 

Subroutine  Required: 

None 

I 
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1. 

Subroutine  Name: 

S4310 

2. 

Purpose: 

Computes 

1st  partial 

of 

go  with  respect  to  X 

3. 

Equations  and 
Procedures: 

Each  row 

of  the  NE  x 

matrix  is 

4.  Input  Arguments: 

5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence: 

8.  Input  Tapes: 

9.  Output  Tapes: 

10.  Scratch  Tapes: 

11.  Storage  Required: 

12.  Subroutine  User: 

13.  Subroutine  Required: 


A  .a* 
i  i 


3Yi 

3Si 


(4.3.10) 


A, SIG , Y, S , NTYPE , AREA.XC, YC, NE,NSE,NDT ,NX , IR, IC, 
NTRAN  S , NR , NC , NCOL , NDTNX , NROW , NBOU  -  see  glossary 

W  work  storage 

Each  element  computed  is  inserted  in  ARAY  by 
S4310A. 


None 

Call  S4310  (A,SIG,Y,S,NTYPE,AREA,XC,YC,NE, 
NSE , NOT , NX, IR, IC , NTRANS , ARAY , NR , NC , NCOL , 
NDTNX, NROW, NBOU, W) 

112 


None 


None 

(2600  bytes)  650  words 
OPTFR , S4 51 , S4 61 , EQSOL 
SQRT, S4  310A, TRPRNT 
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1. 

Subrout ine  Name : 

S4310A 

2. 

Purpose : 

Output  of  calculated  matrix  element  to  a 
temporary  file. 

3. 

Equations  and 
Procedures: 

Write  row  number,  column  number,  transpose 
control,  element  value  and  summing  control 
to  a  temporary  data  file. 

4. 

Input  Arguments: 

1B1  Row  number 

IB2  Column  number 

NT  Transpose/symmetry  control 

ANS  Element  data 

NS  Summary  control 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S4310A(A,NR,NC,IBl, IB2, NT, ANS , NS) 

8. 

Input  Files: 

None 

9. 

Output  Files: 

114 

10. 

Scratch  Files: 

None 

11. 

Storage  Required: 

(378  bytes)  95  words 

12. 

Subroutine  User: 

S4614 , S4621 , S4622 , S4  32 , S4  35, S432  3, S4  310, 

S4  325 , S466 

13. 

Subroutine  Required: 

None 
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1. 

Subrout ine  Name : 

S4 31 OB 

2. 

Purpose : 

Insert  calculated  element 

into  analysis  matrix 

3. 

Equations  and 
Procedures : 

Reads  element  (ANS) ,  row 

and  column  positions 

(IB1.IB2),  transpose/lower  triangle  control 
(NT),  and  summing  control  (NS)  from  input 
file.  Stores  elements  in  output  matrix  "A" 
as  specified  by  control  variables 
NS=0  =1 

NT 

+1  A( IB1 , IB2) =ANS  A(IB1,IB2)=A(IB1, 132)+ANS 

+2  A( IB2 , 1B1 ) =ANS  A ( IB2 , 1B1) =A ( IB2 , 1B1 ) +ANS 

-1  A ( IB1 , IB2 ) =ANS  A(1B1, I B2 ) =A ( 1B1 , 1B2) +ANS 

(for  all  IB1<IB2) 

-2  A(IB2, IB1)=ANS  A( 1B2 , 1B1 ) =A( IB2 , 1B1)+ANS 

(for  all  IB1<IB2) 


4. 

Input  Arguments: 

A 

Matrix  to  be  assembl 

M 

Number  of  rows 

N 

Number  of  columns 

5. 

Output  Arguments: 

A 

Assembled  matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call 

S4310B  (A,M,N) 

8. 

Input  Files: 

114 

9. 

Output  Files: 

None 

10. 

Scratch  Files: 

None 

11. 

Storage  Required: 

(798 

bytes)  200  words 

12. 

Subroutine  User: 

EQSOL , OPTFR,S451,S461 

13. 

Subroutine  Required: 

None 

1  bla 


1. 

Subroutine  Name: 

S432 

2. 

Purpose : 

Generates  1st  partial  of  gA  with  respect  to  A. 
Also  GA  is  computed  if  requested. 

3. 

Equations  and 
Procedures : 

*  IF 

gAx  =  Xi  “  1  IGA  (4-3-1) 

A  A1 

Ai* 

^  ?  6 i j  (diagonal  matrix)  (4.3.2) 

j  j 

4. 

Input  Arguments: 

AMIN,A, NE, IR, IC ,NTRANS , NR , NC, IGA  -  see  glossary 

5. 

Output  Arguments: 

ARAY  matrix  for  storage  of  computed  partial 

VECT  Vector  for  storage  of  gA 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S432  (AMIN,A,NE,IR,IC, NTRANS , ARAY, NR, NC, 
VECT, IGA) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(794  bytes)  199  words 

12. 

Subroutine  User: 

OPTFR, EQSOL 

13. 

Subroutine  Required: 

S4310A 
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1.  Subroutine  Name: 


S4323 


2.  Purpose: 

3.  Equations  and 

Procedures : 


4.  Input  Arguments: 


5.  Output  Arguments: 

6.  Error  Returns: 

7.  CaLLing  Sequence: 

8.  Input  Tapes: 


Generates  4  oartials  3g  3gA 

3A  '  3x  ’ 


,  4^  and  gA  and  gx 

a  A  o  X 


gA  =  A*  ({tPi}T  (x)  +  [n1] 
gx  =  [$]  (x)  +  [^J  (p) 


(4.3.21) 

(4.3.35) 


SjL_l  =  _  1  -  j )  .  X  +  n(j)miPi](4.3.23) 

3A .  VT  A.  2  1  V  7  1 


14  -  **« 


3xk  A; 


(4.3.23) 


*£*  =  -  ±o  ([$  ]  (x)  +  [$.]  { P) ) 5th  eol. 

3A  AJ  3  J  (4.3.38) 


iLSi  =  a 

3xk  T ik 


(4.3.40) 


A,NE,P,NDN,NDC,NX,NR,NC,ND,X,NTRANS-see  glossary 

IR1  Vector  containing  displacement  con¬ 

straint  row  labels 

XR2  Vector  containing  redundant  constraint 

row  labels 

Id  Vector  containing  unknown  design 

variable  column  labels 

IC2  Vector  containing  redundant  variable 

column  labels 

DSTARU  Vector  containing  upper  displacement 

limits 

DSTARL  Vector  containing  lower  displacement 

limits 

W  Work  storage 

IGX  Control  to  compute  gx 

IGD  Control  to  compute  gA 

ARAY  is  modified  by  the  partials 

VECT  Storage  for  gA  and  gx 

None 

Call  S4323  (A,NE,P,NDN,NDC,NX,IR1  ,IR2, 1C1,  IC2, 

ARAY , NR , NC , DSTARU , DSTARb , ND , X , W , NTRANS , VECT , 

ICX, IGD 

113 
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9.  Output  Tapes: 


None 


10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(3483  bytes)  871 

words 

12. 

Subroutine  User: 

OPTFR, S451 , S461 , 

EQSQL 

13. 

Subroutine  Required: 

S4310A 
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1.  Subroutine  Name: 

2.  Purpose: 


S4325 

Generate  Frequency  Constraint  Optimizing 
Factors, 


2  2^2 
3&v  .  Sw  »  ISiC—  summed  into  V7  TT 

9A  »A,  ^  3V  K 


3.  Equations  and 
Procedures : 


4  •  ^  ^<WV.  +  r  ^  <V„<Vp  +  4  v  (V,  <Vr 


(P„v 


(4.3.31) 


(v*)2  Ag  ffig  Bgs  K  s 


|aii  .  *  2  cxk)t  t  ]  fpK>  +  ( vT  [  nj  5V] 


3  3 


non 


$  A  „_,A 


+  .i-  7  ,'P  }  rz.]fP„}  Where  Z=  [Z]  .Jsmj 

+  (W  )2  ’•  >-  jJ-  K;  NE  M/ 

1  s 


(4.3.32) 


^2  1  ? 

3.  Stf  =  — ~~~t  r 

9A  3Ak  A^ j  [ 


>[  ^/Vm  +  24,I  (Xi)n(Pi)p  +  (pi)q(Pi)r]6jK 


,2  (p.) 

(«i*>2  (BgrfBg  Ag)3  1 


4.  Input  Arguments: 


A,WS,NE,NDN,NX,NR,NC  -  See  Glossary  for 
definitions 

PK  -  Vector  containing  inertial  load  mode 


XK  - 


"  redundant  mode 


NWC  -  Number  of  modes 
W  -  Work  storage 

Z  -Matrix  containing  mass  contrioutions  ot 

elements /DOF  , 

I  MU  -  Vector  Identifying  active  freq.  constraint 
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MUW  -  Vector  containing  frequency  p's 
IRI  -  Assembly  vector  containing  row 
labels  (3gw2/3A3A) 

IC2  -  Assembly  vector  containing  column 
labels  Ogw2/3A3A)  ,3gw/3A 
IGP  -  Assembly  vector  containing  row 
labels  (dw/9A) 

NGW  -  Control  switch  to  compute  GW 
NDW  -  "  "  "  "  3GW2/3A3A 

NT  -  Transpose  control 

5.  Output  Arguments:  ARAY  -  External  matrix  used  for  linear  pro¬ 

gramming  and  Newton  Raphson  solutions  GW  - 
vector  to  store  computed  gw  values 

6.  Error  Returns:  None 

7.  Calling  Sequence:  Call  S4325  (A^S.PK.XK.NWC.NE.NDN'.NX.W.GW, 

ARAY,NR,NC,Z , IMU,MUW, 1R1 , IC2 , IGP , NGW, NDW, NT) 

8.  Input  Tapes:  113 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (3228  BYTES)  807  words 

12.  Subroutine  User:  OPTFR,S451,S46l,EQSOL 

13.  Subroutine  Required:  S4310A,S4325A .MRPRNT 
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1.  Subroutine  Name:  S4325A 

2.  Purpose:  Compute  f\  for  gw  and  both  partials  of  gw 

3.  Equations  and 

Procedures : 

T  T 

f.  =  (Xv)  [PHl]  {X„}  +  2.0  (x  )  [PSI]  (p  1 

IK  K  K.  K 

+  (PV)T  [OMECA]  (p  >  (4.3.26) 

K.  K. 

4.  Input  Arguments:  PHI , PSI , OMEGA, NX, NON  -  See  Glossary  for 

def inition 

X  -  Vector  containing  redundant  mode 
P  -  Vector  containing  inertial  load  mode 
W  -  work  storage 

5.  Output  Arguments:  ANS  -  f  is  returned  to  calling  program 

in  this  variable 

6.  Error  Returns: 

7.  Calling  Sequence:  Call  S4325A(PHI, PSI, OMEGA, X,P, ANS, NX, N  N,W) 

8.  Input  Tapes:  None 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  1096  bytes  (274  words) 

12.  Subroutine  User:  S4325 

13.  Subroutine  Required:  None 
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1 


Subroutine  Name: 


S433 


2.  Purpose:  Compute  stresses  and  V.  for  an  element 

and  print  stresses  if  requested. 

3.  Equations  and  Stresses  are  computed  as  follows. 

Procedures: 


Bar  element  Sx  =  S^  Sy  =  0  TXy  =  0 

T~ 


Shear  Panel  Sx  =  0  S  =  0  T  -  —  Lr  ■  - 

y  xy  a/area 


Triangle 

sx  =  -si  S  =  -1 - 

aJ  AREA  ’  y  A  J  AREA 

T  =  Jj _ 

'  " '  A  /AREA 

Quad 

S  =  -ll_  +  S2  yc  s 

S3 

S4  X(: 

-  — 

A  AREA  A*AREA  ’  V 

A J  AREA 

A*AREA 

168 


10. 

Sciatch  Tapes: 

None 

11. 

Storage  Required: 

(1344  BYTES)  336  words 

12. 

Subroutine  User: 

BASIC,  BAST CD 

13. 

Subroutine  Required: 

None 
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1. 

Subroutine  Name: 

S435 

2. 

Purpose : 

Generates  partial  of  go  with  respect 

to 

A  and  go 

3. 

Equations  and 
Procedures : 

1  .  Y1 

i  ;  aso1  =  -  Yi  . 

AjO 

3A-  a  2a  .  6ii 

3  Ai  Oj* 

a)  For  each  unknown  design  and  active  stress  constraint  compute 

3go^/9Aj.  Assemble  in  coefficient  matrix  using  routine  S4310A. 

b)  For  each  active 
constant  vector 

stress  constraint  compute  g,.  and  assemble  in 

4. 

Input  Arguments: 

A,SIG,Y,NE,IR,IC,NTRANS,NR,NC, IGS 

See  Glossary  for  definitions. 

5. 

Output  Arguments: 

ARAY  -  Assembled  matrix  containing 
VECT  -  Vector  containing  gs  values 

9Gs 

/3A 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

Call  S435(A, SIG,Y,NE,IR,IC,NTRANS, ARAY, NR, NC, VECT 
IGS) 

8. 

Inpur  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(842  BYTES)  211  words 

12. 

Subroutine  User: 

OPTFR,  S451,  S461 ,  EQS<{>L 

13. 

Subroutine  Required 

:  S4310A 
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Subroutine  Name: 


Purpose : 

Equations  and 
Procedures : 


Input  Arguments: 


Output  Arguments: 
Error  Returns: 

Calling  Sequence: 


Input  Tapes: 


S4S1 

Perform  a  partial  Newton  Raphson 

a.  Compute  matrix  size 

b.  Print  design  variables  and  perform  a  basic 
analysis  and  update  dynamic  variables.  Call 
PR.INI ,  BASIC ,  REGEND 

c.  Iteration  loop 

1)  Initilize  working  arrays 

2)  Compute  and  assemble  the  partials  G*S/A+, 
G*S/X,  G*D/A+,  G*D/X,  GX/A+ ,  GX/X,  G*W/A+ 
into  array  WARAY  and  the  constraints  G*S, 
G*D,  GX  and  G*W  into  VECTOR  "G"  (Call 
ASSEM,  S435,  S4310,  S4323,  S4325) 

3)  Print  intermediate  data  if  requested 

4)  Solve  for  Delta  A  and  Delta  X.  Call  SIMQ 

5)  Check  for  convergence  of  all  A  and  X  and 
update  using  New  =  Old  +  Delta 

6)  Print  updated  design  (call  PRINI) 

7)  Perform  a  static  analysis  and  update  dynami 
variables  (Call  BASIC , REGEND) 

8)  If  any  variable  did  not  converge  and  design 
variables  did  not  fall  below  lower  limits 
and  less  than  MAXIT  iterations  have  been 
performed  then  repeat  loop. 

d.  Check  all  constraints  for  violations.  If 
violations  occur  set  error  switch,  raise 
design  variables,  print  new  design  and 
perform  a  static  analysis  and  update  dynamic 
variables  (call  PRINI, BASIC, REGEND) 

PRD, PRS , PRR, NASTAR, NSSTAR, NDSTAR,NXSTAR,NWSTAR,W, 
G, WARAY , IYI , A,SIG, NTYPE, AREA,XC, YC,NE, NSE,NDT,NX, 
NCOL, NOTNX, NBOU, ISI, NDC, IMUSL,IMUDL, I AREA, IMU, 

P , NDN , UP , LOW , I X , I DYN , NW , IMUXL , I MUWL , I CC , ALAMBD , 

AMI N ,WT , KL , NODE , ND , COND , IDEL . NSNL , NDC NL , NXNL , 

NWNL , I CHECK , WDYN , PRI , WS , MAX T ' ' , I PK , I XK , I Z 
See  Glossary  for  definition. 

W, A, WDYN,  See  Glossary  for  definition 

Write:  MATRIX  is  singular  IN  S451  -  Partial  N-R. . 
STOP  IERR:  0  =  converged,  1  =  did  not  converge, 

2  =  constraint  violation 

Call  S451 (PRD, PRS, PRR, NASTAR, NSSTAR, NDSTAR,NXSTAR, 
NWSTAR,W,G, WARAY, IYI , A, SIG, NTYPE , AREA, XC , YC, NE, 

NSE , NDT , NX , NCOL , NDTNX , NBOU , I S I , NDC , I MUSL , I MUDL , 

I AREA, IMU, P, NDN, UP, LOW, IX, IDYN, NW, IMDXL, IMUWL, 

IERR , ICC , ALAMBD , AMI N , WT , KL , NODE , ND , COND , IDEL , NSNL , 
NDCNL , NXNL , NWNL , 1 CHECK , WDYN , PRI , WS , MAXIT , I PK , I XK , 
IZ.ARAY) 

None 


9.  Output  Tapes: 


None 


10.  Scratch  Tapes:  None 

11.  Storage  Required:  (7246  BYTES)  1812  words 

12.  Subroutine  User:  OPTFR 

13.  Subroutine  Required:  SIMQ, S435 .ASSEH  ,  BASIC, PRINI ,S4310 ,S4 323 ,  S4325, 

TESTR, MRPRNT , REGEND , S4  3 1  OB 
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1.  Subroutine  Name: 

2.  Purpose: 

3.  Equations  and 

Procedures : 


I 


|  i 

I 

> 

k 

■ 

4.  Input  Arguments: 


5.  Output  Arguments: 


S461 

Performs  a  full  Newton-Raphson . 

a)  Iteration  ' oop 

1)  Compute  matrix  dimensions  and  initialize 

working  storage  -a 

2)  Compute  and  ass  a  the  >artials  G*S/A+, 

C*S/X,C*D/A+,G*.  'X/A-*  ,GX/X,G*W/A+,C2*W/A-t-A+, 

and  the  constraints  *1),GX,G*W  and  store 

in  matrix  WARAY . 

(CALL  ASSEM,S435,S4j10,oA323,S4325>) 

3)  Compute  and  assemble  the  partiais  G2L/A+A+, 
G2L/A+X.G2L/XX  and  store  in  matrix  WARAY. 

(CALL  ASSEM, S46 14 , S4621 , S4622 ) 

4)  Compute  and  assemble  the  partiais  CL/A+  and 
GL/X  and  store  in  matrix  WARAY. 

5)  Print  lower  triangle  matrix  WARAY  if  requested. 

6)  Solve  for  DELTA  A+,  DELTA  X  and  DELTA  U*  using 
routine  GELS. 

7)  Check  for  convergence  of  all  design  variables 
and  redundant s  and  update  to  new  values.  Set 
error  switch  if  any  did  not  converge. 

8)  Check  for  convergence  of  all  U*,  update  to  new 
values,  verify  all  U*  are  positive,  if  not  posi 
tive  set  to  0.0  and  turn  constraint  OFF. 

9)  Verify  that  all  design  variables  are  greater 
than  0.0.  If  not  set  to  minimum  and  turn  con¬ 
straint  "ON". 

10)  Print  new  design,  perform  a  static  analysis 

and  update  dynamic  variables.  CALL  PR1N1, BASIC 
REG END. 

.'!)  If  any  variable  did  not  converge  and  the  number 
of  iterations  is  less  than  maximum,  repeat  iter¬ 
ation  loop. 

b)  If  design  converged  check  all  constraints: 

1)  Check  all  displacement  constraints  for  viola¬ 
tions  and  turn  "ON"  if  violated. 

2)  Check  all  stress  constraints  for  violations 
and  turn  "ON"  if  violated. 

3)  If  any  violations  occurred,  repeat  Step  A. 

c)  Compute  the  LaGrange  multipliers  for  all  design 

constraints  which  are  "ON"  using  routines  ASSEM, 

S435,S4323,S4325,  and  S432. 

d)  Compute  total  weight  from  LaGrange  multipliers. 

PRD , PRS , PRR , NASTAR , NSSTAR , NDSTAR , NXSTAR , NWSTAR , NE , NX , 

W , WARAY , I AREA , IMU , IMUAI. , IMG  SL , IMUDL , IMUXL , IMUWL , A , S I G , 

IX, IYI , ISI , IDEL ,NTYPE, AREA, XC ,YC ,NSE,NDT,NCOL ,NDTNX , 

NBOU , NDC , P , NDN ,  UP , LOW, ND , NW, ALAMBD, WT , COND, AMI N , KL, 

NODE , ICC , NSNL , NDCN1. , NXNL , NWNL , WDYN , PRT , WS , MAXIT . IPK, 

IXK, IZ, IDYN 

See  Glossary  for  Description. 

W, A, WDYN  -  See  Glossarv  for  Definitions. 
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6.  Error  Returns: 


Write:  Matrix  is  singular  in  S461  -  full  N.r-... 

I ERR  STOP 

7.  Calling  Sequence:  CALL  S461  (PRD,  PRS,  PRR,NASTAR,NSSTAR,NDSTAR,NXSTAR, 

NWSTAR , NE , NX , W , WARAY , 1AREA , I MU , IMUAL , IMUSL , IMUDL , 

IMUXL , IMUWL ,A,SIG,IX,1Y1,IS1, IDEL ,NTYPE, AREA,XC , YC , 

NS E ,  NDT ,  NCOL ,  NDTNX ,  NBOil,  NDC,  P,  NDN ,  UP ,  LOW ,ND,NW,  ALAMB I) , 
WT, IERR.COND,  AMIN, KL, NODE, ICC, NSNL , NDCNL , NXNL , NWNL , 
ICHECK, WDYN , PRI , WS , MAX1T , 1PK , IXK, IZ , 1 DYN, ARAY) 


8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11  . 

Storage  Required: 

(14040  bytes)  3510  words 

12. 

Subroutine  User: 

OPTFR 

1  1. 

Subroutine  Required: 

CELS, S4 35, ASSEK, BASIC, PRIN1 , S4310, S4 323, S4 325, 

S4614 , S4621 , S4622 , TESTR.MRPRNT, REGEND, SHIFTA, S4310B 
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I.  Subroutine  Name:  S4614 


2.  Purpose: 

3.  Equations  and 

Procedures : 


4.  Input  Arguments: 


Second  partials  with  respect  to  A  A  . 

•  i  K 


2  Y. 


9  A  ,  3A^ 


A  . 


3  *  *d  1R 

A-i  0  j 


(4.6.4) 


42g  1 


3  * 

\)  Ai 


XY  +  i2(J)mi  P"5  UK 


(4.6.8) 


r-Tg1 

SX  2  -i  1 

~^t  '  ,  ~  _  l<$'.  „Xv  +  .  ;  Pm! 

>A .  ’A  3  l?  *m  1 K 

'  k  A I 


(4.6.12) 


i  * 


+•  f 

JAJ  'ak 


•  r,.. 

-t-  + 

,AJ  ,ak 


'  K 


+  + 
,AJ 


■»x 

+  4- 

,Ai  ak 


2  i  * 

•  *  ’  8W 

+  ,,i  - — - 

W  +  -r 

,aj  'ak 


(s .6. 20) 


NOTE:  Tlie  4lh  term  in  Eq .  4.6.20  is  summed  into  the 
coefficient  matrix  bv  routine  S4325. 


a)  Each  term  ^  Eq .  4.6.20  is  computed  and  summed 
into  the  lower  triangular  coefficient  matrix 
by  routine  S4310A. 


NR , NC , A ,  S IG , X , P , UP , LOW , ND , NE , NDC , NX , NDN , W 


MUS  =  Vector  containing  uc 

MUD  =  Vector  containing  lia, 

MUX  =  Vector  containing  i:x 

YI  =  Element  force  vector 

IRA  =  Output  ROW  and  COL  vet  ter 

TRS  =  Stress  constraint  active  vector 

IRD  =  Displ.  constraint  ictive  vector 

IRX  =  Redundant  constraint  active  vector 

W  ~  Working  storage  vector 


5 . 

Output  Arguments: 

ARAY  -  Coefficient  matrix 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALI.  S4614  (ARAY ,  NR ,  NC , MUS , MU D, MUX ,  Y I ,  A .  S 1 0 , } 

i.P.UP, 

■OW , NL>, NE , NDC , NX , NDN ,  IRA,  IRS,  IRD.  1  RX . W) 
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I 


iD-A  1 28  235 


FORCE  METHOD  OPTIMIZATION  II  VOLUME  II 
(U)  BELL  AEROSPACE  TEXTRON  BUFFALO  NY 
NOV  82  AFWAL-TR-82-3088-VOL-2  F33615-8 


USER' S  MANUAL 
J  R  BATT  ET  AL. 
-C-3214 

F/G  9/2 


UNCLASSIFIED 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  OF  STANDARDS  I96J-A 


8. 

Input  Tapes: 

Unit  13  containing  <J>  ^  &  H. 

9. 

Output  Tapes: 

None 

10: 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(2338  bytes)  577  words 

12. 

Subroutine  User: 

S461 

13. 

Subroutine  Required: 

S4310A 
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1.  Subroutine  Name: 


S4621 


2 .  Purpose : 

3.  Equations  and 

Procedures: 


4.  Input  Arguments: 


5.  Output  Arguments: 

6.  Error  Returns: 

7.  Calling  Sequence : 

8.  Input  Tapes: 

9.  Output  Tapes: 


Compute  3L 


i* 


3AJ3XK 


.2  i* 

3  8o 

3AX 


-2  1*  .2  i 

i*3  gA  .  i  3  *x 
+  +  V 

1  saT^  x 


j 


j 


Each  partial  in  the  expression  is  computed  and  summed 
into  the  lower  symmetric  analysis  matrix  ARAY 


.2  i 
3  8o 

3aj3xk 


3Yi  bA„  6ij  (4.6.6) 


A.  V*  3S  <*>  1  *K 

l  i  i 


J2  1 
3  8A 

3aj3xk 


2  *  ^(j)Ki 

AJ  Ai 


(4.6.9) 


.2  i 

3  8X 

3Aj3XK 


~2  hi 


(4.6.13) 


NR,NC,A,SIG,ND,UP,LOW,NE,NX,NDN,S,NSE,NTYPE,NDC,NCOL, 

NOT, NDTNX.NBOU.XC.YC, AREA,  -  See  Glossary. 

W  *  Working  Storage 

MUS  «  Vector  containing  stress  p's  (length  *  NSNL) 

MUD  =  Vector  containing  displacement  p's  (length  =  NDCNL) 
MUX  =  Vector  containing  redundant  p's  (length  ■  NXNL) 

IRS  =  Vector  identifying  the  active  stress  constraints 
and  row  label. 

IRD  *  Vector  identifying  the  active  displacement  con¬ 
straints  and  row  label. 

IRX  *  Vector  identifying  the  active  redundants  and  row 
label. 

IRA  «  Vector  identifying  the  unknown  design  variables 

and  column  label.  , 

YI  •»  Mise8-Hencky  criteria  for  the  element. 

Each  partial  is  summed  into  the  ARAY  analysis  matrix 
by  S4310A. 

None 

CALL  S4621  (ARAY,NR,NC,MUS,MUD,MUX,A,SIG,YI,ND,UP, 

LOW , N3 , NX, NDN , S , NSE , NTYPE , W , NDC , IRS , IRD , IRX , IRA , NCOL , 

NDT , NDTNX , NBOU , XC , YC , AREA) 

113 

None 


Scratch  Tapes: 


None 


Storage  Required:  (2312  bytes)  578  words 

Subroutine  User :  S461 

Subroutine  Required:  S466,  S4310A 


1.  Subroutine  Name: 


2.  Purpose: 

3.  Equations  and 

Procedures : 

4.  Input  Arguaents: 


5.  Output  Arguments: 

6 .  Error  Returns : 

7.  Calling  Sequence: 

8.  Input  Tapes: 

9 .  Output  Tapes : 

10.  Scratch  Tapes: 

11 .  Storage  Required : 

12 .  Subroutine  User : 


13.  Subroutine  Required 


S4622 


Compute 


„  1  _ 3  V  i  .  i  .  i  ,  _ 

9Xj8XK  3Sl<S->3S1C“>  1  «  1  ‘I  ' 

A,SIG,Y,S,NTYPE,AREA,XC,YC,NE,NSE,NDT,NX, IR, IC, 
NTRANS , NR , NC , NCOL , NDTNX , NROW , NBOU 

W  -  Working  storage 

MUS  -  Vector  containing  stress  y’s  (length  -  NSNL) 
IRS  -  Vector  identifying  this  active  u*s. 

Partial  after  multiplication  by  pc1  that  corresponds 
to  it  and  is  active  is  summed  into  the  ARAY  analysis 
matrix  by  S4310A. 

None 

CALL  54622  (A,SIG,Y,S,NTYPE,AREA,XC,YC,NE,NSE,NDT, 
NX, IR, IC, NTRANS , ARAY, NR.NC, NCOL, NDTNX, NR0W,NB0U,W, 
MUS, IRS) 

NTAPE-I12. 

None 

None 

(3268  bytes)  817  words 
S461 

S4310A, S4622A,TRPRNT 
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1. 

Subroutine  Name: 

S4622A 

2. 

Purpose: 

2 

To  compute  components  of  3  Yi  for  triangles 

and  quads.  (m) 

3. 

Equations  and 
Procedures: 

The  bar  and  shear  elements  have  a  zero  value  for 
this  quantity. 

The  quantities  output  from  this  routine  are  used 
to  form  a  3  x  3  matrix  in  S4622  for  the  triangular 
element. 

The  quad  element  uses  AA  to  FF  with  some  modification 
to  form  a  5  x  5  matrix  in  S4622. 

4. 

Input  Arguments: 

S1,S2,S3,A,Y  -  Mises  Hinckey  criteria. 

\rea 

5. 

Output  Arguments: 

AA,BB,CC,DD,EE,FF  -  components  of  stress 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  S4622A  (S1,S2,S3,A,Y,AA,BB,CC,DD,EE,FF) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(672  bytes)  168  words 

12. 

Subroutine  User: 

S4622 

13. 

Subroutine  Required: 

None 
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1.  Subroutine  Nome: 


S466 


2 .  Purpose : 


Calculate  2nd  partial  of  go  with  respect  to  A.  and 

V  3 


3.  Equations  and 
Procedures: 


±JSl 

aA.j9xK 


*  2  * 

Aj  at 


bl  *k  6i^ 


(4.6.6) 


4.  Input  Arguments:  A,SIG,Y,S,NTYPE,AREA,XC,YC,NE,NSE,NDT,NX,IR,IC, 

NTRANS , NR , NC , NCOL , NDTNX ,  NROW , NBOU ,  -  See  Glossary. 

W  -  Working  Storage 

MUS  -  Vector  containing  stress  y's  (length  -  NSNL) 
IRS  -  Vector  identifying  the  active  y's. 

5.  Output  Arguments:  The  element  of  the  partial  multiplied  by  the  corres¬ 

ponding  yo  is  put  into  ARAY  the  analysis  matrix. 

6.  Error  Returns:  None 

7.  Calling  Sequence:  CALL  S466  (A,SIG,Y,S,NTYPE,AREA,XC,YC,NE,NSE,NDT, 

NX , IR , 1C , NTRANS , ARAY , NR , NC , NCOL , NDTNX , NROW , NBOU , W , 
MUS, IRS) 

8.  Input  Tapes:  112 

9.  Output  Tapes:  None 

10.  Scratch  Tapes:  None 

11.  Storage  Required:  (2698  bytes)  677  words 

12.  Subroutine  User:  S4621,S4621A 

13.  Subroutine  Required:  S4310A.TRPRNT 
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1.  Subroutine  Name: 


TAP Ell 


2.  Purpose: 


Generate  files  Ill, and  structure  cutter  Initialization 


3.  Equations  and  File  NSS1  (unit  1)  is  read  and  the  structure  cutter 

Procedures:  matrix  is  generated.  This  matrix  is  processed  by 

routine  AA.  File  Ill  (General  Problem  Information) 
is  written  for  further  processing  in  the  program. 
Routine  S241  is  called  to  generate  the  basic  normalized 
structural  matrices  4>  and  ft. 


4. 

Input  Arguments: 

All  calling  arguments  set  up  working  storage  for 
this  routine  and  its  subroutines. 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  TAPE11  (A,A7 ,A8,KL, IELT, SIGU,DISPU,DISPL,AREA, 

ALL, DELTA, ALL2.A1.A2, A3, A4,A5,A6,WT,XC,YC,W, BETA, VA,NR,NC) 

8. 

Input  Tapes: 

NSS1 

9. 

Output  Tapes: 

Ill, 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(5396  bytes)  1349  words 

12. 

Subroutine  User: 

LINK2 

13. 

Subroutine  Required: 

S241,MRPRNT,AA 
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1. 

Subroutine  Name: 

TESTR 

2. 

Purpose : 

Test  for  divergence  in  design. 

3. 

Equations  and 
Procedures: 

a)  Compute  the  ratio  of  (new  design/old  design). 

b)  Compare  ratio  to  a  criteria  and  set  error  switch 
if  ratio  is  greater  than  criteria  or  less  than 
1/criteria. 

4. 

Input  Arguments: 

A,G,W,IAREA,NE,NXNL,PRI  -  See  Glossary  for  defi¬ 
nitions. 

5. 

Output  Arguments: 

I  ERR 

6. 

Error  Returns: 

IERR 

7. 

Calling  Sequence: 

CALL  TESTR  (A,G,W,IAREA,NE,NXNL,IERR,PRI) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

None 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(964  bytes)  241  words 

12. 

Subroutine  User: 

S451.S461 

13.  Subroutine  Required:  None 


1. 

Subroutine  Name: 

TRPRNT 

2. 

Purpose: 

Transpose  matrix  print. 

3. 

Equations  and 
Procedures: 

Print  transpose  of  a  matrix. 

4. 

Input  Arguments: 

NR  -  Number  of  rows  in  transposed  matrix. 

NC  -  Number  of  columns  in  transposed  matrix 

AKTRIX  -  Matrix  to  print 

5. 

Output  Arguments: 

None 

6. 

Error  Returns: 

None 

7. 

Calling  Sequence: 

CALL  TRPRNT  (NR.NC.AMTRIX) 

8. 

Input  Tapes: 

None 

9. 

Output  Tapes: 

101,  Unit  10 

10. 

Scratch  Tapes: 

None 

11. 

Storage  Required: 

(490  bytes)  123  words 

12. 

Subroutine  User: 

FSD, S4622 , S241 , S231 

13. 

Subroutine  Required: 

None 
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A. 5  PROBLEM  SIZE  LIMITATIONS 


The  key  limitations  are  brought  about  by  the  total  number  of 
elements ,  number  of  load  conditions  and  the  number  of  degrees  of  free¬ 
dom  in  the  problem. 

Due  to  dynamic  storage  allocation  techniques  the  problem  size 
can  be  controlled  by  changing  two  cards  in  the  'MAIN'  routine.  Both 
the  dimension  of  the  'WORK'  array  and  the  size  of  the  variable  'NWORK' 
must  be  equal.  The  delivered  size  is  20,100  words,  but  may  be  adjusted 
to  your  system. 

If  there  is  insufficient  storage  space  defined  for  a  problem  the 
program  will  print  a  message  indicating  the  amount  of  storage  required 
for  the  problem  and  the  amount  of  storage  reserved  by  the  above  mentioned 
cards  in  the  'MAIN'  routine.  To  execute  the  problem  modify  the  dimension 
and  'NWORK'  variable  to  be  what  the  problem  needs  or  reduce  the  number  of 
elements,  number  of  load  conditions  and/or  number  of  degrees  of  freedom 
indicated  on  the  input  sheet. 
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APPENDIX  B 

RAPID  REANALYSIS  COMPUTER  CODE 

S.  Gellin 
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B . 1  INTRODUCTORY  REMARKS 


An  interactive  FORTRAN  computer  code  was  written  to  perform  static 
and  dynamic  analyses  of  damaged  structures  using  the  rapid  reanalysis 
theory,  procedures  and  terminology  discussed  in  Section  2.4  of  Volume  1, 
Ref.  1.  Structural  idealizations  are  limited  to  the  use  of  bar  elements. 
Both  large  scale  damage  (d^i=l.  Type  A)  and  small  scale  damage  (d,  i<l. 
Type  B)  cases  are  included  in  the  computer  code. 


B.2  PROGRAM  CAPABILITIES  AND  LIMITATIONS 


The  computer  program  described  herein  performs  basic  static  and 
dynamic  analyses  and  rapid  reanalyses  for  dama^..  types  A  and  B  on 
plane  truss  structures.  These  can  either  be  subjected  to  a  given 
applied  load  or  be  in  a  state  of  free  vibration,  but  not  both  simul¬ 
taneously.  Several  damage  cases  of  either  type  may  be  considered  within 
a  given  analysis  case. 

The  maximum  number  of  elements  is  limited  to  50.  The  maximum 
number  of  grid  points  is  25,  thus,  the  maximum  number  of  degrees  of 
freedom  is  50.  In  addition,  the  maximum  number  of  reaction  forces  is 
10  and  the  maximum  number  of  applied  loads  (global  X  and  Y  directions 
only)  is  10.  For  a  given  structure,  the  maximum  number  of  redundants 
(which  equals  the  number  of  elements  plus  the  number  of  reactions 
minus  the  number  of  degrees  of  freedom)  is  20.  The  minimum  number  of 
redundants  is  1;  thus,  truss  structures  that  are  statically  determinate 
cannot  be  analyzed. 

Roller  supports  are  assumed  to  act  in  one  of  the  coordinate  direc¬ 
tions.  Caution  is  advised  when  performing  a  type  A  damage  analysis,  as 
the  sophisticated  grid  point  removal  scheme  as  discussed  in  Volume  1 
is  not  included  in  the  computer  program.  The  degrees  of  freedom  are 
numbered  according  to  grid  point  number  and  direction.  Thus,  grid  point 
j  has  d.o.f  2j-l  in  the  X-direction  and  2j  in  the  Y-direction. 
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B . 3  DESCRIPTION  OF  INPUT  DATA 


The  so-called  "free  format"  is  used  to  input  data  into  the  program. 
Each  "card"  or  "line"  has  several  numbers  separated  by  commas .  A  de¬ 
tailed  description  follows.  [Bracketed  paragraphs  are  only  needed  if  a 
dynamics  case  is  being  run.] 


Croup  Number  of  Cards  Quantities  Read 

1  1  Number  of  grid  points,  number  of  elements,  number 

of  applied  loads. 

[Comment:  If  a  dynamics  analysis  is  desired,  specify  the  number  of  applied 
loads  as  zero.] 

[2  1  Number  of  desired  modes  and  frequencies,  maxi¬ 

mum  number  of  iterations  used  in  eigenvalue 
routine. ] 

3  Number  of  grid  X-coordinate  of  grid  point,  Y-coordinate  of 
points,  in  grid  point,  support  code, 

numerical  order. 


Comment:  Support  code  =  0  for  unsupported  grid  point;  “  1  for  constrained 
in  the  X-direction  only;  =  2  for  constrained  in  the  Y-direction  only;  ■  3 
for  pinned  condition. 

U  Number  of  ele-  First  grid  point  number  of  element,  second  grid 
ments,  in  point  number  of  element, 

numerical  order. 


Comment:  "First"  and  "second"  grid  point  numbers  are  in  no  particular 

order,  i.e.,  they  may  be  reversed. 

5  See  comment  First  element  number,  last  element  number,  cross 

below.  sectional  area. 


Comment:  The  form  of  this  "card"  helps  to  specify  element  properties  in 
an  efficient  manner,  particularly  if  most  of  the  elements  have  the  same 
value  for  the  particular  property.  As  an  example,  suppose,  for  a  12-bar 
truss,  that  all  the  elements  have  a  value  for  a  given  property  of  1.,  ex 
cept  element  8,  which  has  a  value  for  the  property  of  2.  A  series  of 
cards  which  would  specify  this  circumstance  are:  1,  7,  1. 

8,  8,  2. 

9,  12,  1. 

The  card  assigning  the  last  element  to  its  given  value  must  be  the  last 
in  the  list.  The  others  may  be  rearranged  in  any  order. 

6  See  comment  First  element  number,  last  element  number, 

after  5  above.  Young's  modulus. 
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Group  Number  of  Cards 


Quantities  Read 


[7  See  Comment 

after  5  above. 


First  element  number,  last  element  number, 
mass  density.] 


[Skip  group  8  and  proceed 


to  group  9.] 


8 


Number  of  applied  Degree  of  freedom  number,  value  of  applied 
loads .  load . 


Comment:  The  sign  of  the  value  determines  the  sense  of  the  load  for 
the  specified  degree  of  freedom. 


9  1  Number  of  damage  cases. 

Comments:  If  the  number  of  damage  cases  is  specified  as  zero,  then  the 
program  ends  after  completing  a  basic  static  or  dynamic  analysis.  No 
further  input  cards  would  be  needed. 


10  Number  of  damage  See  comment  below, 
cases. 


Comment:  What  follows  is  a  description  of  several  subgraipings  of  input 
data  listed  as  10.1,  10.2,  etc.  The  collection  of  these  subgroups  form 
the  group  type  10.  The  number  of  group  10  type  units  is  the  number  of 
damage  cases. 

10.1  1  Damage  type,  number  of  damaged  elements. 

Comment:  For  type  B  damage,  set  type  -  1;  for  type  A  damage,  set  type  =  2. 
If  damage  type  A  is  specified,  proceed  to  subgroup  10.3. 

10.2  Number  of  damaged  Number  of  an  element  that  has  sustained  damage, 

elements.  stiffness  damage  level. 

[Proceed  to  subgroup  10.4.] 

End  of  group  type. 

10.3  Number  of  damaged  Number  of  an  element  that  has  sustained  damage, 
elements. 

[Proceed  to  subgroup  10.4.] 

End  of  group  type. 

[10.4  Number  of  damaged  Mass  damage  level.] 
elements. 

[Comment:  Order  cards  in  reference  to  the  order  of  elements  presented 
in  subgroups  10.2  or  10.3.] 


[End  of  group  type.] 
End  of  input. 
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B.4  DESCRIPTION  OF  OUTPUT  DATA 


The  results  of  the  basic  analysis,  whether  static  or  dynamic, 
is  printed  first,  followed  by  the  results  of  the  various  damage  cases 
stipulated. 

For  static  analysis  only,  three  classes  of  output  are  given  and 
labeled  as  "ELEMENT  FORCES",  "REACTION  FORCES",  and  "DISPLACEMENTS". 

"ELEMENT  FORCES"  Number  of  elements.  Element  number,  corresponding 

force. 

Comment:  Positive  quantities  imply  tension,  negative  quantities  imply 
compression. 

"REACTION  FORCES"  Number  of  reactions.  Degree  of  freedom  number  that 

reaction  corresponds  to,  re¬ 
action  force. 

Comment:  Reactions  positive  in  negative  directions. 

"DISPLACEMENTS  Number  of  degrees  Degree  of  freedom  number, 

of  freedom.  displacement. 


For  dynamic  analyses  the  above  headings  and  outputs  are  still  used, 
but  they  are  taken  to  mean  force,  reaction  and  displacement  mode  shapes, 
respectively.  The  three  sets  are  preceded  by  the  frequency  corresponding 
to  the  mode  shape,  and  a  set  of  these  four  subsets  are  presented  for  the 
number  of  modes  specified  in  group  2  of  the  input. 

If  a  set  of  output  data  is  presented  for  a  damage  case,  the  set  is 
preceded  by  a  title  specifying  whether  it  is  for  Type  A  or  B  damage.  In 
addition,  there  are  diagnostic  messages  to  indicate  trouble  spots;  for 
instance,  if  the  convergence  criteria  for  a  given  calculation  was  not  met 
in  a  certain  prescribed  number  of  calculations. 
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B.5  ILLUSTRATIVE  EXAMPLE 


A  seventeen  bar  truss  was  used  to  demonstrate  program  capabilities. 
Figure  B.l  depicts  the  truss  and  lists  the  damage  cases.  Cases  Dl-2  are 
the  dynamic  counterparts  to  static  cases  Sl-2.  The  input  data  for  the 
statics  problems  is  shown  in  Figure  B.2,  and  the  analysis  results  in 
Figure  B.3.  The  input  for  the  dynamics  problems  is  shown  in  Figure  B.4, 
and  the  output  in  Figure  B.5. 

The  estimated  compile  time  for  this  program  is  about  6.2  seconds  on 
the  IBM  3031  computer  using  the  VSPC  system;  the  estimated  execution  time 
is  .2  seconds  per  analysis  (i.e.,  3  mode  shape  determinations  is  counted 
as  3  analyses). 

A  program  listing  is  given  in  Figure  B.6. 


Figure  B.l  Seventeen  Bar  Truss  Example 
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Figure  B.2  Input  for  Statics  Cases 
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Figure  B.3  Output  for  Statics  Cases 


FigureB. 3.1  Basic  Analysis 


ELEMENT  FORCES 


1 

0 • 3S01E+06 

2 

0.70&2E+05 

3 

-0.3499E+06 

4 

—0. 48S2E+03 

1 

5 

0.2495E+06 

] 

A 

0.7149E+05 

1 

7 

-0.2505E+06 

1 

8 

0.4A86E+04 

S 

9 

0 . 1552E+06 

4 

10 

0.A331E+05  _ 

7 

11 

-0. 1448E+06 

i 

12 

-0.4477E+05 

i 

13 

0. 1414E+0A 

10 

14 

-0. 1000E+06 

n 

15 

-0.7080E+05 

u 

16 

-0.6993F+05 

U 

u 

17 

-0 . 781 1E+05 

IS 

REACTION 

FORCES 

l» 

17 

1 

-0.4000E+06 

II 

2 

-0 . 500AET05 

>• 

3 

0.4000E+06 

20 

21 

4 

-0.4994E+05 

n 

DISPLACEMENTS 

is 

u 

1 

0.0 

» 

2 

0.0 

H 

3 

0.0 

17 

4 

0.0 

It 

5 

-0.2333E+00 

s 

6 

-0.3275E+00 

N 

7 

0, 2334E+00 

11 

8 

-0.3278E+00 

11 

9 

-0.4003E+00 

11 

10 

-0. 1057E+01 

M 

11 

0.3997E+00 

U 

12 

-0.1054E+01 

U 

13 

-0.4968E+00 

11 

14 

-0.2035E+01 

M 

15 

0.5032E+00 

M 

16 

“0 . 2064E+01 

It 

17 

-  -0.5635E+00 

It 

18 

-0.3320E+01 

"  a«  d»NiM  •*(«*(  mmoj 


Figure  B.3  (Continued) 


Figure  B.3. 2  Damage  Case  SI 


i 

i 

] 

4 

J 

t 

I 

I 

1 

I* 

II 

II 

I) 

14 

II 

I* 

II 

II 

l» 

» 

>1 

» 

II 

14 

»[ 


TYPE  B  DAMAGE 

ELEMENT 

FORCES 

1 

0.34B4E+06 

2 

0. 7291E+05 

3 

-0 . 3516E+06 

4 

-0 . 3336E+04 

5 

0.24B2E+06 

6 

0.7323E+05 

7 

-0.2318E+06 

8 

0. 1835E+04 

9 

0.1536E+06 

10 

0 . 6559E+05 

11 

-0.1464E+06 

12 

-0.4638E+05 

13 

0« 1414E+06 

14 

-0.1 000E+06 

15 

-0 . 6851E+05 

16 

-0.6819E+05 

17 

-0 . 7583E+05 

REACTION 

FORCES 

1 

-0 • 4000E+06 

2 

-0 . 4845E+05 

3 

0. 4000E+06 

4 

-0.5155E+05 

_ DISPLACEMENTS  _ 

1 

0.0 

2 

0,0 

3 

0.0 

4 

0.0 

5 

-0 . 2344E+00 

6 

-0. 3316E+00 

7 

0.2323E+00 

B 

_ -0.3338E+00 

9 

-0.4022E +00 

10 

-0.1066E+01 

11 

0.3978E+00 

12 

-0.1065E+01 

13 

-0.4998E+00 

14 

-0.2050E+01 

IS 

0.5002E+00 

16 

-0.2081E+01 

...  17 

-0.5665E+00 

1  18 

-0.3336E+01 

197 


Figure  B.3  (Continued) 


Figure  B.3.4  Damage  Case  S2 
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Figure  B.4  Input  for  Dynamics  Cases 
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Figure  B.5  Output  for  Dynamics  Case 


Figure  B.5  (Continued) 
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Figure  B.5 


Figure  B.5. 2 
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Figure  B.5. 2  (Continued) 
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Figure  B.5  (Continued) 


Figure  B.5. 4  Damage  D2 
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Figure  B.5.4  (Continued) 
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